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1 Introduction 



1.1 The motivation 

The Deep Inelastic Scattering (DIS) of leptons off hadrons is interpreted as a knockout 
of one of the charged partons of the target by hard Rutherford scattering followed by 
a complete shattering of the target nucleon or nucleus. One of major discoveries at 
the electron-proton collider HERA at DESY was the observation that the large rapidity 
gap events, in which the target nucleon emerges in the final state with a loss of a very 
small fraction of its energy-momentum, constitute a substantial and approximately scaling 
fraction of high-energy/ small-x DIS of electrons and positrons on protons [1,2]. Although 
the major features of such events and their cross sections have been correctly predicted 
within perturbative QCD [3] , the very existence of large rapidity gap events for nuclear 
targets is nearly paradoxical: as well known the deposition of dozen MeV energy is already 
sufficient to break up the target nucleus, still the theory predicts that for a sufficiently 
heavy nucleus and for the Bjorken variable x < 10~^ the fraction of rapidity-gap DIS 
with retention of the target nucleus in exactly the ground state must be exactly 50 per 
cent [4] and there is a direct evidence for that from the E665 Fermilab experiment [5] . The 
discovery of rapidity gaps at HERA has led to a renaissance of the physics of diffractive 
scattering in an entirely new domain, in which the large momentum transfer from leptons 
provides a hard scale. It also vindicated the early suggestions of Bjorken to look for hard 
diffraction in hadronic interactions [6] and stimulated a revival of the rapidity gap physics 
with hard triggers — large-p_L jets, H^^-bosons, excitation of heavy flavors, — at the 
proton-antiproton collider Tevatron (for the recent review see [7-9] and references therein). 
Whether the existence of such a hard scale makes the diffractive DIS tractable within the 
perturbative QCD or not has been a subject of intense theoretical and experimental 
research during the past decade or so. A good summary of the pre-1997 status of the 
vector meson production physics is found in the monograph of Crittenden [10], the pre- 
1999 status of theoretical ideas on diffractive DIS was reviewed by Hebecker [11]), for 
the general introduction into the physics of diffractive scattering see the recent books of 
Barone and Predazzi [12] and Forshaw and Ross [13]. 

The subject of this review is a special case of diffractive DIS — the exclusive production of 
vector mesons. One disclaimer is in order: we focus on the high-energy and/or very small-x 
regime dominated by the pQCD pomeron exchange and, facing the size limitations, don't 
discuss very interesting low to moderate energy data from the HERMES collaboration 
which are strongly affected by the non- vacuum exchanges (for the review and references 
see [14]). The past decade the topic of high-energy diffraction has been dominated by new 
fundamental data coming from the ZEUS and HI experiments at HERA. The interest in 
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the exclusive electroproduction of vector mesons is multifold. From the purely experimen- 
tal point of view, the HERA experiments offer a prime example of diffractive scattering at 
energies much higher than were attainable before. Furthermore, the self-analyzing decays 
of spin-1 vector mesons allow one to unravel the mechanism of diffraction in full com- 
plexity. Specifically, the HERA experiments for the first time gave an unequivocal proof 
that the s-channel helicity non-conservation persists at highest available energies [15,16]. 
On the theoretical side, starting from the seminal papers on the color dipole approach 
by Kopeliovich, Zakharov et al. [17-21] and the related momentum space approach by 
Ryskin [22] and Brodsky et al. [23], it has been understood that the exclusive diffractive 
production of vector mesons in DIS is a genuinely hard phenomenon, whose major fea- 
tures can be described by pQCD. This can be understood in terms of the shrinkage of the 
photon with the increase of the hard scale [17, 18,24], and because of this shrinkage the 
diffractive production probes the hadronic properties of the photon and vector mesons 
at short distances. One of the direct manifestations of this shrinkage of the photon is 
a decrease of the diffraction slope with the increase of the hard scale [25,26], which has 
for the first time been observed at HERA [27,28], for the earlier evidence from the NMC 
experiment see [29]. Finally, the presence of the hard scale enables one to test the mod- 
ern theoretical ideas on the mechanism of the t-channel exchange with vacuum quantum 
numbers, i.e., the QCD Pomeron. The way the QCD Pomeron is probed in diffractive 
vector meson production is similar to, but still different from, that in the conventional 
inclusive DIS. For instance, large-t diffractive production of vector mesons probes the 
QCD Pomeron in a hard regime [30] inaccessible in inclusive DIS. 

1.2 From inclusive DIS to DVCS to exclusive vector meson pro- 
duction 

To this end recall the basics of inclusive DIS of leptons off nucleons 

e{k) p{P) e{k') X. 

To the lowest order in QED it is treated in the one-photon exchange approximation. 
The leptons serve as a source of virtual photons of energy u and virtuality = —q^ 
(the scattering kinematics and the 4- momenta are shown in Fig. 1) and the fundamental 
process is the virtual photoabsorption 

In the fully inclusive DIS only the scattered lepton is detected and one sums over all the 
hadronic final states X. Then the observed inclusive DIS cross section is proportional to 
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Figure 1: The kinematics of DIS 

the absorptive part of tlie forward, at vanisliing momentum transfer A, virtual Compton 
scattering amplitude T^u{^, Q}, Qfn^ A = 0) shown in Fig. 2 




Figure 2: The unitarity relation between DIS and forward Compton scattering 



i\Q') p{P) ^ i\Q') p{P) (1) 

and, invoking the optical theorem, can be cast in the form of the flux of virtual transverse 
(T) and scalar (longitudinal) (L) photons times the total photoabsorption cross sections 
(Tt and ai. 

Now take a closer look at the Compton scattering amplitude as a function of the virtuality 
of the incident {in) and final (/) state photons, Qf^ and respectively. In fully inclusive 
DIS this amplitude is accessible only for Qf^ = = and at vanishing momentum 
transfer A = 0. When continued analytically to = the amplitude 7^i.(z/, 0, Q^) 
will describe the exclusive real photon production often referred to as the Deeply Virtual 
Compton Scattering (DVCS) [31] 

7*{Q')p{P)^lp{P'), (2) 

while the further continuation to = —fny gives the amplitude of the exclusive vector 
meson production 

l*{Q')p{P)^V{v)p{P'). (3) 
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Both DVCS and exclusive vector meson production can be studied experimentally by 
selecting a special final state X = 7 p{P') ot X = V p{P'), respectively. Furthermore, 
both reactions can be studied at the non- vanishing momentum transfer A, i.e., t = — 7^ 
0, for the definition of the kinematical variables see Fig. 3. 



Figure 3: Schematic diagram of exclusive vector-meson production in ep inter- 
action, ep eVp. 

The point that inclusive DIS, DVCS, and exclusive vector meson production are described 
by the same analytic function taken at different values of suggests from the very outset 
the complementary probe of high energy pQCD in the three reactions (1), (2), and (3). For 
instance, in the forward, A = 0, Compton scattering probed in inclusive DIS the helicity 
flip amplitudes vanish for the kinematical reason. In contrast to that, the inclusive vector 
meson production at A 7^ enables one to determine the full set of helicity-conserving and 
helicity-flip amplitudes and investigate the spin properties of hard (generalized) Compton 
scattering to full complexity. 

1.3 When vector meson production is dominated by small color 
dipole interactions? 

The intimate relationship between inclusive DIS, DVCS, and exclusive vector meson pro- 
duction is still better seen in the lightcone color dipole picture of small-x DIS which 
illustrates nicely the interplay of the scattering mechanism and the (partonic) structure 
of particles. It is needless to recall the outstanding role of the photon-matter interactions 
in the conception and formation of the quantum mechanics and quantum field theory. In 




P(P') 
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the early years of the nonrelativistic quantum mechanics the photon has been regarded 
as structureless and the focus of the theory was on spectral lines, photo-effect and the 
related phenomena. With the advent of the first quantum field theory — the Quantum 
Electro Dynamics (QED), — it has become clear that the fundamental transition 



7 4^ e^e 



(4) 



between bare particles gives rise to a concept of a dressed physical photon that contains 
all bare states to which it couples via (4) and still higher order QED processes. At 
low energies, the virtual vacuum polarization gives rise to the well known Uehling-Serber 
radiative correction to the Coulomb potential; at higher energies the familiar Bethe-Heitler 
e^e~ pair production in the Coulomb field of a nucleus can be viewed as materialization 
of the e^e~ component of the physical photon (see Bjorken, Kogut, Soper [32]). The 
Compton scattering which is behind inclusive DIS at very small values of the Bjorken 
variable x can be viewed as (i) the transition of the virtual photon to the qq pair (the 
color dipole) at a large distance 

1 



(5) 



upstream the target (here is the nucleon mass), (ii) interaction of the color dipole with 
the target nucleon, and (iii) the projection of the scattered qq onto the virtual photon 
(Fig. 4a). Notice the very special choice of the stage (iii): if one lets the scattered color 
dipole materialize as hadrons, one ends up with the large rapidity gap DIS — the diffrac- 
tive excitation 7* p{P) p{P')X. Here the production of continuum hadronic states X 
is modeled by the continuum qq states (Fig. 4b), whereas the projection of the scattered 
qq color dipole onto the vector meson gives the exclusive (diffractive, elastic) vector me- 
son production, and projection onto the real photon gives so-called DVCS (Fig. 4c). The 






1^ 

C) 



■p 



Figure 4: The unified picture of Compton scattering, diffraction excitation of the 
photon into hadronic continuum states and into the diffractive vector meson 



amplitude of the transition of the photon into the qq state, alias the qq wave function 
of the photon, and the amplitude of scattering the color dipole off the target are the 
universal ingredients in all the processes. The wave function of the virtual photon is well 
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known [24] , and different processes probe the color dipole scattering amplitude at different 
dipole sizes [18]. 

For instance, irrespective of the photon's virtuality Q^, the inclusive diffractive DIS into 
the continuum states is controlled for the most part by interaction of large color dipoles [3]. 
The scaling violations in the proton structure function (SF), dF2p{x,Q^) / dlogQ"^ , come 
from small color dipoles of size 

2 4 

~Q2 + 4^2' (6) 

whereas the absolute value of F2p{x, Q^) receives contributions from large to small dipole 
sizes [33,34], 

^ < ^2 ^ 

Q2 + 4m2 ^2 ^ ^ 

In contrast to the inclusive DIS and inclusive diffractive DIS, the amplitude of the exclusive 
vector meson production is dominated by the contribution from small dipoles of size [18,19] 

6 

r rs 



m 



often referred to as the scanning radius (formula (8) is applicable only if rs is smaller 
than the typical hadron size). This, exclusive vector meson production offers a cleaner 
environment for testing transition from soft to hard scattering. 

The color dipole formalism is entirely equivalent to the BFKL formalism of the (trans- 
verse) momentum dependent gluon distributions in the leading log ^ approximation [35- 
37]. Within this formalism, often referred to as the A;_L-factorization, Eq. (8) suggests 
that the vector meson production probes the gluon density of the target at pQCD hard 
scale [19-22,38] 

Q'^hQ^ + ml) = \, (9) 

which is large for heavy quarkonia ( J/\E', T, ...) or for large Q^. In the hard regime of 
small scanning radius, the vector meson production amplitudes will only depend on the 
wave function of vector mesons at a vanishing quark- ant iquark separation in the two- 
dimensional transverse, or impact-parameter, space. There still remains a certain sensi- 
tivity to the separation of quarks in the longitudinal direction, which nonrelativistically 
is conjugate to the longitudinal Fermi motion of the quark and antiquark in the vector 
meson or the partition of the longitudinal momentum of the vector meson between the 
quark and antiquark in the relativistic lightcone language. As a result, the vector meson 
production amplitude is not calculable from the first principles of pQCD, still the sen- 
sitivity to the soft input can to a large extent be constrained by the decay V — ^ e~^e~, 
which proceeds via the short- distance annihilation qq —>■ e^e~ . Then Eq. (8) suggests 
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that, upon factoring out the emerging V e+e~ decay amphtude, the vector meson pro- 
duction amphtudes will depend on the hard scale Q in a universal manner. Finally, the 
energy dependence of the vector meson production amplitude offers a more local probe 
of the properties of the hard pQCD Pomeron than the inclusive DIS. 

1.4 The scale for the onset of hard regime 

Before opening the issue of hard production of vector mesons, one needs to define the 
typical soft production. Here a brief comment on the venerable Vector Dominance Model 
(VDM) is in order. Because of the obvious dominance by the vector meson pole contribu- 
tion, the point that at = —friy the amplitude of the production of the timelike virtual 
photon 7*((5j) will be proportional to the appropriate vector meson production amplitude 
times the '-)*{— m\r)V transition amplitude is a tautology. Experimentally, the timelike 
photons are produced in the e^e~ annihilation and the '-)*{— m\r)V transition amplitudes 
are measured at the e^e~ colliders and, of course, in the decay Vq e^e~ . The assump- 
tion that the ground state vector meson pole contribution dominates the photoproduction 
amplitudes, and the 7*(Qy.)l^ transition amplitude does not vary substantially from the 
vector meson pole = — down to = is the basis of the very successful VDM 
as formulated by Sakurai [39], Gell-Mann, Zachariasen, Scharp and Wagner [40,41] (for 
the comprehensive review of foundations and tests of the VDM, see Bauer et al. [42]). 

From the color dipole point of view, the success of the VDM in real photoproduction 
derives from the proximity of the distribution of color dipoles g/g/ in the ground state 
vector mesons and in the real photon. So, the qjqf component of the physical photon can 
be approximated by the corresponding vector meson (quarkonium) and the amplitude of 
interaction of the color dipole with the nucleon can be approximated by the vector meson- 
nucleon scattering amplitude, for an illustration see Fig. 5. From the naive quark model 



Figure 5: The VDM amplitude for the vector meson photoproduction 

viewpoint, the p-meson is the hyperfine partner of the vr-meson and the 2-dimensional 
charge radius, R.,^, of the vr^ sets the relevant scale. 
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One comes to the same conclusion from the experimental! observation that to a very good 
accuracy the t-dependence of elastic vrA^ scattering, real Compton scattering 7p — > 7p 
and real photoptoduction 7p — > pp is the same [42,43]. Indeed, within the VDM which 
is a very good approximation for real photons, the differential cross sections of the latter 
two processes are proportional to the differential cross section of elastic pN sacttering. 
Then the equal t-dependence of the ttN and pN elastic scattering entails an equality of t 
he radii of he p-meson. 

Experimentally, the charge form factor of the pion is well described by the VDM p-pole 
formula and [44,45] 

2 

(r2)V2 ~ 0.55 fm ^ — . (10) 
nip 

The onset of the hard regime in diffractive vector meson production requires that the 
scanning radius is smaller than any other hadronic scale. First place, one needs r| -C 

» ^mj ^ 1 (GeV)^ (11) 

The corrections in the small parameter rs/ (ry))^^'^ depend on the wave function (WF) 
of the vector meson. With the soft, Gaussian, wave function, ipvi^) ~ exp(— r^/(ry)), 
in order for the WF-dependent corrections not to exceed < (20 -r- 30)% one needs > 
(2 4) (GeV)^. For the hard. Coulomb-like, wave functions, ^vi''") ~ exp(— r/ (r^)^/^), a 
still higher > lO(GeV)^ is needed for a similar insensitivity to the shape of the wave 
function, for the related discussion of the onset of pQCD see [46]. Even for the heavy 
J/\l/ the scanning radius at g^ = is large, 

r, ^ — ~ (r}/^)i/2 ^ 0.4 fm (12) 

(for the charmonium parameters see [47-49]), so that for the onset of the short-distance 
regime insensitive to the shape of the wave function of the J/\E' one needs > m'^j/q,- 
In the realistic QCD there is still another scale - the propagation radius for perturbative 
gluons which is small, Rc ~ (0.2 0.3)fm (for the lattice QCD evaluations of Rc see [50], 
for the origin of Rc in the instanton models of QCD vacuum see [51], the analysis of 
heavy quarkonia decays is found in [52]). The color dipole cross section is of true pQCD 
origin only for dipoles r < Rc, i.e., the fully pQCD description of diffractive vector mesons 
requires rs <Rc, i-e., 

3fi 

Q' ^ QIqcd = ^^{20^ 30)GeVl (13) 

One must not be discouraged, though: the r-dependence of the dipole cross section 
does not change any dramatically from the pure pQCD domain of r < Rc to the non- 
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perturbative domain of r > R^., and the fundamental concept of the scanning radius 
remains viable up to < 1 fm, see also the discussion in Sect. 3.3.2. 

The large momentum transfer, \t\ ^ 1 GeV^, is still another way to probe the structure of 
the photon and vector meson at short distances, r ~ ^ 1 GeV~^. It is generally 

2 2 

believed [46,53] that \t\ supersedes Q as a hard scale if |t| > Q • The caveats of t as 
large scale and of single BFKL pomeron exchange dominance will will be discussed to 
more detail in Sect. 4.11. 

1.5 The structure of the review 

In this review we focus on the onset of hard pQCD regime in exclusive vector meson 
production at HERA. The presentation of the experimental data and of theoretical ideas 
goes in parallel, and an intimate connection between the vector meson production and 
the inclusive DIS will be repeatedly underlined. For this reason the presentation of the 
theoretical ideas on vector meson production will be heavily biased towards the color 
dipole picture and its momentum-space counterpart — the so-called /c^-factorization. 

The brief description of the HI and ZEUS detectors, the kinematics of DIS and of the 
vector meson production, the event selection, the definition of major observables and of 
the spin density matrix of virtual photons is presented in Section 2. The subject of Section 
3 is an overview of basic theoretical ideas on the vector meson production. Here we discuss 
briefly the Regge theory of the soft photon and hadron interactions, the QCD approach 
to the vacuum exchange (the Pomeron), the flavor dependence, the connection between 
the vector meson production and the leptonic decay of vector mesons, the origin of s- 
channel helicity non-conservation (SCHNC) and the exclusive-inclusive duality connection 
between inclusive diffractive DIS and vector meson production. We also introduce the 
color dipole approach to DIS and vector meson production and explain how the shrinkage 
of virtual photons makes the vector meson production pQCD tractable. The unified 
microscopic QCD approach to small-x DIS and diffractive vector meson production — the 
/c^-factorization approach, which is equivalent to the color dipole approach, — is presented 
to more detail in Section 4. Here we discuss both the small-t production within the 
diffraction cone and major ideas on large-t proton dissociative reaction. This section can 
be skipped in the first reading, but is essential for understanding the status of theoretical 
calculations of the vector meson production. 

In Section 5 we start the presentation of the physics results with the helicity structure of 
the vector meson production. This includes the definition of the spin observables, an in- 
troduction into the important subject of the s-channel helicity non-conservation (SCHNC) 
and the comparison of the experimental data on the spin density matrix of produced vector 
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mesons with the theoretical expectations from the color-dipole/Zc^-factorization approach 
^ In Section 6 we discuss the Q^-dependence of the vector meson production cross sec- 
tions as well as the longitudinal-to-transverse cross sections ratios Ry = (Jl/cjt- We put 
special emphasis on the flavor dependence of cross sections, emphasize an importance of 
(Q^ -|-my) as the hard pQCD scale and comment on the sensitivity of Ry = o'l/o't to the 
short distance wave function of vector mesons. In Section 7 we review the experimental 
data on the energy dependence of the cross sections and its theoretical interpretation in 
terms of the Pomeron exchange. We show how the the change of the energy dependence 
from light to heavy flavors and from photoproduction to DIS is controlled by [Q^ + m^) 
as the hard pQCD scale. We comment on tricky points in comparison of hard scales and 
energy dependence in inclusive DIS and diffractive vector meson production and on the 
impact of interplay of the scanning radius rg and the position of the node of the radial 
wave function for the \1/(2S') production cross section. The focus of Section 8 is on the 
t-dependence of the cross sections, both in the low-t and high-t regimes. The discussion of 
low-t data centers on the Q^, fny and ly-dependence of the slope of diffraction cone. The 
recurrent theme is a universality of diffraction slopes as a function of the scanning radius 
and/or (Q^ + my) as the hard pQCD scale. The properties of the Pomeron trajectory 
aipit) extracted from the vector meson production data are discussed in detail: the exper- 
imentally observed shrinkage of the diffraction cone for the J /'^ production gives a strong 
evidence for ajpit) which decreases with t approximately linearly at < 1 GeV^, but 
then starts rising up to a]p{t) ~ 1.3 in the hard regime of large Finally, in Section 9, 
we summarize the principal findings from HERA experiments on diffractive vector meson 
production and list open issues in the pQCD interpretation of these data. 

It is important to mention here that not all currently available HERA data are always 
shown in each plot where they may belong to. This is because sometimes the published 
plots from HI, ZEUS and other authors are used without any modifications. This is 
especially valid for the "preliminary" HI and ZEUS plots, that have been shown to the 
conferences and are not yet submitted in form of "official" papers. Such plots are just 
taken as they are. If for some compilation and figures only very recent data are used it is 
explained in the correspondent caption. 



Throughout this review, the numerical results shown for the /s^-factorization are either taken from the 
PhD thesis [54] or performed specially for this review [55]. 
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2 The experimental overview 



2.1 HERA 

HERA(Hadron Electron Ring Anlage) is the world's first lepton-proton collider located 
at Deutsches Elektronen Synchrotron (DESY) site in Hamburg, Germany (see e.g. [56] 
and references therein). The HERA ring has a circumference of about 6.3 km with two 
separated synchrotron rings for electrons (positrons) and protons. It runs 10-30 m below 
ground level and has four experimental areas. In two of them the beams are made to 
collide to provide ep interactions for the experiments HI and ZEUS. The remaining two 
areas are used by the fixed target experiments: HERMES [57], which scatters longitu- 
dinally polarized electrons off stationary polarized targets, and HERA-B [58,59], which 
investigated CP-violation in the B^B^ system by scattering beam halo protons off wire 
targets (was shut down in 2000). HERA was commissioned in 1991 with the first ep 
collision observed by HI and ZEUS in the spring 1992. A major HERA upgrade took 
place during 2000-2002 break. A significant luminosity increase should be achieved by 
stronger focusing of both the electron and the proton beams, see Tab. 1 where the design 
and achieved HERA values, as well the values of HERA after upgrade, are summarized. 
Further information about HERA luminosity upgrade can be found in [60,61]. 



HERA Parameters 


Design 


2000 


Design after upgrade 


p/e beam energy (GeV) 


820/30 


920/27.5 


920/30 


p/e beam current (mA) 


160/58 


>100/>50 


140/58 


Number of bunches proton/electron 


210 


180/189 


180/189 


Time between crossings (ns) 


96 


Proton /3-function x/y (m) 


10/1 


7/0.5 


2.45/0.18 


Electron /3-function x/y (m) 


2/0.7 


1./0.7 


0.63/0.26 


Specific luminosity [c'm~'^ s'^mA'"^) 


3.4x10^9 


8x1029 


1.6x10^0 


Luminosity [cmr'^s^^) 


1.5x10^1 


2x10^1 


7x10^1 



Table 1: HERA parameters. 



2.2 The detectors HI and ZEUS 

The HI and ZEUS are general purpose detectors with nearly hermetic calorimetric cov- 
erage and a large forward-backward asymmetry to accommodate the boost of the ep 
center-of-mass in the direction of the proton beam. 
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Figure 6: The HI detector 
The main components are: 

1,7,11 - Beam, Compensating and Muon toroid magnets ; 

2 - Central tracking detector ; 

3 - Forward tracking and Transition radiators; 
4,5 - Liquid Argon Calorimeter; 

6- Superconducting coil; 
9 - Muon chambers ; 

12,13 - Warm electromagnetic and Plug calorimeters. 



12 



Ovefvis'x of the ZEUS Detector 
( hngitu(Jinol cut ) 




— I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 h 

W fn -5 m 

Figure 7: The ZEUS detector 
The main components are: 

VXD - Vertex Detector, after 2000 upgrade Silicon Microvertex Detector ; 

CTD - Central Tracking Detector ; 

FDET - Forward Detector ; 

RTD - Rear Tracking Detector ; 

F/RMUON - Forward/Rear Muon Chambers; 

BMUOI/0 - Barrel Muon Inner/Outer Chambers; 

F/B/RCAL - Forward/ Barrel/ Real Calorimeters; 

BAC - Backing Calorimeter. 
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The HI and ZEUS detectors are described in details elsewhere [62,63]. The detectors 
are shown in Figs. 6, 7. The main difference between the HI and ZEUS detectors is the 
choice of the calorimetry. In the HI case the main liquid argon calorimeter with different 
tracking detectors inside is surrounded by a large diameter superconducting solenoid thus 
minimizing the amount of inactive material in the path of the particles between the in- 
teraction point and the calorimeter. In the ZEUS case only tracking chambers are placed 
inside a superconducting solenoidal magnet, surrounded by a uranium-scintillator sam- 
pling calorimeter with equal response to the electromagnetic and hadronic components. 
Both detectors are surrounded by muon chambers. Some of the components most relevant 
for the vector meson analysis are outlined below. 

2.2.1 Tracking detectors 

Charged particles are measured both for HI and ZEUS by the central tracking detectors 
operating in magnetic field of 1.15 T and 1.43 T respectively. Both trackers are build 
mainly of drift, jet and proportional chambers. The part closest to the beam pipe in HI 
case uses silicon detectors (Central Silicon Tracker). During 2000-2001 shutdown ZEUS 
has also installed a Silicon Micro Vertex Detector that should significantly improve the 
resolution of the tracking system and the vertex reconstruction. 

The polar angle coverage is 15° < 6 < 164(165)° for Hl(ZEUS) correspondingly. The 
relative transverse-momentum resolution is (j{pt)/pt ~ 0.006pT with pt in GeV for both 
experiments. Charged particles in the forward direction are detected in the forward 
tracking detector covering the polar angle range 7° < 9 < 25° and 7° < 6 < 28° for 
the HI and ZEUS respectively, the backward part (172° < 9 < 176°) is covered by the 
backward silicon tracker, BST in the HI, and by the Small Rear Tracking Detector in the 
ZEUS cases. 

Charged tracks measured by the tracking system are used to reconstruct the interaction 
vertex for each event. 

2.2.2 Calorimetry 

The tracking detectors of HI are surrounded by a liquid argon calorimeter (LAr, 4° < 
6 < 154°, a/E : 0.12/v^ and 0.50/v^ for electromagnetic and hadronic showers corre- 
spondingly, E in GeV) and a scintillating fiber calorimeter (spaghetti calorimeter, SpaCal, 
153° < 6 < 178°, 0.075/^/E for electromagnetic showers). 

The central tracking detector of ZEUS is placed inside of a thin super-conducting coil. 
Surrounding the solenoid is the high-resolution uranium-scintillator calorimeter (CAL) 
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which covers the angular range 2.6° < 6 < 176.2° with equal response to the electromag- 
netic and hadronic components and with energy resolution of 0.18/a/E and 0.35/ \/E for 
the electromagnetic and hadronic components correspondingly. 

In 1998-2000 a Forward Plug Calorimeter (FPC, lead-scintillator sandwich calorime- 
ter [64]) was installed in the 20x20 cm^ beam hole of the forward part of the CAL with 
only a small hole of radius 3.15 cm in the center to accommodate the beam pipe. It 
extended the pseudorapidity coverage of the forward calorimeter from t] < 4.0 to rj < 5.0. 
A similar device — Beam Pipe Calorimeter (BPC) — was installed in the rear region 
of the ZEUS detector, 294 cm away from the nominal ep interaction point, mainly to 
measure the lepton scattered at very small angle. Both these calorimeters were removed 
during 2000-2001 shutdown because of the changed beam-pipe geometry for the HERA 
luminosity upgrade. 

2.2.3 Muon detectors 

The HI muon system consists of an instrumented iron return yoke (Central Muon Detec- 
tor, CMD, 4° < ^ < 171°) and a Forward Muon Detector (FMD, 3° < ^ < 17°). 

The ZEUS muon system covers the polar angles between 10° < 9 < 171°, in additional 
the forward part has additional drift chambers for high-momentum muon reconstruction 
for polar angles between 6° and 30°. 

2.2.4 Forward detectors and proton taggers 

Both HI and ZEUS have very forward detectors, placed along the beam-line in the direc- 
tion of the proton beam, 20-90 m away from the nominal interaction point. 

Forward Proton Spectrometer (EPS in HI) and Leading Proton Spectrometer (LPS in 
ZEUS) consist of movable " Roman Pots" forming together with the magnets of HERA a 
kind of magnetic spectrometer. Scattered protons with a different energy and/or angle 
compared to the nominal beam protons are separated from the beam and are detected at 
appropriate positions. 

Also in the proton direction experiments have placed simple scintillation counters (5.15 m 
and 23-24 m from the nominal interaction point in the ZEUS case, five stations between 9 
and 92 m in the HI case) that are used as Proton Remnant Taggers (PRT). These taggers 
cover very high region of pseudorapidity (e.g. 4.3 < rj < 5.8 for ZEUS) and are used to 
tag the events where the proton dissociate. 
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2.2.5 Luminosity detectors and electron taggers 

The luminosity is determined from the rate of the Bethe-Heitler bremsstrahlung process 
ep e'-fp, where the high-energy photon is detected in a lead-scintillator calorimeter 
(LUMI) located at Z = —107 m in the HERA tunnel in the ZEUS case or by a crystal 
Cherenkov calorimeter (PD) located at Z = —103 m in the HI case. 

In the lepton direction the experiments have Photoproduction Taggers (PT) at 8 and 44 
m from the nominal interaction point for ZEUS, Electron Taggers (ET) in the HI case. 
They detect leptons scattered under very small angle (less than few mrads) . The leptons 
measured in the PT (ET) are used to tag photoproduction events, thus significantly 
reducing the background. 



2.3 Kinematics and cross sections 
2.3.1 Kinematics of DIS 

Because of the small electromagnetic coupling ~ 1/137, the deep inelastic scattering 
of leptons off protons is treated in the one-photon exchange approximation. The generic 
diagram for DIS e{k) p{P) ^(A;') X is shown in Fig. 1. The relevant kinematic 

variables are: 

• = —q^ = —{k — k'Y , the negative squared four-momentum of the virtual photon; 

• W"^ = {q + P)"^ = 2mpU + — Q"^, the squared center-of-mass energy of the photon- 
proton system; 

• y = {P ■ q)/{P ■ k), the fraction of the positron energy transferred to the photon in the 
proton rest frame. 

• X = Q^/2{P ■ q), the Bjorken variable, which in the parton model interpretation of 
DIS has a meaning of the fraction of the proton's lightcone momentum carried by the 
struck charged parton. 



2.3.2 The flux and polarization of photons 

The amplitude of DIS equals 



T{e{k) p{P) e{k') X) = ^^{e{k')\lMk)) ■ 9,. ■ {X\JMP)) , (14) 



where and stand for the electromagnetic current of leptons and hadrons. The leptons 
serve as a source of photons and the physical process is the virtual photoabsorption 



^*{q)p{P) ^ X, 
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The virtual photons have three polarization states: the two spacelike transverse ones with 
helicities A^* = ±1, 

e(±) = --^(±e. + ze,)-e±^* (15) 

and the timelike scalar state (often misnamed the longitudinal one, hereafter we follow 
this tradition) 



P 



{P ■ qY + P2g2 

For the purpose of future convenience, here we choose the 2;-axis along the photon's 3- 
momentum, the x-axis in the 7p Vp reaction plane, and $ is the azimuthal angle 
between the reaction and the (e, e') scattering planes (for more details see below section 
5.1). The complete set includes still another spacelike vector 

Making use of the expansion 

g,. = e;(L)e.(L) - e;(+)e,(+) - e;(-)e,(-) - el{S)e,{S) , 
and of the current conservation, {e{S) • J) = 0, one can write down 

{e{k')\l,\e{k)) ■ g,, ■ (X| J.b(P)) = (e(A;')|/;.|e(A;))e;(L) ■ (X| J.b(P))e.(L) 

- (e(A;')|/,|e(A;))e;(+).(X|J.b(P))e.(+) 

- (e(A;')|/,|e(A;))e;(-)-(X|J.|p(P))e.(-). (16) 

Now notice that 

T(7*(A^,;g)p(P) ^ X) = v^i^ ■ (X| J,b(P))e,(A^,) (17) 

is precisely an amplitude of the photoabsorption for the photon of polarization A^^,, 
and (e(A;')|/^|e(/c))e* (L), — (e(A;')|/^|e(A;))e* (+), — (e(A;')|/^|e(A;))e* (— ) define the emission 
by leptons of photons of appropriate polarization, which is quantified by the spin density 
matrix of the photon p\i\. Then, making use of the expansion (16) the differential cross 
section for the leptoproduction of the specific final state X can be expressed through the 
photoabsorption cross sections as 

"^""^dy/^^^ ^^^^^^^^''^^ ^ PA'ArfavA(7*P - ^) , (18) 
^ ^ A',A=+,-,L 

where drx is the element of the appropriate phase space, 

rfaA'A(7> ^ X)= ^ =r-(7^P ^ X)r(7> ^ X)dTx (19) 
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and 



rTiQ\y) = ^-(i-y + ly') (20) 

is the flux of transverse photons. With this normalization the spin density matrix of the 
photon equals 







P+l\ 




( \ 












P-L 






1 






\PL+ 


PL- 


PllJ 








e 


1 



where 

^ ^ 2(1 -y) 

is the longitudinal polarization of the virtual photon. We also indicated the small-x 
approximation for e which is appropriate for DIS at HERA. 

Notice, that because of the current conservation one can define the longitudinal photon 
interaction amplitude in terms of the current component J^-, which is customary in elec- 
tronuclear physics, for instance, see [65-67]. It does not affect the observed cross section 
(18) because the different normalization of the amplitude (16) for longitudinal photons is 
compensated for by the change of the relevant components of the spin density matrix of 
the photon. 



2.3.3 The transverse and longitudinal cross sections for DIS 

In the fully inclusive DIS one integrates over the whole phase space of the state X and sums 
over all states X. Then by virtue of the optical theorem one can relate the photoabsorption 
cross section to the absorptive part of the Compton forward scattering amplitude 



y^(Jx'x{l*P ^ X)= =lmTyx{l*P l*p). (23) 



The crucial point is that for the unpolarized target the helicity-fiip. A' ^ A, amplitudes 
vanish in the forward scattering. Then the virtual photon-proton cross section, a^*^, can 
be determined from the measured positron-proton cross section: 
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where a^^ = o"++ = a and ^ are the transverse and the longitudinal virtual photo- 
production cross sections, respectively. The often discussed total inclusive cross section, 
'^tot = ^ + (Tl ^ can be determined from cr'^*^ through the relation: 

. 1±Rdi^ * . . 

l + eRois ^ ' 



where 



7 P 
L 



Rdis = ^ • (26) 



(Because R is heavily used for different ratios, we supply it by the subscript DIS.) In the 
kinematic range of most of the discussed measurements, the value of e is close to unity, 
and because Rdis is small, a"'*^ differs from a^^^ by less than one percent. 



2.4 Kinematics of diffractive vector meson production 

Diffractive vector meson production corresponds to the special two-body final state which 
contains only the vector meson and scattered proton 

e{k) p{P) e{k') V{v) p{P'), 

where \^={ci;,p,0,J/^/',\l/',T} and k, k', P, P', and v are the four-momenta of the incident 
lepton (positron or electron), scattered lepton, incident proton, scattered proton and 
vector meson, respectively, see Fig. 3. 

The new kinematic variable is t = (P — P')^ = — qY = — + tmm; the squared four- 
momentum-transfer at the proton vertex. At high energies the longitudinal momentum 
transfer = mp(Q^+my) /W"^ is small, tmin = can be neglected, and t ^ t' = — A^. 
Besides t, the new important variables are the orientation of the production plane with 
respect to the electron scattering plane and the appropriately defined polar and azimuthal 
angles of the decay pions, which will be discussed in Section 5. 

The major background process is the proton-dissociative reaction ep — eVY, and in 
addition to the above quantities. My, the mass of the diffractive excitation of the proton, 
is used. 



2.5 The event reconstruction 

For the photoproduction events, ~ 0, ranged from the kinematic minimum, = 
Mly'^/{1 - y) ~ 10"^^ GeV^, where Me is the positron mass, up to Ql^ax ~ 1 GeV^, the 
value at which the scattered positron starts to be observed in the calorimeter, with a 
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median of approximately 5 ■ 10^^ GeV^ (differs slightly for ZEUS and HI and from 
year to year with modifications in calorimeter geometry). Since the typical is small, 
it can be neglected in the reconstruction of the other kinematic variables. 

For the DIS events the kinematic variables are reconstructed using the momenta of the 
decay particles and the polar and azimuthal angles of the measured scattered lepton. 
Neglecting the transverse momentum of the outgoing proton with respect to its incoming 
momentum, the energy of the scattered positron can be expressed as: 

E,, ^ [2E, - {Ev - p^)]/il - cosOe'), 

where Eg is the energy of the incident lepton, Ey and py are the energy and longitudinal 
momentum of the vector meson V, and 6^/ is the polar angle of the scattered lepton. The 
value of was calculated from: 

=2Ee'Ee{l + C0See') . 



The photon-proton center-of-mass energy, W, can be expressed as ~ 2Ep{E —pz)v + 
Q^, where Ep is the laboratory energy of the incoming proton and {E — pz)v is the 
difference between the energy and the longitudinal momentum of the vector meson. The 
fraction of the positron momentum carried by the photon is calculated from y = {E — 
Pz)v/'2E,. 

The squared four-momentum transfer at the proton vertex is given by \t\ = (pe' +Pv')x + 
{Pe' +Pv)y- 



2.6 Data samples and event selection 

The kinematic region for each particular data sample can be found in Tabs. 2, 3. The 
tables summarize all the recent data discussed in this paper, for the overview of the 
pre-1997 experimental data, see [10]. 
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VM 


Mode 


Decay Ch. 


Q\ GeV2 


W, GeV 


t, GeV2 


Year 


Lumin.,p6 1 


Exp. 


Ref. 


PHOTOPRODUCTION WITHIN THE DIFFRACTION CONE 




Elastic, P.-Diss. 




4-10-6 


50-100 


0.-0.5 


1994 


2.17 


ZEUS 


[68] 




Elastic 


TT + TT- 


10-4 


25-70 


0.073-0.45 


part of 1999 


3.0 


HI 


[69] 




Elastic 




0.05 


26-285 


0-1.2 


1996-97 


20.5 


HI 


[70] 




Elastic 




5-10-5 


20-170 


0-1.8 


1996-97 


38.0 


ZEUS 


[71] 




Elastic 




5-10-5 


20-290 


0-1.2 


1998-2000 


55.2 


ZEUS 


[71] 




Elastic 


l+l~, J/ipTT+n- 


10-4 


40-160 


n.d. 


1993-94 


6.3 


HI 


[72] 


^(25) 


Elastic, P-Diss 


1+1-, J/ipn+TT- 


0.055 


40-150 


0-5 


1996-2000 


77.0 


HI 


[73] 


^(25) 


Elastic 


e~^e~ 


5- 10-5 


50-125 


n.d. 


1999-2000 


55.3 


ZEUS 


[74] 


T 


Elastic 




0.11 


70-250 


0-1.2 


1994-97 


27.5 


HI 


[70] 


T 


Elastic 




5- 10-5 


80-160 


n.d. 


1995-97 


43.2 


ZEUS 


[75] 


PHOTOPRODUCTION AT HIGH t 




Elastic, P.-Diss. 


TT+TT-, K+K- 


7-10-6 


85-105 


0-3 


1995 


2.0 


ZEUS 


[76] 




P.-Diss. 


TT+TT-, K+K- 

/i+yU", e+e" 


7-10-6 


80-120 


0-12 


1996-97 


25.0 


ZEUS 


[77] 




P.-Diss. 




n.d. 


50-160 


1-21 


1999 


19.1 


HI 


[78] 




P.-Diss. 




0.06 


50-200 


2-30 


1996-2000 


78.0 


HI 


[79] 



Table 2: The recent HI and ZEUS photoproduction measurements discussed 
throughout this paper. An overview of older data can be found in [10]. Information 
that is not available is labeled by "n.d.". The determination of the diffraction slope 
is based on the data at \t\ < 1 GeV^. 



VM 


Mode 


Decay Ch. 


GeV2 


W, GeV 


t, GeV2 


Year 


Lumin.,p6 1 


Exp. 


Ref. 


DEEP INELASTIC SCATTERING 


P 


r .-UlSS. 




(-OO 


fin 1 fin 


n 1 
U-i.O 


1 00/1 

iyy4 


fi 

z.o 


rli 


[finl 
[oUJ 


P 


Elastic 


TT^TT 


U.zO-U.oO 


on on 


n n fi 


1 ooc; 
iyyO 


Q fi 

o.o 




[fill 
[81J 


P 


Elastic 


TT^TT 


O-OU 


QO 1 fi7 


n n fi 
U-U.o 


1 ooc; 
iyyO 


fi n 
O.U 




[fill 
[81J 


P 


Elastic 


TT^TT 


n OK n fit; 


on on 


n n fi 
U-U.o 


1 ooc; 
iyyO 


Q fi 

o.o 




[15J 


P 


Elastic 


TT^TT 


o-oU 


/in 1 on 


n n fi 
U-U.O 


1 ooc; 
iyyO 


fi n 
O.U 




[15J 


P 


Elastic 


TT^TT 


i-DU 


'in 1 /in 
oU-i4U 


n n 
U-U.o 


1 00 ofi 
iyyO-yO 


A n 

4.U 


rli 


fi fil 


rfi 
P 


Elastic 


TX^Tl 


o K fin 


/in 1 on 
4U-izU 


n Q n 

U-O.U 


1 007 

iyy ( 


fi n 

O.U 


rli 


[fiOl 

[ozj 


rfi 
P 


Elastic 


TT^TT 


o fin 


QO 1 fin 


C\ C\ P. 

U-U.o 


1 OOfi 07 

iyyO-y / 


Qfi n 
oo.U 


ZjII/UO 


[fiQl 

[ooj 


rfi 
P 


Elastic and P.-Diss. 




o fin 


c;n 1 /in 
0U-i4U 


n o 
U-z 


1 OOfi 07 

iyyO-y / 


Qfi n 
oo.U 


ZjII/UO 


[fi/il 
[o4J 


P 


Elastic 


TT^vr 


fi Rn 


/in 1 fin 
4U-ioU 


n n 
U-U.o 


onnn 
zUUU 


/I O /I 
4z.4 


rli 


[fi K\ 
[oOJ 




Elastic 




Q on 

O-ZXj 


/in 1 on 

4U-iZU 


n n fi 
U-U.o 


1 OOfi 07 

iyyO-y ( 


Q7 7 
( . ( 




ffifil 
[oOJ 


9 


Elastic 


i\ i\ 


fi on 


/10 1 Q/1 


n n fi 
U-U.o 


1 00/1 

iyy4 


fi 

z.o 


rli 


[finl 
[oUJ 


9 


Elastic 




lie; 
i-io 


/in 1 Qn 


n n c; 
U-U.O 


1 ooc; Ofi 
iyyo-yo 


Q 1 
O.i 


rli 


[fi7l 

[8 (\ 




Elastic 




2-70 


35-145 


0-0.6 


1998-2000 


66.4 


ZEUS 


[88] 




Elastic 




2-40 


50-150 


n.d. 


1995 


6.0 


ZEUS 


[81] 


J/i, 


Elastic 




2-80 


25-180 


n.d. 


1995-97 


27.3 


HI 


[89] 




Elastic 




0.15-0.8; 2-100 


30-220 


0-1 


1998-2000 


69.0; 83.0 


ZEUS 


[90] 


^(2S) 


Elastic 




1-80 


40-180 


n.d. 


1995-97 


27.3 


HI 


[89] 



Table 3: The recent HI and ZEUS electroproduction measurements discussed 
throughout this paper. Overview of older data can be found in [10]. Notation is 
the same as in Tab. 2. 



3 An overview of theoretical approaches to diffrac- 
tive scattering 



3.1 The rudiments of the Regge theory of strong interactions 

As Bjorken emphasized, the foundations of the Regge theory are as sohd as QCD itself 
[91]. Because the physics of diffractive scattering is permeated by ideas and concepts 
from the Regge theory of strong interactions, a brief introduction into this subject is 
in order. For the more rigorous treatment and for technicahties one must consult the 
textbooks [12,13,92], the review papers [93-95] and the collection of reprints [96]. 



3.1.1 The s-channel asymptotics from the i-channel exchanges: 
spin and energy dependence 

There is a deep connection between the high-energy behavior of a binary reaction ah ^ cd 
and the spin, J, of the elementary particle with mass M exchanged in the t-channel: 

Aab-.cdiW\t) = 3'^f^3M{t\ wy (27) 

dt {t-M'^Y ' ^ ' 

Although it follows in a straightforward manner from the analysis of Feynman diagrams, 
it is instructive to look at (27) from the t-channel point of view. In the crossed channel 

a c —^h d 

the total c.m.s. energy squared is t = {pa + Pe)^ = {Pa — PcY = M"^, the momentum 
transfer squared is {pa — Pi)^ = s = {pa +PbY, and the exchanged particle emerges as a 
resonance at t = in the partial wave J. The angular dependence of this contribution 
to the scattering amplitude is given entirely by the Legendre polinomial 

Kl^JW',t) = Aj{t)Pj{cos9,) = ^^^f ^ + {-iy]PA-cos9,) , (29) 



cos^, = l + ^— ^. (30) 



where (for the sake of simplicity we take = = rric = = /i) 

The so-called signature = ±1 separates the crossing-even and crossing-odd amplitudes; 
for instance, in the crossing-even vr^vr'' scattering cr^ = +1 and the contribution from the 
odd-partial waves to (29) vanishes identically. 
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The amplitude (29) depends on only through the Legendre polinomial and can readily 
be continued analytically into the high energy domain of 3> M^, /i^, \t\, which amounts 
to — cos ^ 1 and Pj{— cos9t) oc (— cos^j)"' oc (W^)'^ , i.e., we derived the asymptotics 
(27) by analytic continuation from the t-channel to s-channel scattering. 



3.1.2 The Regge trajectories 

On the one hand, the existence of high-spin resonances is an ultimate truth of the physics 
of strong interactions; on the other hand, the exchange by an elementary particles of spin 
J > 1 would conflict the fundamental Froissart bound [97] 

A{W^,t) KWHog^W . (31) 

The sole known way out of this trouble is offered by the Regge theory: one must improve 
upon the above mock-up analytic continuation going from the sum over integer (or half 
integer) partial waves to the Sommerfeld- Watson integral over the complex angular mo- 
mentum J with which the analytic continuation to — cos Of ^ 1 must be complemented 
with the appropriate deformation of the integration contour on the complex- J plane [98]. 
The key point is that the asymptotic behaviour of the s-channel amplitude will be con- 
trolled by singularities of the partial wave Aj{t) in the complex- J plane. If the singularity 
is the (Regge) pole, 

Mt) oc ^ (32) 
J - an{t) 

then one obtains precisely the amplitude of the form (27) with J = aii{t). The t-channel 
unitarity dictates [98,99] that the Regge pole must be a moving one, i.e. it must have a 
finite slope a'jiit). Experimentally, the Regge trajectory anit) for the s-channel scattering 
at i < can be extracted from the energy dependence of the differential cross sections, and 
can be linked to the resonance mass spectrum by extrapolation of the mass-dependence 
of the spin of t-channel resonances, J„ = aR^M"^). Such Chew-Frautschi plots are well 
approximated by straight lines, 

aR{Ml) ^ aR{Q) + a'^t . (33) 

For instance, for the p,u!,A2,f2 families of resonances with non- vacuum quantum numbers 
such extrapolations suggest the intercept a_R(0) ~ 0.45, in very good agreement with the 
results from the scattering experiments. To cite few examples, the p-trajectory is best 
studied in the charge-exchange 7r~p n^n, the y42-trajectory is probed in 7r~j9 rjn, 
the cu-trajectory is probed in the regeneration Ks oii the isoscalar target, the 

TT-trajectory is probed in the charge-exchange np — >■ pn, etc. For classic reviews on the 
Regge trajectories see [95], a more recent discussion of the Chew-Frautschi plots is found 
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in [100]. The high-lying Pomeron, p,uj,A2,f2 are the natural spin-parity exchanges, i.e., 
the spin J and parity P of particles lying on the corresponding Regge trajectory are 
related by P = (—1)"^. The un-natural spin-parity vr, exchnages, P = — (— 1)"', have 
much lower intercepts, a.„{0) ~ «Ai(0) ~ 0. 



3.1.3 The universality aspects of the Regge exchange 

The Ansatz (27) bears all the salient features of the realistic reggeon-exchange ampli- 
tude: 

1. The trajectory J = aR{t) is universal for all beam and target particles, it only depends 
on the t-channel quantum numbers. 

2. Dependence on the initial and final state particles has a factorized form. 

3. If one parameterizes the t-dependence of the near forward differential cross section by 
the so-called slope parameter B, 

^^eM-B\t\), (34) 

then the factorization property entails 

B{ab^cd) = Bac + Bm + Br, (35) 

where Bac and Bm come from the form factors of the a — > c and h ^ d transitions, 
and Br characterizes the exchanged reggeon. 

4. Notice that 

(iy2)afl(t) ^ (|^2)a«{o) . (|^2)a^t ^ (H/2)"«{o) . exp[-a^|t| \og{W% (36) 

what entails Gribov's growth of the slope parameter with energy, alias the shrinkage 
of the diffraction cone [99]: 

BR = 2a'^\og( — ]. (37) 



V ■50 

The slope of all the non-vacuum Regge trajectories is about the same, 

ai, ^ ^ ^ 0.9 GeV"^ , (38) 

for the recent summary see [100]. 

5. The phase of the reggeon exchange amplitude is uniquely fixed by the analytic con- 
tinuation of the signature factor r]{at,t) = at + (— 1)"«W = at — exp[—i7r{aji{t) — 1)]: 

i?ev4(iy2^t) _ [tan[f(ai?(t) - 1)], if a* = +1, ^^^^ 



ImA{W^,t) 



cot[f(a^j(t)-l)], iiat = -l. 
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3.1.4 The vacuum exchange: the Pomeron trajectory from hadronic 
scattering 



Elastic scattering is driven via unitarity by strongly absorptive inelastic mult ipro duct ion 
processes, which is nicely illustrated by the impact parameter representation — the high 
energy version of the partial wave expansion. In high energy elastic scattering the mo- 
mentum transfer, A, is the two-dimensional vector transverse to the beam momentum. 
The elastic scattering amplitude can be cast in the form of the Fourier transform 

^A{W', A) = 2zjd'b[l- S{b)] exp(-^bA), (40) 

where S{b) = exp(2i5(b)) is the S-matrix for elastic scattering at an impact parameter 
b and the angular momentum / = |p| ■ |b|. The total elastic and inelastic cross sections 
equal 

aei = j d^b\l-S{b)\^ , 
a,^ = j d''b[l-\S{b)\^] . (41) 

Strong absorption implies the predominantly imaginary scattering phase. One often uses 
the so-called profile function r(b) = 1 — S{b). The small momentum transfer expansion 
in (40) gives 



(42) 



so that the diffraction slope B is determined by the mean impact parameter squared 

1 1 Jd'bb'Tjb) 

^-2^^^ = 2- Jd-bTib) ■ ^^2) 

The extreme case is the scattering on the absorbing black disc of radius R for which 
|S'(b)| = 6{R — |b|), which is a good approximation for the scattering of nucleons off 
heavy nuclei. Then 

f^ei = = -(^tot = vri?^ (44) 

and the diffraction slope equals 

Bel = ^R' . (45) 

Such a flat, energy independent, elastic scattering must be contrasted to the two-body 
reactions with the non-vacuum exchange which constitute a tiny fraction of high energy 
inelastic collisions of hadrons and have cross sections that vanish at high energy, 

a(ab cd) oc ^ — < . (46) 
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The importance of strong absorption for high-energy hadron interactions is evident form 
the proximity of central partial waves oipp scattering to the unitarity limit, r(b)<l, [101], 
although the periphery of the nucleon is still gray, and for all the hadrons aei is still 
substantially smaller than am, see the plots in the Review of Particle Physics [102]. As 
emphasized first by Pomeranchuk, the particle-antiparticle cross section differences vanish 
at high energy, see [102] , and from the t-channel viewpoint the elastic scattering is domi- 
nated by the vacuum exchange. In 1961 Chew and Frautschi conjectured that the vacuum 
channel too can be described by the reggeon — first dubbed the Pomeranchukon, later 
shortened to Pomeron, — exchange with an appropriate spin-2, C-even, isoscalar, posi- 
tive parity resonance lying on the Pomeron trajectory (the early history of the Pomeron 
is found in [96]. 

If the Pomeron were a simple Regge pole, it would have been utterly distinct from the 
non- vacuum reggeons: 

• For all hadrons and real photons the total cross sections rise with energy and the 
phenomenological Pomeron trajectory has ajp(O) = 1 + Ajp ^ 1.1 > 1 (notice that 
from now on the A is still used for the four-momentum exchanged, but Ai^dex is used 
to define the variation of the intercept of the coorrespondent index trajectory from 
unity). Such a rise of the vacuum component of the total cross section, 

can not go forever, though. At asymptotic energies it would conflict the Froissart 
bound. Furthermore, the partial waves of elastic scattering would overshoot the uni- 
tarity bound. Indeed, in the often used exponential approximation, -^AiW"^, A) oc 
exp(— |i?A^), and neglecting the small real part of the small-angle scattering ampli- 
tude, one finds 

r(b) = ^-exp(— ^) (48) 

and with the unlimited growth of Otot one would run into r(b) > 1. The unitarity 
(absorption, multipomeron exchange,...) corrections, which must eventually tame such 
a growth of r(b) and of Otot with energy, were shown to be substantial already at 
moderate energies [103-105]. The multipomeron absorption affects substantially the 
determination of A^: the first estimate 

Ajp ~ 0.13 

with the perturbative treatment of absorption based on Gribov's reggeon field theory 
[106, 107] goes back to the 1974-75 papers by Capella, Tran Thahn Van and Kaplan 
[103,104]. Within a more realistic model for absorption, the ITEP group [105] found 
the equally good description of the hadronic cross section data with substantially 
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larger Ap ^ 0.23. If one follows the Donnachie-Landshoff suggestion [108] to ignore 
the absorption corrections altogether and stick to the simplified pole terms, then the 
Particle Data Group finds Ajp = 0.095 [109]. However, according to the 2002 edition 
of the Review of Particle Properties [102], still better fit the the experimental data is 
provided by the parameterization [110] 



which is consistent with the Froissart bound and from the Regge theory viewpoint 
corresponds to the triple-pole singularity at j = 1, i.e., = 0! 

• The shrinkage of the diffraction cone in elastic scattering suggests very small slope 
of the Pomeron trajectory a]p{t): the combined analysis of the experimental data on 
elastic pp,pp, vr^p, K^p scattering at the CERN SPS/FNAL and CERN ISR energies 
gave a'jp ^ 0.13 ± 0.025 GeV~^ ( [111], for the review see [112]). The extrapolation 
of these fits under-predicts the pp diffraction slope at the Tevatron, which call for 
a'jp ~ 0.25 GeV~^. Incidentally, the last value of a'jp has been used by theorists ever 
since 1974-75 [103,104], but it must be taken with the grain of salt: the observed growth 
of the diffraction cone can to a large extent be due to the unitarity/absorption driven 
correlation, cf. Eqs. (44) and (45), between the total cross section and the diffraction 
slope so that the Tevatron data can well be reproduced with the still smaller values 
of a'jp [113]. 

To summarize, the Donnachie-Landshoff (DL) parameterization [108] 



must only be regarded as a convenient short hand description of the local, W < 1 TeV, 
energy dependence of the vacuum component of the elastic scattering of hadrons. 

3.1.5 The diffraction slope: variations from elastic scattering to 
single to double diffraction excitation 

The variation of the diffraction slope (35) from elastic scattering to single (SD) to double 
(DD) diffraction excitation exhibits certain universal features [93,94,114]. An excellent 
guidance is provided by a comparison of elastic proton-nucleus, pA pA, to quasielastic, 
pA p'A*, scattering. The latter reaction, in which one sums over all excitations and 
breakup of the target nucleus without production of secondary particle, must be regarded 
as diffraction excitation of the target nucleus. 

The crucial point is that at a sufficiently large {p,p') momentum transfer such that the 
recoil energy exceeds the typical nuclear binding energy, which can viewed as hard scatter- 
ing, the t-distribution of scattered protons in quasielastic (nucleus-dissociative) pA —>■ p'A* 



a,ac{ab) = a,ac{ab) = Zab + 2B\og'{W/Wo), 



(49) 



a,oft{t) = l.l + 0.25GeV-2.t 



(50) 
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is the same as in elastic pp scattering, Bdiss{pA — > p'A*) ^ Bpp [115, 116]. The quasielas- 
tic pA p'A* becomes sort of a deep inelastic scattering with quasifree bound nucleons 
behaving as partons of a nucleus and quasifree pN — > p'N scattering being a counterpart 
of the Rutherford scattering of leptons off charged partons in DIS off the proton. The 
summation over breakup of a nucleus into all continuum excitations is important, for 
excitation of the specific discrete state A* of a target nucleus, pA — > p'A*, the diffraction 
slope will still be large. 

Baa* ~ B,i ^ ^R^. (51) 
Now define the ratio of differential cross sections 

RaUo{d^ss/elm = ^^^-(^^^ ^ ^^-^^) / daApA^p'A)^ ^^2) 

Elastic scattering: Ratio{diss / el){t) ^ 1 is the dominant process within the diffrac- 
tion cone, RWt\ <^ 1. However, elastic scattering dies out rapidly for RWt\ > 1, where 
quasielastic scattering takes over: Ratio{diss / el){t) ^ 1. This point is clearly illustrated 
by the experimental data [117, 118] on elastic and nucleus-dissociative p^'^C scattering 
shown in Fig. 8. Notice the diffractive dip-bump structure, familiar from optical diffrac- 
tion, in the differential cross section of pure elastic scattering. For a sufficiently hard 
scattering, \t\ = > 0.06 (GeV) , the sum of the elastic and nucleus-dissociative cross 
sections, dasc = daei + da diss is clearly dominated by the nucleus-dissociative dadiss- 

In the regime of strong absorption the integrated cross section of quasielastic or nucleus- 
dissociative scattering is small [115], 

(Td^ss{pA ^ p'A*) « a,i{pA ^ p'A) ^ iaf„t (53) 

Exactly the same considerations apply to elastic scattering and diffraction excitation of 
hadrons and real photons, a = p,7r, K,'~f on the free nucleon target, b = p. Let -Bat be the 
contribution to the diffraction slope of electric pp scattering from the Pomeron-proton- 
proton vertex, so that 

= 2Bn + Bjp. (54) 

In the single or target-dissociative (SD) reaction, pp — ^ pY, and double dissociation 
(DD), pp XY, one must distinguish the low-mass (LM = resonances, low- mass con- 
tinuum states, ...) and high-mass (HM) states X,Y. The boundary between the low- 
mass (exclusive low mass states, resonances, ...) and high-mass continuum excitations is 
Mx,Y ~ 2 GeV. The case of small-mass excitation is an exact counterpart of excitation 
of discrete nuclear states in pA — > p'A*. Then (51) suggests that the contribution to the 
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Figure 8: A comparison of elastic (points, the lower set of data points and the 
lower curve) and combined elastic plus nucleus-dissociative (triangles, the upper set 
of data points and the upper curve) p^'^C scattering data [117, 118]. The theoretical 
calculations are from Czyz et al. [116]. 
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diffraction siope from tlie pY transition Bpy ~ -Btv? so tliat in SD and DD into low-mass 
states 

Bsd{LM) = + Bn + Bjp^ BooiLM, LM) = B^y + B^x + Bjp ^ B,i , (55) 

in good agreement with the experimental data from in the CERN ISR and FNAL exper- 
iments [119-123]. The SD into high-mass (HM) continuum, pp ^ pY{HM), corresponds 
to the complete breakup of the target proton and the reaction can be viewed as elastic 
scattering of the beam proton on one of the constituents of the target. Consequently, the 
dependence on the size of the target proton vanishes, Bpy ~ 0, and in SD into high-mass 
states (often referred to as the triple-Pomeron region) and mixed low&high mass DD 

Bsd{.HM) ^ BdD{LM, HM) ^ B^ + Bp ^ ^B,i ^ 6 GeV"^ . (56) 

In DD pp X{HM)Y{HM) with excitation of high-mass states from both the target 
and beam Bpx ~ Bpy ^ and only the t-channel exchange Bjp contributes to diffraction 
slope. Experimentally, this component is abnormally small [122, 123] 

Bdd{HM, HM) ^ Bp ^{1-2) GeV-2 . (57) 

Finally, although in irp, Kp, pp scattering only the central partial waves are close to the 
strong absorption limit, and the ratios <7ei/<Jtot ~ (0.15 ^ 0.25) are still substantially 
smaller than ^ for the strongly absorbing nuclear target, the strong inequality (Tdis{pp 
p'Y) <^ CTeiiPP pp) holds in close similarity to (53). Typically, in pp interactions 
Rpp{diss/el) = aso/o'ei ^0.3, for the review see [93,94,112]. 



3.2 The Regge theory and QCD 

In the realm of DIS the high energy limit amounts to the small-x limit. The SF's of 
small-x DIS are related to the total cross sections as 

i^T,L(x,g^) = -^aT,L(x,Q^). 

Instead of the transverse SF one usually discusses F2(x, Q^) = Ft{x,Q'^) + Fl{x,Q'^). 
The QCD parton model decomposition of the proton SF into the valence and sea quark 
contributions 

F2{x,Q') = a;5^eJ[g/(a;,Q^) + g/(x,g2)] 
/ 

= ^x-u„(x,g2) + la;.rf^(a;^g2) + 2a;^eJg/(x,Q2) (58) 
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must be viewed as a decomposition of the photoabsorption cross section into the non- 
vacuum (non-single) and vacuum (singlet) components. From the viewpoint of the QCD 
evolution, the valence component corresponds to slowing down of the valence quarks to 
a; ^ 1 and depends on the target. At small x, the sea evolves from glue and will be the 
same for the proton and neutron as well as antinucleon targets, i.e. it must be associated 
with the Pomeron exchange. The density of small-x gluons exceeds greatly the density of 
charged partons, which entails that (i) one can model high energy inelastic interactions by 
production of the multigluon final states and (ii) to the so-called leading-log ^ the small-x 
evolution is driven by the splitting of gluons into gluons, with the splitting g ^ qq only 
at the last stage of the evolution. As a result, the QCD vacuum exchange is modeled by 
the tower of color-singlet two-gluon exchange diagrams of Fig. 9, which is described in 
terms of the so-called unintegrated or differential gluon density 

O log At^ 

where n is the gluon transverse momentum. 




Figure 9: (a,h) The subset of two-gluon tower pQCD diagrams for the Pomeron 
exchange contribution (c) to the Compton scattering (DIS) and diffractive vector 
meson production. Not shown are two more diagrams with q q. 

At not so small x, the Q^-dependence of the parton densities is governed by the DGLAP 
evolution [124-126]. Here the evolution goes from smaller to larger Q^, so that once the 
boundary condition is taken at a sufficiently large Qo then one stays in the perturbative 
domain. However, in the language of inelastic multiparticle states the DGLAP evolution 
amounts to summing only the final states with strong ordering of transverse momentum 
and as such, it accounts to only a small part of the available transverse phase space. This 
restriction on the transverse phase space becomes excessively prohibitive and must be 
lifted at very small x. The practical method of summing the leading-log - contributions 
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to the unintegrated gluon density JF(a;, k^) without restrictions on the transverse momenta 
of partons has been developed in 1975 by Fadin, Kuraev and Lipatov [35,36] and refined 
further by Bahtsky and Lipatov ( [37, 127], for the review see [128]). One has to pay a 
heavy price, though: the BFKL evolution receives a substantial contribution from soft, 
nonperturbative transverse momenta of final state partons, where the running strong 
coupling as is not small and the sensitivity to models of infrared-regularization can not be 
eliminated ( [129-133] and references therein). Although the fully satisfactory quantitative 
solution to this problem is as yet lacking, many of the properties of the QCD vacuum 
exchange must be regarded as well established: 

• Discard the asymptotic freedom, i.e., make the approximation as = const and allow 
the infinite propagation range for gluons. Such a model is free of a dimensional param- 
eter and possesses the scale-invariance property, which allows for an exact solution. 
The j-plane singularity of the model is a fixed cut (branching point) [35-37] at 

121og2 , , 

-oo < j<l + I^BFKL = 1 + —as 59) 

TT 

with vanishing a'^pj^]^ = 0, which is natural in view of the lack of any dimensional 
parameter in the model. 

• One can cope with the asymptotic freedom within the BFKL approach only at the 
expense of a certain regularization of the infrared growth of a^. One only needs to 
account for the finite propagation length, Rc, of perturbative gluons as suggested, for 
instance by the lattice QCD studies [50-52]. In their 1975 paper Fadin, Kuraev and 
Lipatov remarked that in this case the branching point is superseded by a sequence 
of moving Regge poles [35]. The positions of the poles were estimated in 1986 by 
Lipatov [127] 

A ^ ^BFKL 

A.^-^. (60) 

Herebelow, when discussing the pure Pomeron amplitudes, we shell refer to A„ as the 
intercept, which must not cause a confusion. Within the color dipole approach the 
poles differ by the number of nodes in the eigen-cross section as a function of the dipole 
size r [134]. The rightmost pole has a node-free eigen-cross section, the nodal structure 
of the eigen-cross sections and the n-dependence of the intercept of subleading vacuum 
poles found in [134] are very close to the quasiclassical approximation results of Lipatov 
[127]. The intercept of the rightmost pole Ajp, the slopes of the emerging Regge 
trajectories and positions of nodes in the eigen-cross sections depend on the infrared 
regularization ( [25,130,131] and references therein). 
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3.3 Poor man's approximations to the QCD Pomeron 



3.3.1 The Q^-independence of the Pomeron intercept 

For each and every pole the intercept does not depend on the probe. In apphcation to 
DIS that means an independence of intercepts on [129-131,135], only the residues can 
depend on Q^, so that the x-dependence of structure functions will be of the form 

F,{x,Q^) = 5^F(")(g2) (-Y" + F-f\Q') . (61) 

Examples of such a BFKL-Regge expansion for the proton and photon SF's with energy 
independent soft contribution F2°-^*(Q^), i.e., Asoft = 0, are found in [134,136-138]. If 
one reinterprets the soft contribution in terms of the soft, nonperturbative, unintegrated 
gluon density, then similar Regge-BFKL expansions hold for the integrated gluon density, 
G{x,Q'^), and the unintegrated gluon density 

J^{x, K ) — 



d log 



An example of the decomposition of J^{x, k^) into the soft and hard components is found 
in [34, 139] and is shown in Fig. 10. 

From the viewpoint of the energy dependence, the Regge cut also can be viewed as an 
infinite sequence of Regge poles. One can approximate the local x-dependence of the 
BFKL-Regge expansion (61) by 

F,{x,Q') = F{Q')i^-j , (62) 

which must not be interpreted that the Pomeron is a Regge pole with Q^-dependent 
intercept, for such an warning see, for instance, Bjorken [91]. An example of how the 
effective intercept A(Q^) changes with the range of x is found in [135, 143, 144], the vari- 
ations of the effective intercept from the unintegrated gluon density jF(x, k^) oc {^Y^'^ 

to the integrated gluon density G{x,Q'^) oc {^^^'^ and to the proton SF F2{x,Q'^) is 
found in [34,139], see Fig. ll,where we show separately the intercept for the hard com- 
ponents of J^{x, K?), G{x, Q"^), F2{x, Q^) and for the same quantities with the soft contri- 
butions included. These intercepts parameterize the local x-dependence for 10^^ < x < 
10"^. The striking finding is that while Thardii^"^) and XhardiQ"^) exhibit a very strong 
scale-dependence, i.e., the contributions form the subleading BFKL poles are large, the 
AhardiQ'^) is about independent one, AhardiQ"^) ~ 0.35-0.45. 
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F(x,/c'), D-GRV 
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k' (GeV^) k' (GeV^) 



Figure 10: The differential gluons structure function of the proton determined in 
[34, 139] from the k±-factorization analysis of the experimental data on F2p{x, Q"^). 
Notice the transition from the x-independent soft component at small k"^ < 1 GeV^ 
(shown by the dashed curve) to the hard component (the dotted curve), which con- 
verges to the derivative dGncLApix, n"^) / OXogn"^ of the integrated gluon density de- 
termined from the LO DGLAP fit to F2p{x, Q'^) . This particular example is for the 
GRV LO parameterization [I40], very similar results are found for the MRS [141] 
and GTEQ [I42] parameterizations. 
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Figure 11: The effective intercepts t{k,'^), X{Q'^), A((5^) of the local x- dependence, 
10~^ < a; < 10^^, of J^{x,K,'^),G{x,Q^),F2{x,Q'^), respectively for the k±- 
factorization analysis [34,139] with the large-K,"^ behavior o/ jF(x, k^) tuned to 
the GRV LO parameterization [I4O] as described in the text. In boxes (a)-(c) the 
dashed lines are for the hard components, the solid lines are found if the soft com- 
ponents are included. The box (d) shows how the intercepts change from jF(x, /t^) 
to G{x,Q^) to F2{x,Q^). The very close results are found for intercepts of pa- 
rameterizations tuned to converge at large to the MRS LO [I4I] o-nd CTEQ 
LO [142]. 
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3.3.2 The contributions from the soft region beyond pQCD 



Here one faces three major questions: (i) is the rise of soft hadronic cross sections driven 
by small dipoles in hadrons, (ii) what is the mechanism of interaction of non-perturbative 
large dipoles and (iii) is the soft contribution relevant to the large-Q^ DIS? 

The first question can be answered in the affirmative: the somewhat model-dependent 
estimates suggest strongly that the rise of the hadronic and real photoabsorption cross 
sections receive a large if not a predominant contribution from the interaction of small-size 
color dipoles in hadrons [34,136]. This suggests a weak energy dependence of the genuine 
soft vacuum exchange: Asoft ~ 0. The discussion of the potential importance of hard 
contributions to hadronic cross section was initiated in [113], for the recent work along 
these lines see [145]. 

From the color dipole viewpoint, the pure pQCD considerations stop at the dipole size 
^ -Rc ~ (0.2 0.3) fm and can not describe the bulk of the hadronic cross sections. It is 
plausible that at such large dipole sizes the color dipoles spanned between the constituent 
quarks do still remain the important degrees of freedom, but the corresponding soft dipole 
cross section remains a model-dependent phenomenological quantity, for which we only 
have constraints from soft hadronic diffractive scattering or from real or moderate-Q^ 
photoabsorption [24,34,136,138]. Such a soft dipole cross section can be modeled ei- 
ther by the non-perturbative two-gluon exchange [24, 34, 135, 136] or within the closely 
related model of the stochastic QCD vacuum suggested by the Heidelberg group [146]. 
Purely phenomenological attempts to guess the shape of this soft cross section and its 
continuation into the hard region [147] should not be disregarded as well. 

From the practical point of view, the available models for the dipole cross section suggest 
a smooth r-dependence of the dipole cross across r ~ i?c up to r ~ 1 fm. Because the 
pQCD BFKL component of the dipole cross section rises with energy much faster than 
the energy-independent soft dipole cross section, at higher energies the dominance of the 
pQCD component of the dipole cross section will extend beyond r ~ Rc, for which reason 
the lower boundary for the pQCD dominance will be lower than given by eq. (13). One 
can come to the same conclusions from the smooth K^-dependence of the unintegrated 
gluon density from soft to hard region and the dominance of the hard component at large 
i which is clearly seen in Fig. 10. 

Regarding the question (iii), even at very large the virtual photons contain the hadronic 
size qq components and the SF's receive a non-vanishing, even substantial at x ~ 10~^, 
contribution from the interaction of soft dipole. Within the more familiar DGLAP ap- 
proach such a contribution is hidden in the input parton densities; the sensitivity of the 
DGLAP evolution to the input partons is an old news, although eventually the rising 
perturbative QCD component would take over at very large [134-136]. Recently there 
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were many suggestions to start with the Regge parameterization of photoabsorption at 
small to moderate < Ql and take F2^^^^'^\x, Ql) as a boundary condition at = Ql 
for the DGLAP evolution at large Ql ( [148-150] and references therein). 

3.3.3 The two-Pomeron approximation 

The transition from the unintegrated gluon density, /t^), to the conventional, inte- 
grated one, G{x,Q'^), involves an integration, G{x,Q'^) = J {dti^ / k'^)T{x, h?). Simi- 
larly, to the DGLAP approximation the small-x SF involves an integration, F2(x, Q^) oc 
J*^ {dn? / k'^)G{x, H?). Each integration shifts the nodes to larger value of and, fur- 
thermore, enhances the relative contribution from the node-free rightmost eigen-function. 
The model-dependent estimates within the color dipole model show that the QCD vac- 
uum exchange contribution to DIS is numerically dominated by the rightmost Pomeron 
pole plus the energy-independent soft exchange contributions ^ because the subleading 
Pomeron pole contributions have a node in the practically important region of ~ 10-^40 
GeV^ [130,135,136,138,144]. This is the reason behind the remarkable flat dependence 
of ^hard{Q'^) shown in Fig. 11. Consequently, within the kinematical range of HERA, the 
hard contribution to the proton SF can be well approximated by a simple Regge-pole 
formula with the intercept ^hard ~ 0.35-0.45 [34, 139]. This finding is a dynamical jus- 
tification of the two-pole approximation [145,154,155]. The specific models [34,136,138] 
give the concrete Q^-dependence of the residues; on general grounds there are no reasons 
for decoupling of the effective hard Pomeron from soft amplitudes, including the real pho- 
toproduction. We emphasize that the two-Pomeron parameterization only holds in the 
limited range of x and should not be extrapolated far beyond the kinematical range of 
HERA. 

3.4 The basics of the theory of diffractive vector meson produc- 
tion 

Here we comment briefly on properties of diffractive vector meson production starting 
with the nonrelativistic quark model in conjunction with the vector dominance model. It 
offers a useful insight into such fundamental issues as the flavor dependence, the relation 
between the vector meson production and — > e^e~ decay and the way the short 
distance wave function of vector mesons is probed in vector meson production. Then we 
qualify those properties in the color dipole approach. 

^ To this end it is instructive to recall the early doubts in the necessity of the hard Pomeron contribution 
for description of the observed cross sections [151-153] 
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3.4.1 The flavor dependence, the relation to the decay e+e 
and VDM 

On the one hand, the V'^ e^e~ decay amphtude can be parameterized in terms of the 
matrix element of the electromagnetic current 



(0| J,JV) = -y/A^,gvCvV^ , (63) 
where is the vector meson polarization vector, so that the decay width equals 

Y(v'^e^e-) = ^^^^f4^- (64) 

Here the charge-isospin factors cy are Cp = ^(e„-ed) = c^j = ^(en+e^) = = 
e<j = — |, cj/* = Cc = |, cx = Cft = — |. One also often uses the parameter 

J_ _ gycy 
fv m\ ■ 

On the other hand, in the nonrelativistic quark model the vector meson is the weakly 
bound spin-triplet, S-wave qq state, and the decay —>■ e^e~ proceeds via annihilation 
qq — ^ e+e", 

T{V' e+e-) = \Rym\v,,a{qq^e^e-)) = ^^|i?^(0)^ (65) 

my 

where Vqq is the relative velocity of the quark and antiquark in the vector meson and 
Rv{^) is the radial wave function at the origin [156]. This gives a useful relationship 



gv = Rv{^)f-^. (66) 
which amounts to the nonrelativistic calculation of the Feynman diagram of Fig. 12. 




Figure 12: The decay of the vector meson into the lepton pair via annihilation 
qq — >• e~^e~ . 
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Consequently, in the simplified VDM approximation, for transverse photons. 





(67) 



31(^0 e+e-) 




^r{Vp Vp) 



(68) 



myaem. 



Q"^ + m' 



'V 



cvRv{0)y/12aemmv 



(69) 



Precisely the same result is found if one computes the vector meson production amplitude 
through the diagrams of Fig. 13 and applies the additive quark model. 



Figure 13: The Additive Quark Model approximation for the vector meson pro- 
duction amplitude. 



In the case of the pp and ujp final states a very good parameter free description of the 
E401-FNAL measurements of the differential cross section of photoproduction is found 
if one takes isoscalar elastic vrA^ scattering amplitudes for TiVp — > Vp) ( [43, 157] and 
references therein). The sp, sp amplitudes needed for the ipp state can be extracted from 
the ttA^, KN., KN elastic scattering amplitudes 



which gives a perfect description of the t-dependence of the E401-FNAL data on photo- 
production of 0p [157]. Specifically, (71) correctly reproduces the experimentally observed 
change of the diffraction slope from B{^p — ^ ujp) = 12.6 ± 2.3 GeV~^ to B{'yp — »• (pp) = 
6.8 ± 0.8 GeV~^. In terms of the discussion in Section 3.1.3, see Eq. (35), this inequality 
of diffraction slopes suggest that the spatial size of the made of the heavier strange 
quarks is substantially smaller than the spatial size of the u made of the light u, d quarks. 
However, the observed differential cross section is only a half of what is predicted by 
(68) and (71). Within the color dipole approach the culprit is the oversimplified VDM 
approximation (71): the interaction of the quarkonium is controlled by not the number 
and flavor of quarks in the state but rather its size [158]. 




TiVp Vp) = T{qp qp) + T{qp — > qp) . 



(70) 



T{(f)p (pp) = Ti^K'^p K'^p) + T{K p ^ K p) — T^n p ^ n p) 



(71) 



40 



3.4.2 Vector meson production in the color dipole approach 



In the color dipole approach, thanks to Lorentz dilation of time at high energies, the 
partonic fluctuation (to the lowest order, qq pair) of the incident photon is frozen in 
transverse (impact parameter) space during the interaction with the target. This allows 
one to cast the photoproduction amplitude in a quantum-mechanical form [17,19,21] 

T = {^v\^dip\^^) = J dzcfr <iJ*y{r) adipix,r) ^^(r) , (72) 

where z and (1 — z) are fractions of the photon's lightcone momentum carried by the 
quark and antiquark, respectively. The basic quantity here, the cross section of the color 
dipole interaction with the target adip{r), can be calculated for the forward scattering 
case through the unintegrated gluon distribution, 

adipix, r) = ^ J ^J'(x, k) as[max{K'^, A/ r^)] [1 - exp{iKr)] , (73) 



where 



A ~ 9 ^ 10 (74) 



follows from the properties of Bessel functions [33] . Eqs. (72) and (73) sum to the 
leading log^ the towers of two gluon exchange diagrams of Fig. 9, as manifested by the 
unintegrated glue J^{x, k) in the integrand of (72). The x-dependence of the dipole cross 
section is governed by the color dipole BFKL equation ( [129,159,160], see also [161,162]), 
for the discussion of the choice 

x = x,-0.4.^^±^ (75) 

see below Section 4.6. 

In due turn, the unintegrated glue of the proton can be extracted from the experimental 
data on the proton structure function [34, 139], so that there is a microscopic QCD link 
between inclusive DIS and vector meson production, if the vector meson is treated in 
the qq Fock-state approximation. For small dipoles there is a useful relationship to the 
integrated gluon structure function of the proton [33, 163] 

adipix, r) = yr^as (^^^ G (^x, . (76) 

Now comes the crucial point: the lightcone wave function of the virtual photon shrinks 
with Q^, namely, "^-yir) oc exp(— er), where [17-19] 

= z{l - z)Q^ + , (77) 
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where is mass of the quark of the flavor /. Then, for a sufficiently large Q^, the 
dominant contribution to the virtual photoproduction amplitude will come from r ~ 
= 3/e, so that 



where in the integrand of (78) the ^-dependent factors coming from the photon wave 
function have been suppressed. 

Note that the " quark mass" term mj here must not be omitted even for the light flavors. 
This " quark mass" serves as an effective parameter that bounds from above the transverse 
size of the qq state in a real photon. One can discuss the large-size properties of the photon 
only under certain assumptions on the color-dipole cross section for large dipoles or the 
unintegrated gluon density for nonperturbative soft gluon momenta: the early choice has 
been rriu^d ~ 0.15 GeV [134, 135], the more recent /c_|_-factorization analysis [34] of the 
low-Q^ data suggests m^ ^ ^0.22 GeV. 

The result (78) has all the properties of the amplitude (69) subject to important QCD 
modifications: 

• The color-dipole cross section is fiavor independent, and the charge-isospin factors are 
precisely the same as in the VDM. 

• For rs ^ Rv the vector meson production is obviously short distance dominated and 
tractable within pQCD ( [17-19,22,23], for refinements on the applicability of pQCD 
see Collins [164]). The amplitude is proportional to the vector meson wave function at 
vanishing transverse qq separation, \E'y(z,0), which is closely related to the so-called 
vector meson distribution amplitude [165,166]. 

• To the nonrelativistic approximation, 2; ~ 1/2 and my ^ 2mq, one has ~ |(Q^ + 
my), and the factor r| oc 1/{Q'^ + ^v) reproduces the dependence dictated by the 
vector meson propagator. 

• However, atotiVp — > Vp) which enters (69), is substituted for by 



where we used (74) and (9) by which A/vg ^ Q . For large dipoles, rs ~ -Ry, which 
dominate in real photoproduction, adip{Rv) ~ crtot{Vp Vp), but for small dipoles, 
rs <C Rv, which dominate electroproduction, o'dip{xg,rs) ^ o-totiVp Vp) and the 
simplified VDM is bound to fail. 

• For small scanning radii, rs -C Rv, such that 'ifv{z,rs) ~ const, the dependence on 
and the mass of the vector meson my only enters through the scanning radius rs- 




(78) 



Cdipixg, r s) 



Q 



r 2 2 

:^as{Q )G{xg,Q ) , 



(79) 
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Hence the fundamental prediction [21] that cross sections for different vector mesons 
taken at the same value of rs, i.e., the same value of (Q^ + my), must exhibit similar 
dependence on energy and {Q^ + niy). 

Vector meson production probes the integrated gluon SF of the target proton at hard 

2 2 

scale Q given by (8) ( [20-22], for a more accurate definition of Q for light vector 
mesons see [38]). 

Notice an inapplicability of the simplified VDM to heavy quarkonia, for which by 
virtue of small as the Bohr radius 

4 

Rv = as ^ > rs ■ 

myas 

Finally, as far as the t-dependence is concerned, rs can be regarded as the transverse 
size of the 'j* —>■ V transition vertex, so that for the fixed value of x the diffraction 
slope is predicted [25, 26] to decrease with (Q^ + my 



const 



5(g2) r^BN + Crl ^Bn+ ^;7^% . (80) 

+ m^ 

Because the color dipole cross section and the unintegrated gluon SF are related by 
the Fourier transform, all the above results can be rederived in the momentum space 
representation, often referred to as the /c^-factorization or impact factor representation. 
The relevant formalism goes back to the 1978 seminal paper by Balitsky and Lipatov [37], 
although the term " /c^-factorization" has been coined much later on by several groups 
[167-169]. The detailed application of the /c^-factorization to the vector meson production 
is found in [54,170-172] and will be reviewed in the following section 4, the first momentum 

2 

space derivation of the leading log Q approximation is due to Ryskin [22] and Brodsky et 

2 

al. [23] , some corrections to the leading log Q approximation were discussed by Levin et 
al. [173]. Referring to Section 4 for a detailed discussion of the helicity amplitudes within 
/c_L-factorization, here we only cite the gross features of the longitudinal and transverse 
cross sections: 



1 



asiQ')Gixg,Q') , (81) 



asiQ')G{xg,Q') . (82) 



Here the factor ~ Q'^/my in the is a generic consequence of the electromagnetic gauge 
invariance, as has been understood in early 70's [174,175]. 

The Heidelberg group [176] starts with the soft color dipole cross section evaluated within 
the stochastic QCD vacuum model [146]. It shares with other color dipole models the 
predictions for the dependence, but the energy dependence does not follow from the 
first principles of the model and needs to be introduced by hand [177]. 
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3.4.3 Production of excited vector mesons 

The p^,uj^,(f)^ and J/ip are the ground state vector mesons. The \l/'(3686) is the well 
established radial excitation 25'-state, the \I^"(3770) is a solid candidate for the orbital 
excitation D-wave state [47,48], the radial vs. orbital excitation assignment in the p,uj,(j) 
family is not definitive yet [102]. 

The salient feature of the 25* radial excitations is a node of the radial wave function, 
'^2s{z,r), at r = Vnode ~ R{1S) = Rv, which suppresses the V'{2S) production amplitude 
in comparison to the corresponding ^^(15*) production amplitude [17,18,49,178]. The 
strength of the node effect depends on the proximity of the scanning radius rs to the 
node position Vnode- At rs <^ rnode (in the under-compensation regime), which can take 
place at high or for very heavy mesons, the contribution from r > Vnode is small and 
suppression is weak. The under-compensation regime is relevant to the \1/'(2S') production 
where the color dipole model predicts the rise of the ratio a{'^'{2S))/a{J/ip{iS)) with 
rising Q^. For light vector mesons at small the over- compensation scenario of > Vnode 
and strong cancellation is not excluded [26,49,178]. In this scenario the ^'(25*) and V^(1>S') 
production amplitudes will be of the opposite sign, which can be tested experimentally 
via the Soding-Pumplin effect [179,180], and the differential cross sections da(y' (2S))/dt 
may exhibit a sharp forward dip [26,49,54]. In such a regime even a small shift of 
would strongly alter the cancellation pattern, giving rise to an anomalous dependence 
of the ratio a{V'{2S))/a{V{lS)), of the t-dependence of da{V'{2S))/dt and of the ratio 
ai/cTT for the ^^'(25) - the latter effect is due to a slightly different impact of the node 
effect on different helicity amplitudes. A subsequent discussion of sensitivity of the node 
effect to the wave function of vector mesons is found in [181-184], the change of numerical 
results for the \l/'(25') from one model to another must be regarded as marginal. 

The case of the orbital excitation V"{D) is quite different [171]: here the radial wave 
function vanishes at the origin, and the dependence of the V"{D) production will be 
smooth. There are some subtle changes in the helicity amplitudes: in both the V{1S) 
and V"{D) the qq pair is in the spin-triplet state, but the total spin of the pair is along 
in V{1S), and opposite to in V"{D), the spin of the meson. 

The node effects echoes in the hard scale for the ^'(25*) production. In the under- 
compensation regime of relevance to the \1/'(2S') the contribution to the production am- 
plitude from large color dipoles, r > Vnode, is canceled by the contribution from small 
dipoles, r < Vnode- As a result, the \1/'(2S') production amphtude is dominated by color 
dipoles of smaller size than it is the case for the J/\l/(15') and color dipoles models pre- 
dict the hierarchy of hard scales g^(*'(25)) > g^(J/*(15)). Consequently, the ^'(2S) 
production amplitude must grow with energy faster than the J/\I/(1S') production ampli- 
tude [49]. Furthermore, the (negative valued) contribution to the production amplitude 
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from large dipoles, r > r^ode, has a steeper t-dependence than the (positive valued) con- 
tribution from small dipoles, r > Tnode- As a result, the diffraction slope in the ^'{2S) 
production is predicted to be smaller than in the J/'^{1S) production [26,49]. 

3.4.4 Unitarity and saturation in the color dipole language 

The unitarization of rising scattering amplitudes in QCD remains one of the hot and 
as yet unsolved issues. As emphasized in Section 3.1.4, the unlimited growth of the 
model partial waves must be tamed and the unitarity bound r(b)<l must be met in a 
consistent treatment of high energy scattering. The theory is still in the formative stage, 
though. Some of the early works on unitarization have been mentioned in Section 3.1.4, the 
problem of unitarity is most acute for interactions with nuclei, in which case the impulse 
approximation partial waves ro(b) oc A^/^. For the nuclear targets the presence of a new 
large parameter — the optical thickness of a nucleus — leads to certain simplifications 
like the applicability of the eikonal approximation for the color dipole-nucleus scattering 
[24,185-189]. The recent development in imposing the unitarity on nuclear amplitudes, 
often referred to as the color glass condensate, is summarized in [190-193], for a review 
of the early works see [194,195]. A review of the enormous literature on the subject goes 
beyond the scope of this review, we rather present a brief introduction into major ideas. 



Figure 14: (a) The multigluon t-channel exchange diagram contribution to the 
color dipole scattering amplitude and (b ) its approximation by multiple exchange by 
the two-gluon Pomerons. 

Let ro(r, b) be the profile function for the color dipole-nucleon scattering evaluated in 
the single-Pomeron exchange approximation of Fig. 9. The Gaussian approximation (48) 
is not imperative but convenient for the sake of illustration. In the full fledged QCD 
one needs to sum all multigluon t-channel exchanges between the color dipole and nu- 
cleon, including interactions between all exchanged not shown in Fig. 14a. A poor man's 
approximation to this as yet unsolved problem ( [196, 197] and references therein) is the 
multiple exchange by bare Pomerons, in general case the interactions between gluons from 



dipole 




b) 
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different Pomerons must be included. When the multipomeron exchanges are evaluated 
in the eikonal approximation, one obtains the "unitarized" profile function [163] 



r(r,b) = l-exp[-ro(r,b)], (83) 
whereas the so-called /^-matrix unitarization gives [163] 

The latter has been suggested also from the consideration of the so-called fan diagrams 
( [198] , similar results are found from different approximate non-linear evolution equations 
[199,200]), for the so-called tZ-matrix approach see [201]. 

The principal point is that the unitarized partial waves do always respect the unitarity 
bound ro(r, b)<l. The partial waves saturate at the black-disc limit, r(r, b) ^ 1, for all 
impact parameters such that To{r, b) ^ 1, i.e., 

b2<2fi(r)logro(r,b = 0). (85) 

The two unitarized forms (83) and (84) only differ by the rate of approach to the black disc 
limit, the i^-matrix unitarized dipole cross section takes a particular simple form [163] 

a{x, r) = ATTBir) log (^1 + , (86) 

which shows clearly how the power like small-x growth of the bare Pomeron cross section 
ao{x,r) oc is superseded by the oc log ^ behavior, or oc log^ ^ if one allows the 

Regge growth of the diffraction slope B{x, r) [202]. 

Finally, the unitarization alters dramatically the r-dependence of the dipole cross section 
from (76). At asymptotically small x the unitarization is at work already for small dipoles, 
where the r-dependence of diffraction slope B{x,r) can be neglected, see (80). so that 
the dipole cross section would saturate, a{x, r) ^ AirB, for dipoles 

>r^ = ]^ f87) 

~ 7ras{r)G{x,q^ = A/r^) ' ^ ^ 

The smaller is x, the larger is the gluon SF in the numerator in the r.h.s. of (87) and the 
smaller is the saturation scale r^^^. Recently, specific parameterizations for the saturating 
dipole cross section without an explicit reference to the unitarity properties of partial 
waves have been proposed [147,203]. (In principle, the saturation rate and the saturated 
cross section must be adjusted to describe the diffractive hadronic scattering and real 
photoproduction [3,24], which has not been done in the model [147,203].) However, with 
the realistic dipole cross sections the unitarization effects for DIS [163] and for vector 
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meson production [21] were found to be marginal. The extraction of the S'-matrix for 
the color dipole scattering from the vector meson production data by Munier, Mueller 
and Stasto also shows that the dipole-nucleon scattering is not yet close to the strong 
absorption regime [204]. Similar conclusion follows from the impact parameter extension 
[205] of the saturation model [147,203]. Those findings are not surprising, though: as 
shown in [4] in the limit of strong saturation the diffractive rapidity gap DIS must make 
precisely 50 per cent of the total DIS cross section, whereas experimentally the fraction 
of diffractive DIS is about 10 per cent [1,2]. 

To summarize, as soon as impact parameter dipole model has been adjusted to fit the 
experimental data on DIS structure functions and the total cross section and the t- 
dependence of diffractive vector meson production, it is expected to have partial waves 
consistent with the unitarity constraints in the energy and range in which the exper- 
imental data are available. The same must be true of the unintegrated gluon SF of the 
proton extracted in [34] from the DIS data. Applying unitarity corrections to the vector 
meson production amplitudes evaluated with such an unintegrated gluon SF would be the 
double counting. 

3.4.5 Color dipole model and Generalized VDM 

The simplified VDM must be regarded as the leading term of the mass-dispersion relation 
calculation of the dependence of the virtual photoproduction amplitude. The impor- 
tance of contributions from the more distant singularities — the higher vector states and 
the continuum, which we denote generically as Vi, — rises with [174,175,206]. Within 
the resulting Generalized VDM (GVDM) for DIS the calculation of 7*p — > 7*p must allow 
for transitions of photons to all higher vector states, 'y* ^ Vi, followed by the diagonal 
and off-diagonal scattering Vip Vjp and the transition Vj — 7*, see Fig. 15. 



Figure 15: The Generalized Vector Dominance Model diagrams for Compton 
scattering (DIS) and diffractive vector meson production. 

Similarly, the transitions 7* — ^ followed by the off-diagonal scattering ViP — ^ Vp would 
contribute to the vector meson production 7*p — > Vp. If viewed as the mass dispersion 
relation, the GVDM can not fail, but the practical application requires the knowledge 
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of all the diagonal and off-diagonal amplitudes Vip — > Vjp and of the dependence of 
transitions 7* — > VJ. The color dipole model provides the QCD input for the GVDM 
analysis [18,207], the equivalence of the two approaches emphasized in [24] has been 
elaborated by Schildknecht et al. [208]. 

3.4.6 The s-channel helicity non-conservation (SCHNC) 

The scattering of the qq dipole on the target via exchange of the two-gluon tower exactly 
conserves the s-channel helicity of the quark and antiquark (for the QED case see [209, 
210]). This does not imply the conservation of the helicity of photons in the off-forward 
Compton scattering. As a spin-1 particle, photon is similar to the deuteron. In the non- 
relativistic case the pure S-wave deuteron with spin up consists of the spin-up proton and 
spin-up neutron, the longitudinal deuteron consists of the spin-up proton and spin-down 
neutron and vice versa. 

The perturbative QED transition of the photon to the qq pair is described by the familiar 
vertex Cfq'-y^^qA^. The longitudinal (scalar) virtual photon with helicity A.^ = consists 
of the qq Fock state with A + A = = 0, in close similarity to the S'-wave deuteron. 
The crucial point is that the transverse photon with helicity A^ = ±1 besides the qq state 
with A + A = A^ = ±1 contains the state with A -|- A = 7^ A^, in which the helicity of 
the photon is carried by the orbital angular momentum in the qq system (see [211] for an 
early discussion of this mechanism in application to the spin- flip in the nucleon scattering) . 
Furthermore, it is precisely the state chiral-even state with antiparallel helicities, A+A = 0, 
which gives the dominant contribution to the absorption of transverse photons and the 
proton SF F2p{x, Q"^) in the Bjorken limit. The perturbative transition of transverse 
photons to the chiral-odd state with parallel helicities, A + A = A^ = ±1, vanishes in the 
massless quark limit. 

The helicity structure of vector mesons is about the same. From the point of view of 
the vector meson production, it is important that the transverse and longitudinal 7* and 
V share the intermediate qq state with A + A = 0, which allows the s-channel helicity 
non-conserving (SCHNC) transitions between the transverse (longitudinal) 7* and lon- 
gitudinal(transverse) vector meson V [170,212]. This mechanism of SCHNC does not 
require an applicability of pQCD. 

Hereafter we only discuss the experimental data from HERA taken with unpolarized 
protons, hence proton can be treated as a spinless particle, see, however, a brief discussion 
in Section 5.5. Depending on the spin-parity of the t-channel exchange, the helicity 
amplitudes satisfy [213] 

T_,,,_,^ = ±(-l)"--^-r,^,^ , (88) 
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where the +(— ) sign apphes to natural (unnatural) parity exchange. As discussed in 
Section 3.1.2, Pomeron and all the highest lying subleading reggeons have the natural 
spin-parity. Under the dominance of the natural spin-parity exchange, the number of 
independent helicity amplitudes is reduced to five: 

L^L; T^T{X^ = Xv) 
T L; L-^T 

T-^r (A^ = -Ay). (89) 

The first line contains helicity-conserving amplitudes. They are predicted and found to be 
the dominant ones. They do not vanish for the forward production, A = 0. The second 
line in (89) contains two single helicity-flip amplitudes. They must be proportional to 
|A| in the combination (e ■ A) or (V* ■ A), since there is no other transverse vector at 
our disposal. The last line contains the double helicity-flip amplitude, which must be 
proportional to (e ■ A)( V* - A). 

One can thus predict that the s-channel helicity conserving amplitudes will dominate 
in the almost forward production of mesons. As t increases, the relative importance of 
helicity-fiip amplitudes will grow, and, at high enough t, they might become competitive 
to the helicity conserving amplitudes. 



3.4.7 Diffractive vector meson production from extended Bloom- 
Gilman duality. 

The Bloom-Gilman inclusive-exclusive duality relates the a: — 1 behavior of DIS to elastic 
ep scattering [214]. Roughly speaking, if one stretches the x-dependence of the DIS cross 
section determined for the continuum masses W to the elastic limit W — nip, then the 
DIS cross section integrated over the interval 

0<l-x< ° ^ (90) 

will, with the judicious choice of the duality interval [m — p, Wq] , be equal to the elastic 
ep cross section (for the recent active discussion of duality in DIS in connection with the 
JLab data see [215]). Genovese et al. argued [216] that similar parton-duality relationship 
must hold between the diffraction excitation of the small mass continuum 

ep — > eXp 

and exclusive vector meson production. In terms of the diffractive Bjorken variable 
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the diffractive structure functions for the transverse and longitudinal photons have the 
large-/? behavior [216-218] 



FT(x,/3,g2)(x(l-/5)2G2(x,g2), (91) 

1 



FT{x,(3,Q^)^—G\x,ql). (92) 



The relevant hard scales equal [216,217,219] 

TTl^ 1 

- J^iQ" + M'x) , ql--{Q' + Ml). (93) 

The integration over the duality interval [Mmin ~ 2mq,Mr^i], i.e., 1 — [3t,l < M^j^^/Q'^, 

yields the correct large-Q^ dependence of aL,T- Furthermore, both hard scales tend to 
— 2 

the scale Q of Eq. (8) so that Eqs. (91), (92) yield precisely the same dependence on the 
gluon structure function as in (80), (81). Motivated by this observation [216], Martin et 
al. suggested to evaluate the vector meson production cross sections (Tt,l from the duality 
integral [220,221]. This way one encounters a very strong sensitivity of such evaluations 
of aT,L to the duahty interval, 

ar cx (M| - M^J , (94) 
aL (X {Ml - MlJ , (95) 

which is especially strong in the case of ax- 

Similar in spirit to the duality is the unorthodox color evaporation model (CEM). In 
its original formulation [222] it simply states that the color of the qq pair produced in 
1*9 ~^ subprocess happens to be bleached by soft final-state interactions leading to the 
rapidity gap events with the probability 1/9. Within CEM the charmonium production 
is described by the formation of colored open charm cc states which masses Mcc<2m£), 
where rriD is the mass of the D-meson [223]: 

(ionium = 7: , (96) 

9 J dMcc 

2m.c 

where 1/9 is the color bleaching probability. Assuming that about 50% of the onium 
goes into the J/\E', Amundson et al. are able to describe the photo- and hadroproduction 
of the J/\E' [224,225]. Gay Ducati et al. find similar agreement with the total cross 
section of elastic charmonium photoproduction [226,227]. Here we only notice that in 
the near-threshold process 7*5' cc produces the spin-singlet S-wave cc pair. Arguably, 
the color bleaching can not flip the spin of nonrelativistic heavy quarks and that must 
lead to strong suppression factor in the estimate (96). In contrast to that, in hadronic 
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collisions the near-threshold open charm can be produced in the spin-triplet state via 
QQi 99 ~^ 9 ~* cc and the spin dynamics of heavy nonrelativistic quarks does not prohibit 
the formation of J/"^ by color bleaching. 

3.4.8 Models which respect the Froissart bound 

Only a limited range of energy, x, and is spanned by the available experimental data. 
We already mentioned of an equally good description of the available soft cross section data 
by the soft Pomeron pole exchange and logarithmic parameterizations, see Section 3.1.4. 
We also recall an observation by Buchmiiller and Haidt [228] that gross features of the 
small-x proton structure function as measured at HERA are reasonably well reproduced 
by a very simple parameterization 

F,{x,Q') = a + m\og^\og%. (97) 

From the Regge theory viewpoint this corresponds to the dipole singularity at j = 1. The 
dipole singularity model for virtual photoproduction of vector mesons has been proposed 
by Fiore et al. [229,230], in this specific example the nonlinear Pomeron trajectory with 
the branching point singularity at the two-pion threshold in the t-channel is used. For 
each and every vector meson a good description of the vector meson production cross 
sections is found at the expense of five free parameters. A very closely related model was 
proposed by Martynov et al. [231,232]. Troshin and Tyurin suggested a parameterization 
of vector meson production amplitudes in which the high energy growth is tamed by the 
tZ-matrix unitarity constraints [201]. Haackman et al. [233] start with the soft Pomeron 
with Ajp > and impose the unitarization by reggeon field theory methods as mentioned 
in Section 3.1.4. The drawback of such models is that the dependence of the vector 
meson production is parameterized rather than predicted from the microscopic QCD. 
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4 The /cx-factorization: unified microscopic QCD de- 
scription of DIS and vector meson production 



4.1 The leading log^ and qq Fock state approximations 

The color dipole and /ct-factorization approaches to small-x DIS are conjugate to each 
other, the technical correspondence is given by Eq. (73). The advantage of the former is 
in its simple quantum-mechanical representation, still some technical issues such as the 
definition of the lightcone wave functions, the separation of the S'-wave and D-wave states 
of vector mesons, and the role of the so-called skewed, or off-diagonal, gluon distribution 
functions are more transparent in the momentum-space representation. 

The starting point is the BFKL diagram for small-x DIS, Fig. 9 and the reference reaction 
is the non-forward Compton scattering 7*p 7*(A)p(— A). The vector meson produc- 
tion is obtained from the Compton diagram replacing the outgoing pointlike photon 7* by 
the non-pointlike vector meson V. To the leading log ^ the effect of perturbative higher, 
qqg, qqgg etc., Fock states in the pointlike photon amounts to the BFKL evolution of 
the color dipole cross section or of the unintegrated gluon SF while retaining the qq Fock 
state approximation [35-37,129,161,162]. Namely, in the DIS counterpart of (72) one 
calculates the photoabsorption cross section as an expectation value of the dipole cross 
section over the lowest qq state of the photon: 

1 

crtot{l*p) = j dz j d^r^!:^,{z,r)adip{x,r)^!^.{z,r). (98) 




4.2 The hehcity and chiral structure of the photon 

In the momentum representation the chiral structure of the q'J^qA^ vertex is as follows. 
The photon polarization vectors are described in Section 2.3.2, here we only notice that 
in the Sudakov representation 



(99) 



where the two Sudakov lightcone vectors are defines as 



P = p' + ^q'; q = q'-xP; q'^ = p'^ = x = « 1. (100) 

Hereafter it will be convenient to use twice the quark and antiquark helicity. A, A = 
±1, which should not cause a confusion. For the transverse photons, A^ = ±1, in the 
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momentum representation the perturbative QED vertex gives the structure 

Qxlf^Qje^iX^) = ^ \-V2mf6x^^xSx,x + 26^ j[z6x„x - (1 - z)6^ j]{k ■ e(A^))| 
Vz{l - z) ^ ^ J 

(101) 

and for the longitudinal (scalar) photons 

qxl^.(txe,{X-y = 0) = -2Qy/z{l - z)6^^^j . (102) 

Here z and (1 — z) are the fractions of the photon's lightcone momentum carried by 
the quark and antiquark, respectively, and k and —k are the corresponding transverse 
momenta. The perturbative chiral-odd component of the transverse photon with parallel 
helicities vanishes for massless quarks. The scaling contribution to the DIS structure 
function F2{x,Q'^) comes from the chiral-even component with antiparallel helicities. 



4.3 The lightcone hehcity and chiral structure of vector mesons 
and rotation invar iance 

In the vector meson, the quark and antiquark are in the spin-triplet state and either S- or 
D-wave. The lightcone wave function \E'y(^;P) is a probability amplitude for expansion 
of the vector meson in qq states with invariant mass 

= — f = 4(m2 + p2) . (103) 

Z[l — z) 

One calculates first the amplitude of production of the qq pair, 

7> ^ {qq)p' , (104) 

and then projects it onto the vector state by weighting with \E'\/(z, Ic) and the relevant 
helicity factors. The use of the 3-dimensional momentum of the quark in the qq rest 
frame, p = (k, {z — |)M), 

d^p dz(Pk 

is helpful to see a link to the conventional quantum-mechanical description. 

The helicity/chiral structure of the vector meson for the widely used Vqq extension of the 
QED vertex of the form 

y^lqn^l(lf^v{z,k) , (106) 

is the same as for the photon subject to the substitution Q — ^ M for the longitudinal 
vector meson. The vertex (106) gives a certain admixture of the S and D waves. The 
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S D-mixing is familiar from the case of the deuteron, where it originates from the pion- 
exchange tensor interaction, the presence of the potential-dependent S'D-mixing in vector 
mesons is a generic feature of potential models (for the review see [47]). 

The rotation-invariant lightcone description of the pure S and Z^-wave states and the 
corresponding vertices S'^ and are found in [54,171], for the related discussion see 
also [234,235]. To generate the pure 5- wave state one needs to add the generalized Pauli 
vertex. Upon applying the Gordon identities, the pure S'-wave vertex can be cast in the 
form 

Ts{z, k)qfS,qfV, = Tsiz, k)qj - ^ ,d Vf " Vj),] qjV, (107) 



(M + 2m) 



with the helicity/chiral structure 

1 



?A5'^gA^/.(0) 



+ 



2{k-V{Xy)) 



M + 2m 



{-v^m;5,,,A,A + ^h,^-x[^K,x - (1 - ^)K-,x]ik . V{Xv)) 
mfil~2z)6,^_j + V2ik-V{-\))5,,^ 



f 



-2M./^(l^6,^_j 

M(l - 2z) 
{M + 2mf)y/z{l-z) 



mjil-2z)6,^_j + V2{k.Vi-X))6,-, 



:i08) 



Note that momenta pf and pj correspond to on- mass-shell fermions, see details in [54], 
which justifies the usage of the Gordon identity. The corresponding vertex functions Ts 
will only depend on the "radial" variable and can be related to the momentum-space 
radial wave functions ips{z, k): 



(109) 



An important part of the rotation-invariant description is that the transversity condition 
must be imposed at the level of the qq pair, which leads to the concept of the running 
longitudinal polaiization vector Vl{M), which has the Sudakov expansion 



(110) 



such that it is orthogonal to the 4- momentum of the on-mass shell qq pair, (f^g ■ Vl(M)) 
0, where 

Vqq - q + -^^f -qq 

This running polarization vector has been used in (108) 



P' v, = M\ (111) 
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The lightcone extension of the considerations in Section 3.4.1 gives the — > e~^e decay 
constant for the S'-wave state [54] 

1 

= / ''-^ I (2.)3t(l-.) '^^<P) ■ jAf (M + 



1 

= 9v J dz(t)v{z) , (112) 



where 4>v{z) is the so-called distribution amplitude for the pure S'-wave vector meson. 
The general phenomenology of distribution amplitudes can be found in [166, 236, 237] 
and [238]. If the vector meson is saturated by the qq state, then (112) is supplemented 
by the normalization condition 

1 - (^/^'P «I^^(P^)P - ^ / (113) 

The calculations with the fixed longitudinal polarization vector defined for fixed M = my 
break the rotation invariance, which is often the case with parameterizations used in the 
literature [21,26,204,239-241]. In technical terms, the fixed polarization vector leads to 
a mixing of the longitudinal spin-1 state and spin-0 states. One of the drawbacks of the 
fixed polarization vector is that the V — > e+e~ decay width would depend on the polar- 
ization state of the vector meson, while the rotation invariant approach with the running 
longitudinal polarization vector guarantees that the decay constants for the transverse and 
longitudinal vector mesons are identical. Quite often, in the fixed-polarization-vector ap- 
proaches, the manifestly different radial wave functions are introduced for the transverse 
and longitudinal vector mesons [181,205]. 

The principal effect of the Pauli vertex in the helicity/chiral expansion (108) is the chiral- 
odd parallel-helicity component of the transverse vector meson which does not vanish for 
massless quark. Going back from vector mesons to real photons, Ivanov et al. [242] argued 
that the related nonperturbative chiral-odd component in the real photon is large. They 
relate the normalization of this component of the real photon wave function to the quark 
condensate and its magnetic susceptibility [243]. 

4.4 The impact factor representation for the heUcity ampUtudes 

The /c^-factorization, or impact factor, representation for the vector meson production 
repeats closely that for the Compton scattering amplitude in the case of DIS [37,168]. 
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The three changes are that now the momentum transfer A 7^ 0, the vertex function for 
the S'-wave vector meson is different from the 7^ vertex (104) for the photon, and the con- 
ventional unintegrated gluon SF which describes the t-channel exchange is replaced by the 
off-forward (skewed) unintegrated gluon structure function of the target, T{x\^ X2-, i^i-, 1^2) ■ 
Here 



+ Ml 

iy2 



X2 



Ml - rr^y 



(114) 



and Ml is the invariant mass of the intermediate qq pair, for the kinematical variables see 
Fig. 16. 



V 



k-(1-z)A 




Figure 16: The kinematical variables entering the k± factorization representation 
(116) for vector meson production amplitudes. 



The imaginary part of the total amplitude can be written as [54, 171] 



Im r = p^2CyV47rae^ f d^n,^ ^_„„r..2 ^2, ,,2 



47r2 



dzd'^k 

a5(max{K'', e'^+k'^}) J'(xi, X2, Ki, ^2) / — —■I{Xv,X-y), 



zl-z) 



(115) 
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where for the pure S'-wave vector mesons the integrands /(Ay, A^) have the form 

;i - 2zf 2m/ 



r{L,L) = AQMz\l-z) 



V — '^f 



Y = — A^ 



I%L,T) 



1 + 



4:z{l - z) M + 2mf 



2^2 , 



+ K + (1 - ^) ](*r • ^1 



(116) 
(117) 



[{k-^*)^2-{2z-l)\k-^^)%] 



M + 2mf 

2z{l-z){^,,^l^-<^^y^l^) (118) 
[{Kn. - kyny)^2 - {2z - lf{k.A,., - ky^iyW^ ; 

1 - 2zf 2m/ 



m/ 



M + 2m/ 
-2Mz(l - z){2z - l)(e$i)^* 



1 + 



Az{l - z)M + 2mf 



Mm f , 
+ T^-7^(22- l)(e**)$2 



M + 2m/ 
-2Qz(l - ^)(2z- 1) 



^' M + 2m/^ . 



(119) 
(120) 



Here r = Ic + (2; — i) A and 



$2 



1 



+ 



1 



+ 



1 



+ 52 (r + A/2)2 + e2 (r - A/2)2 + e 



(121) 



_ r + K r-K r + A/2 r-A/2 

' " ~ (r + k)2 + £2 ~ _ k)2 + ^2 + + ^/2)2 + £2 + _ ^/2)2 + ^2 ' ^ ) 

for the definition of ^ see eq. (77). Here l/(if^ + £^) and ^2 = "^vi^jk) describe 
transitions into the q states with the sum of hehcities of the quark and antiquark A + A = 
X'y*,Xvy whereas k/{k'^ + e"^) and ^'i = kipv{z,k) describe transitions of transverse and 
vector meson into the q states with A + A = 0, in which the hehcity of the photon and 
vector meson is carried by the orbital angular momentum in the qq state. 

In the calculation of the double helicity-flip amplitude (118) the x-axis is chosen along the 
momentum transfer A. The point made in Section 3.4.4 that the helicity flip proceeds 
via the intermediate state with A + A = is manifest in (117)-(119). The corresponding 
integrands for the D-wave states can be found in [171]. 

The 2;-dependence of the integrands shows that the end-point contributions (2; ^ 1 or 
1 — 2; -C 1) are suppressed in the longitudinal amplitude T(L, L) already in the integrands, 
while for the other helicity amplitudes this suppression comes from the wave functions, 
see discussion in Section 4.7 and 4.9 below. The factor (22; — 1) in the integrands of 
the helicity-fiip amplitudes T(L, T) and T(T, L) corresponds to the longitudinal Fermi 
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momentum of quarks in the vector meson, which makes manifest the relativistic origin of 
hehcity flip. The expected hierarchy of the hehcity flip amplitudes is as follows [170,212]. 
Roughly, 



ITn 



Oil 



(123) 



2" ' y^'^^) 

"V 



(125) 





^1-1 1 





\t\ 






|2 + 


^0 


2 mv^/Q'^ + rriy 



For heavy flavour vector mesons, the helicity flip amplitudes are expected to be further 
suppressed by the non-relativistic Fermi motion. 

The real part of the amplitude can be reconstructed from the imaginary part using the 
derivative analyticity relation [244,245]: 

4.5 The off-forward unintegrated gluon density: the A-dependence 
within the diffraction cone and the BFKL Pomeron trajec- 
tory 

Thanks to a large amount of high-precision data on both in the soft and hard regimes, 
the simple, ready-to-use parameterizations for the forward unintegrated gluon density 
T{x^K?) are now available [34]. These parameterizations can be exploited in different 
high-energy reactions and bring the gluon density of the proton under control. 

For the practical application of the formalism of section 4.4, one needs the off-forward 
unintegrated gluon distribution JF(xi,X2,k + |A, — k + |A). Its dependence on the 
momentum transfer A comes from two courses. The first one is the soft quantity that 
can be dubbed the two-gluon form factor of the proton. The second is the A-dependence 
of the BFKL two-gluon ladder. When viewed in the impact parameter space, at each 
splitting of the gluon into two gluons, g — > gg, the hard gluon, which carries the large 
longitudinal momentum of the parent gluon, emerges at the same impact parameter as the 
parent gluon, whereas the soft one, which carries small longitudinal momentum, emerges 
at an impact parameter |Abj| ~ from the parent gluon. Consequently, as illustrated 
in Fig. 17, the splitting of gluons in the process of the log ^ evolution is accompanied by 
the Gribov-Feinberg-Chernavski random walk [246,247] of small-x gluons to larger and 
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larger impact parameters b. The asymptotic freedom, i.e., the running as, enhances the 
role of large random walks of the order of the perturbative gluon propagation radius Re- 
This suggests that for (b^) will rise proportionally to the number of gluon splittings, i.e., 

(b^) oc i?^log- oc i?^logiy^ (127) 




Figure 17: The sequential splitting of gluons in the Feynman diagram (a) for 
production of the multigluon final state viewed as a random walk in the impact 
parameter space of gluons from the qq pair of the photon 7* to the proton target. 

In conjunction with the Regge formula (37) and the definition (43), this entails the finite 
slope a'^pKL of Regge trajectory of the hard BFKL Pomeron. Evidently, the dimen- 
sionfull Ci'^FKL is the soft parameter and as such it depends manifestly on the infrared 
regularization of QCD. The solution of the color dipole BFKL equation with Yukawa-type 
cutoff and infrared freezing of as gave d'^p^L — 0-12 ^ 0.15 GeV~^ [25, 134,248]. The 
results for the shrinkage rate c^'^fkl depend on the admixture of subleading BFKL poles 
and exhibit weak dependence on k^. The quoted value is found for the specific boundary 
condition, which gives a good description of the HERA results on the proton structure 
function [134,136]. 

When viewed in the momentum space, the same Gribov-Feinberg-Chernavski diffusion 
suggests the weakening of the k-A correlation with the number of splittings. Indeed, 
Balitsky and Lipatov have shown that the dependence of the off-forward gluon density 
(k ■ A) corresponds to subleading singularities [37,127], and, arguably, can be neglected 
for HERA energy range. 
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Consequently, for small momentum transfers within the diffraction cone the A-dependence 
can be factored out as 

X2, hi + ^A, + ^A) = T{xi,X2, K, -k) exp ( ^^^^ j , (128) 



We parameterize h^p as 



hip = b2G + '^a'sFKL log ^2 I °2 ' (^29) 



where the soft parameter 62G can be regarded as a slope of the form factor of the proton as 
probed by the color singlet two-gluon state. In principle, one can determine it experimen- 
tally isolating the BFKL contribution to diffractive DIS into high mass states. Strictly 
speaking, this parameter 620 as well as the Pomeron slope a' can change from the soft, 
non-perturbative, to hard, BFKL, gluon density, taking [54,55] the universal parameters, 
62G = Bn = 4 GeV"^ with xo = 3.4 ■ 10""^ and the K^-independent a^'^* = o:'^fkl = 0-25 
GeV^^ is the poor man's approximation. 



4.6 The ofF-forward unintegrated gluon density: the dependence 
on skewness 

Bartels was the first to observe [249] that two gluons enter the amplitude at X2 7^ a;i ~ x, 
because the invariant mass squared Mf of the intermediate qq system is close to ^ rriy 
for the final qq state and is far from the virtuality of the incident photon —Q^. Such a 
skewed unintegrated gluon density can be, in principle, accessed in DVCS [31], but that 
is not yet a practical solution. Shuvaev et al. [250] and Radyushkin [251] argued that at 
small X the skewed distribution can be related to the conventional one: if oc x""^, then 

22A+3 r(A + ^) 

J^{xi,X2 <^ Xi,k.,-k) = Rg ■ J='{xi,K.) ; Rg = — , . (130) 

y/n i(A + 4j 

The above factor Rg can be effectively accounted for in a form of the x-rescaling 

iV / 1 ^ ^ 



where c(A) changes from ^ 0.435 at A = to 0.4 at A = 1. Given this very fiat dependence, 
one can take fixed c = 0.41, so that 

0, K, — k) jF(cx, k) . (132) 
Hereafter we approximate the skewed gluon density by the forward density take at 

Xg = cxi = c — — — . (133) 
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Of course, once this rescaling of x is implemented and the Fourier transform to the color 
dipole representation is performed staring from Eqs. (115)-(120), the color dipole and 
/c^-factorization approaches will be identical to each other. 



4.7 The Ansatze for the wave function 

For the heavy quarkonia a good insight into the functional form of the radial wave function 
(WF) ipvip^) comes from the potential model calculations [47,48]. Here, at least for 
the T(IS'), the role of the QCD Coulomb interaction is substantial. It is less so for the 
charmonium, whereas the gross properties of lighter vector mesons which have a large size 
are entirely controlled by the confining interaction and here one is bound to the model 
parameterizations [54] . The popular harmonic oscillator WF emphasizes the confinement 
property, it decreases steeply at large p^. 



il^is = Cl exp [^^-Y^ ' ^25 = C2 {inode " P^af) exp [^^-^^ ■ (134) 



which emphasizes the contrast between the non-pointlike vector meson and pointlike pho- 
ton for which r^*(2,k) = ^/Ana^ = const. The position of the node, ^node, is fixed 
from the orthogonality condition. The attractive pQCD Coulomb interaction between 
the quark and antiquark enhances the WF at small IZqq and/or large relative momentum, 
the minimal relativization of the familiar Coulomb WF suggests 



V'lS P = . 2 2^2 5 ^25 P = ^ . 2 2^3 ' ^35) 



which decreases as an inverse power of p^, much slower than (134). The factor 1/vM 
in (135) is a model-dependent suppression to make the decay constant (112) convergent. 
Arguably, those two extreme Ansatze give a good idea on the model dependence of vector 
meson production amplitudes. The radius Oi and the normalization ci are fixed by the 
Vo —>■ e^e~ decay constant (112) and the normalization condition (113). The hybrid 
model in which the short- distance QCD Coulomb interaction in light vector mesons has 
been treated perturbatively is found in [21,49]. 

4.8 The hard scale Q : the Unk to the leading log Q^-approximation 
and the exponent of the VF-dependence 

For soft gluons, 

^ < {e^ + k^) = z{l- z) (g2 + M^) , (136) 
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one can expand $2 and $1 as [252] (for the sake of simplicity we consider A = 0) 



$2 ; 
$1 



(£2 + k2)3 



4k 



1 - 
1 - 



2k2 



2(1 -^)(g2 + M2) 

k2 



K 



K, 



(137) 
(138) 



A natural approximation is M'^ ^ my. Then the factor (Q^ + 



m 



2 \-2 

v) 



which emerges in 



(137) and (138) corresponds to precisely the factor r| of the color dipole approach, see 
Eqs. (78) and (79). The determination of the hard scale in the gluon SF is a bit more 
subtle. 

2 

Expansions (137) and (138) define the leading logQ contribution with logarithmic inte- 
gration over k'^: 



2{1-^){Q2+A/2) 



J^ixg, k) = G{xg, z{l - z){Q^ + M^)) 



(139) 



The emerging running hard scale depends on z and M2, for the heavy quarkonia the wave 
function of the vector meson is peaked at 2 ~ ^ and one can take M"^ ^ my, consequently, 
z{l - z){Q^ + M2) = of Eq.(9). 

w(/cVw(o) 





w,(/c^)/w,(o) ; 




w,(^')M(o) - 


W=75 GeV, Q'=100 GeV' 









K, GeV" 



Figure 18: The normalized weight functions W^IQ"^, k'^)/Wl{0) and 
Wt{Q^,k'^)/Wt{0) for the p production calculated at Q2 ^ ^qO GeV2 in 
the fc^-factorization approach [38]. 



The contribution from small dipoles, r < rs, or from hard gluons beyond the leading 
logQ domain, >^Q , is an integral part of the /c^-factorization approach [21,33,34]. 
In this region 

2 2k 
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Figure 19: The scales Qgl gt which the p production maps the gluon density 
as a function of found in the /c^-factorization approach [38]. Shown also is the 
heavy flavor approximation Qq = Q = \{Q'^ + "^^). 

2 

and the correction to the leading logQ result (139) can be cast as (see also [173]) 

/7 2 
^^(x„ ^ Hx„q') ■ \ogC, . (141) 

where logC^ ~ 1 and depends on the exact K,^-dependence of J^{xg, k). Following [33,34] 
one can combine (139) and (141) as 

G{Xg,Q') + J^{Xg,Q') ■ \OgCg ^ G {x g , C gQ^ ) . (l42) 

Consequently, the gluon density is mapped at hard 

QI = CgQ' , (143) 

which is slightly different from Q . As already mentioned above, I^{L, L) of (116) is more 
peaked at z ~ |, whereas {T,T)xv=x^ of (117) extends more to the end points z ~ 

2 2 

and 2 ~ 1. This leads to an inequality Qqi > Qgt implies that the typical dipole 
sizes in the T ^ T amplitude are somewhat larger than for the L ^ L amplitude. 

For a more quantitative analysis the ^^-integrations can be cast in the form 

1 r ^ ~*2 

—ImTLL,TT ^ J -^H^g, ^) ■ WlAQ\ k') ^ WlAQ', 0)G(x„ Qgl.gt) (144) 

The typical behaviour of normalized weight functions Wl,t{Q'^, '?^)/W^l,t(Q^, 0) is shown 
in Fig. 18 and for smooth Ti^Xq, h?) they can be approximated by the step-function 

2 2 2 

^i.QGL,GT-'^'^)^ where Qgl,gt a^^e defined by the median, Wl,t{Q'^, Qgl,gt) = \Wl,t{Q'^, 0), 
the results are close to the ones found in [21]. At moderately large the strong scaling 
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violations in T{xg^ k^) shown in Fig. 10 have a strong impact on Q^l gt shown in 

2 2 2 

Fig. 19. The inequahty Qqi^ > Qgt found in [21] is retained and the hierarchy of Qgl gt 

2 

from hght to heavy flavors is the same as of Q . 

2 2 

This effect of Qgl > Qgt demonstrated in Fig. 20 on an example of G{xg, Qgt) ^^d 
G{xg,QQj^) for the p production: the finding of G{xg,QQi) > G{xg,QQj^) reflects the 
inequality Qgl Qgt equal Q^. For the same reason different helicity amplitudes can 
have a slightly different energy dependence. The effect of different scales diminishes at 
larger Q^ with weakening scaling violations, see Fig. 10. 




Figure 20: The effect of different hard scales on the integrated gluon density (142) 
which enters the dominant SCHC amplitudes T^l o-nd Ttt for the p production as 
found in the k±- factorization approach [55]. 



A tricky point is that the small-K expansion for the double-flip amplitude starts with 
the constant and it is dominated by the soft-gluon exchange, the term oc is of higher 
twist [170,212]. 

2 

Of course, the above separation into the leading logQ and hard regions is redundant 
when either the color dipole and the full fledged /c^-factorization approaches are used. 
The fc_|_-factorization analysis reveals somewhat better the role of in (136) for the 

2 

determination of the leading log Q region. One may call that the inclusion of the Fermi 
motion effects [239], it is automatically contained in the color dipole calculations and 
should not be discussed separately. The specific form of the Fermi motion correction 
suggested by Frankfurt et al. [239] is not borne out by the /c^-factorization analysis, 
though. 

Finally, Eq. (144) clearly shows that the energy dependence of the production amplitude 
is controlled by the x^-dependence of the integrated gluon density, i.e., by the effective 
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intercept ^(Qgl.gt) defined in section 3.3.1: 



ImT{W^,t 



0)ll,tt oc - 

X 



l+\{Qc 



(145) 



4.9 The production amplitude and the vector meson distribu- 
tion ampUtude 

At large there emerges one useful approximation usually referred to as the coUinear 
approximation. Specifically, here the explicit Ic^- dependence in Q is neglected and the 
d^k integration can be factored out. In the color dipole language this amounts to ne- 
glecting the r-dependence of the vector meson wave function and taking it at r = 0, but 
keeping its ^-dependence. The results are best seen in the momentum representation. For 
instance, the factor z^{l — z)"^ in I^{L,L), Eq. (117) cancels the factor z~^{l — z)~^ in 
eq. (137) for $2, and the vector meson wave function enters Tn in the form of an integral 



d^kMip*y{z,k) 



1 + 



:i-2zf 



2m 



f 



(146) 



4z(l-z) M + 2m/_ 

which is very similar to the distribution amplitude defined by (112) at the factorization 
2 

scale Q . We emphasize that the factor in the square brackets in (146) depends on 



the S'D-wave mixing. For the nonrelativistic heavy quarkonia M ^ my, 
proportionality of the two quantities is exact, and one finds, cf. (68), 
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(147) 



(148) 
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with Rlt = 1 if the slight possible difference of the t-dependence for the L and T cross 
sections is neglected. Since the proportionality Ry oc at small is a generic con- 
sequence of the electromagnetic gauge invariance, the true dynamical features of vector 
meson production are manifested by the departure of Rlt from unity. We strongly advo- 
cate to represent the experimental data in terms of this parameter. 



4.10 The ratio R = (Tl/ctt and short distance properties of vector 
mesons 

One word of caution on the ratio (148) is in order [253]. In the Introduction we emphasized 
how the vector meson production 7*p Vp is obtained by analytic continuation from 
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the elastic Compton scattering 7*p 'y*p. In the course of this analytic continuation one 
changes form the pointlike '~f*qq to the non-pointlike Vqq vertex. 

Making use of the optical theorem for the Compton scattering, one finds 

Here we used the prediction [135] for inclusive DIS Rdis = (^L{'y*p)/o'T{'y*p) ~ 0.2, 
which is consistent with the indirect experimental evaluations at HERA [254]. This result 
Rcompton <^ 1 for the elastic scattering of pointlike photons 7*p 7*p must be contrasted 
to Ry ~ Q^/rriy ^ 1 when the pointlike 7* in the final state is swapped for the non- 
pointlike vector meson. Evidently, the predictions for Ry are extremely sensitive to the 
presence in light vector mesons of quasi-pointlike qq component with the WF concentrated 
at short qq separation. The crude model estimates in [21,49] suggest that modifications 
of the wave function by attractive short-distance pQCD interaction do indeed lower the 
theoretical results for Ry- 

Here we just recall the pQCD radiative correction to the Weisskopf-Van Royen non- 
relativistic approximation (65) for the leptonic decay width ( [255], see also [256]) 

r(yO ^ e+e-) = ^^^\^^,y{0)\' (1 - ^asimf)]' , (150) 

which even for the J/\l/ suppresses the decay width by a factor of ~ 2. Remarkably, 
this particular NLO correction is of Abelian nature - it derives from the Karplus-Klein 
QED radiative correction by a substitution Cpasim'j). In the LO approach the 

conservative radii of vector mesons are fixed from Eq. (112) without allowance for the 
pQCD correction (150). Although the formula (150) can not be directly applied to light 
vector mesons, it is reasonable to wonder what will happen to the vector meson production 
phenomenology if the Celmaster & Barbieri et al. correction for the p meson were a factor 
of 3. To a crude approximation, that will enhance the wave function at the origin by the 
factor ^ -\/3 and decrease the radius of the vector meson by the factor ^ S^/^ ^ ^.44, 
which is not off-scale. We expect a substantial reduction of the predicted Ry = ctl/o't 
for such a squeezed p-meson. 

4.11 Perturbative QCD calculations at high-^ 

Vector meson production at large \t\ is believed to be dominated by small impact pa- 
rameters b ~ l/A/|t[. Simultaneously, the large \t\ is expected to select the small-size 
configurations in the 'y*V transition vertex. Consequently, at \t\ such that 1/ ^/\t\ <C r^. 
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i.e., |t| ^ Q , it becomes the hard pQCD scale for the process. Arguably, in the real 
photoproduction of heavy flavour mesons the correct hard scale is \t\ + rriy. 

Which \t\ is large enough for \t\ or |t| + rriy to become the hard pQCD scale? That can 
be decided only a posteriori, the answer depends on the normalization of the hard pQCD 
amplitude and on how fast the soft amplitudes do vanish at large A [46]. The large-t data 
taken at HERA correspond to the Regge regime of |t| <^ W"^. In the real photoproduction 
one starts with the typical hadronic 7 ^ transition, and our experience with large-t 
hadronic reactions is a very discouraging one. At small t within the diffraction cone the 
differential cross sections have the exp(— behaviour, but at larger t > 1 GeV^ the 
slowly decreasing multiple-pomeron exchanges take over: the n-pomeron exchange gives 
the t-dependence cx exp(— -^1^1) (for the review and references see [93]). It is fair to say 
that the unequivocal evidence for hard pQCD mechanism in high-energy elastic proton- 
proton scattering is as yet missing, the soft double-pomeron mechanism dominates for t of 
several GeV^ quite irrespective of the specific model for the soft-pomeron amplitude (for 
the recent fits to elastic pp scattering see [257]). The model-dependent estimates show 
that the rate of decrease of the soft amplitude slows down dramatically at large- For 
instance in tttt elastic scattering at moderate energies the dominance of the hard pQCD 
amplitude requires \t\ > 4 GeV^ [258]. The modern handbag mechanism for large-t two- 
body reactions, as well as electromagnetic form factors of nucleons and pions, relies on 
soft wave functions of hadrons ( [259] and references therein). 

Under these circumstances the single-BFKL pomeron exchange interpretation of the large- 
t vector meson data is at best the poor man's approximation. Under this very strong 
assumption of two-gluon tower exchange, the fc^-factorization formalism expounded in 
section 4.4 is perfectly applicable at large The approximation of section 4.5 for the 
unintegrated gluon density is only applicable within the diffraction cone, A^i?^ < 1, and 
must be modified. 

The available experimental data at large- |t| are for the proton dissociate photoproduction 
7P — s> V{A.)Y, which at large t can be described in the equivalent parton approximation 
of Ginzburg et al. [30], 

where x' = |t|/(my + |t|) is the fraction of proton's lightcone momentum carried by 
the struck parton. It is reminiscent of the familiar coUinear factorization, but in the 
calculation of the hard cross section dav{'^*q — * Vq')/dt the exchange by soft gluons with 
the momentum \k\ <^ |A| must not be included. For instance, in a splitting q q'g with 
the q' —g relative transverse momentum p the transverse size of the q'g pair is r„/„ ~ l/|p|- 
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The exchanged gluons with the wavelength A = 1/|k| ^ Vg/g can not resolve such a pair 
which will act as a pointlike color triplet state indistinguishable from the parent quark q. 

The first application of (151) by Ginzburg et al. was to the process 77 — > VY treated 
in the two-gluon exchange approximation. In the more advanced BFKL approach this 
constraint amounts to endowing the target partons in (151) by a t-dependent non-pointlike 
structure; the practical prescription has been developed by Forshaw and Ryskin ( [53], 
see also [260]). 

Notice that W^^ = x'W^p and the Regge parameter for 7*g — > Vq' is 

exp(Ary) = -^f— , (152) 
my — t 

where Arj is the rapidity gap between the produced vector meson and the hadronic de- 
bris of the proton. The lower limit of the x'-integration, Xmin < x' < 1, is set by the 
experimental cuts. 

Hereafter we focus on real photoproduction. Beyond the leading order pQCD the only 
working approximations for T{xi,X2, k + | A, — 1 A) in the large-A regime is based on 
the Lipatov's solution of the leading order BFKL equation in the scaling approximation 
as = const. In contrast to the cases of DIS or diffractive vector mesons at small A, 
where the BFKL evolution from the proton side starts from the soft scale an becomes 
hard only on the virtual photon end of the gluon ladder, at large \t\ the large momentum 
transfer A flows along the whole ladder which may make Lipatov's scaling approximation 
better applicable at large \t\. Important point is that at large A the A. — k, correlation 
in jF(xi,X2, K. + 1A, —k + |A) becomes very important, technically the dependence on 
the azimuthal angle between A and k is described by the conformal spin expansion. 
We wouldn't go into the technicalities of the formalism, it is fair to say that for the 
leading helicity amplitudes the changes from pQCD two-gluon exchange to scaling BFKL 
approximation are for the most part marginal. The sensitivity to the wave function of 
the vector meson and real photon is dramatic, though. The detailed discussion is found 
in the recent paper by Poludniowski et al. [261], here we summarize the major points. 

The equivalent parton representation (151) makes it obvious that the helicity properties of 
the target dissociative reaction do not depend on the target. In their analysis of 77 — > VY 
to the perturbative two-gluon exchange approximation Ginzburg et al. [30] allowed for the 
Fermi motion of quarks parameterized in terms of vector meson distribution amplitudes 
with broad z-distribution and found the dominance by the spin-flip transition 77- —>■ Vl, 
cf. the large-t extension of eq. (123). The double- flip and non-flip amplitudes have similar 
t-dependence and are suppressed. The found differential cross section is of the form 

'j^lhllf^o.a,,^ra,r(V^e^e-)$. (153) 

at \tr 
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Figure 21: The helicity amplitudes differential in z = u for pQCD two-gluon and 
scaling BFKL approximations for \t\ = 10 GeV^ and asAy ~ 2.4, where Ay is the 
rapidity gap between the vector meson and debris of the proton. The pQCD two- 
gluon results have been multiplied by a factor of 3. From Poludniowski et al. [261]. 

The flavour dependence is typical for tlie cxL-dominance, cf. eqs. (147), (174). Tlie sensi- 
tivity to tlie wave function is obvious from the much discussed unreahstic nonrelativistic 
limit of z = ^. Here the hehcity flip amplitude vanishes and , see eq. (119), the SCHC 
transition dominates ( [262]), also see eq. (119),, the differential cross section (153) ac- 
quires the extra factor my/\t\ and, furthermore, the predicted cross section exhibits an 
accidental, artificial, zero at \t\ = rriy. 

Upon the radiative corrections the electromagnetic vertex acquires the anomalous mag- 
netic moment (Pauli) component oc a^yp1/2mf which, as we discussed in section 4.3.3, 
contributes to the chiral-odd parallel-helicity wave function of the photon. Schwinger's 
classic calculations show that for virtual photons the form factor of such a perturbative 
Pauli vertex vanishes oc rn^/Q'^ (see §117 of the textbook [263]). Ivanov et al argued [242] 
that for real photons the nonperturbative chiral-odd vertex is substantial - they relate it to 
the product of the quark condensate and the magnetic susceptibility of the vacuum [243] 
- and will enhance strongly the SCHC transition 7^ — Vt- 

The ideas of D. Ivanov et al. [242] have been extended to the scaling BFKL approximation 
by Poludniowski et al. [261], where one can find references to early studies. The crucial 
point is an enhancement of the SCHC amplitude by the chiral-odd parallel-helicity compo- 
nent in th photon. The importance of realistic ^-distributions in vector mesons is clearly 
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seen from Fig. 21 which show the hehcity amphtudes for the p production differential 
in z for chiral-even and chiral-odd photon wave function. As we stated above, for the 
leading chiral-even helicity amplitude T^i (M+o in the notations of Poludniowski et al.) 
the changes form the pQCD two-gluon to scaling BFKL approximation are marginal. The 
BFKL approximation enhances further the chiral-odd contribution to the SCHC non-flip 
amplitude and makes an approximate SCHC the dominant feature of large-t vector me- 
son production in the studied region of |t| < 6 GeV^. The double-flip amplitude also is 
enhanced. 

Finally, making use of the scaling BFKL unintegrated glue one is committed to the BFKL 
intercept (59) and prediction of the steep rise of the cross section with energy 



which has been emphasized by Ginzburg et al. already in 1986 [30]. 
4.12 Beyond the leading log^ approximation 

The above described pQCD description of the vector meson production is based on the 
manifestly leading log^ BFKL formalism. Going to the NLO log ^ BFKL remains the 
major challenge to the theory. The principal feature of the leading log^ approximation is 
that adding soft perturbative gluons, i.e., the higher Fock states, to the color dipole can 
be reabsorbed into the x-dependence of the color dipole cross section, which is equally 
true for DIS and vector meson production. Going to the NLO log^ approximation is much 
more tricky. While the nearly decade long efforts have culminated in the derivation [264] 
of the NLO BFKL evolution kernel, the matching calculations of the effect of hard gluons 
in the NLO impact factors are missing. In the case of DIS those hard gluons are of 
perturbative origin, but even so, despite the great progress [265,266], the closed result 
for the impact factor of the virtual photon is not available yet. In the case of vector 
mesons, the evaluation of the NLO impact factor can not be separated from the issue 
of the higher, nonperturbative, qqg, Fock state of the vector meson, in which the gluon 
carries a finite fraction of the vector meson's momentum. With the reference to the non- 
relativistic intuition, one may argue that an admixture of such a non-perturbative qqg 
Fock state is small in the T, but not the lighter quarkonia the non-relativistic treatment 
of which is suspect. The issue of nonperturbative qqg and higher Fock states in light 
vector mesons remains open. Hence we are bound to stay within the qq Fock state and 
leading log^ approximations. 

None of the NLO corrections is expected to change substantially the predicted Xg- and 
Q^-dependences but they affect strongly the predicted cross section. For instance, one 




(154) 
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often includes the correction for skewness (132) for the reason that at small x and large-Q 
it enhances the predicted cross section almost by a factor 2. Similarly, the real part of the 

2 

dipole amplitude (126) is a NLO correction which is substantial at large Q . A consistent 
treatment of the potentially more important Celmaster & Barbieri et al correction (150) 
and its counterpart for the impact factor is not available yet. In their analysis of NLO 
corrections Levin et al. did not consider the Celmaster & Barbieri et al correction but 
argued that by analogy with the Drell-Yan production of lepton pairs the cross section of 
near threshold diffractive cc pairs acquires the i^-factor [267] 

K^l + ya5, (155) 

which, by virtue of duality arguments, shall propagate to the vector meson production 
cross section, while Dremin asserts that NLO Sudakov effects rather suppress the cross 
section [268]. The full fledged NLO log^ fc^-factorization analysis necessary for consistent 
treatment of all these corrections is not yet available. 
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5 Helicity properties of vector meson production 



5.1 General introduction 

The angular distribution of the exclusive production of vector mesons decaying into 
particle-antiparticle final state is usually described in the so-called helicity frame [42,269]. 
For the unpolarized lepton beam, such as that at HERA, the cross section depends on 
three angles explained in Fig. 22 for the specific case of p° tt^tt" production. 




Figure 22: Kinematics of the process ep — > ep^p — > e7r~^7r~p. The angular 
variables used in the determination of the helicity amplitudes from the decay distri- 
butions: the azimuthal angle between the (e,e') scattering plane and the {p^,p') 
production plane; (ph, the azimuthal angle between the production and decay planes; 
and 6h, the polar angle of the positively- charged decay pion defined with respect to 
the direction of the p° momentum in the Yp cms system. Illustration is taken from 
[15]. 

This angular distribution can be expressed in terms of the virtual photon density matrix 
P\'^\.y{^h) (22), the helicity amplitudes of virtual photon transition with helicity into 
the vector meson with helicity Xy, Txv,\^i and of the angular factors that describe the 
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vector meson decay into the final particle-antiparticle state. For the case of decay into 
scalars (p tt+tt" and — K^K~ decays), these angular factors are given by spherical 



harmonics ^i.a^: 
da 

dcosOh d(j)h d^h 



a-W{cos9hAh,^h) (156) 

X~f,X!^;XvX.y 

In the case of the J/\l/ production one measures the angular distribution of leptons in the 
decay J/ip i'^i", then the D^/^(6'/i, 0/i)-functions instead of spherical harmonics will 
appear in (157). Since the photon density matrix Pa;:^At. and the angular factors are known, 
the study of the angular dependence of the cross section reveals the helicity structure of 
the 7*p — > Vp transition. 

The conservation of the s-channel helicity in the scattering of electrons in the Coulomb 
field is known since 1954 ( [209], see also the textbook [210]). Motivated by the early ex- 
perimental data on vector meson photo- and electroproduction. Oilman et al. suggested 
the s-channel helicity conservation (SCHC) as the fundamental feature of diffraction scat- 
tering ( [269], for the discussion of the pre-HERA experimental situation see [42]). Within 
diffraction cone, the helicity of the vector meson Xy coincides approximately with the he- 
licity of the incident photon X^, see Section 3.4.6. It is reasonable therefore to start with 
the case of strict s-channel helicity conservation (SCHC), Xy = X^. 

As emphasized in Section 3.1.2, the Pomeron and all higher lying secondary reggeons are 
natural parity t-channel exchanges. Hereafter we analyze the helicity properties of vector 
meson production assuming natural parity exchange. A good idea on why the possible 
contribution from unnatural parity exchange can be neglected is given by the longitudinal 
double-spin asymmetry. 



5.2 Longitudinal double-spin asymmetry and unnatural parity 
exchange 

One is familiar with the helicity structure function of the proton, gi{x,Q^), which mea- 
sures the mean helicity of partons in the longitudinally polarized proton. It is determined 
experimentally in the DIS of the longitudinally polarized leptons off longitudinally polar- 
ized proton target, where the polarized leptons serve as the source of circularly polarized 
photons. The measured cross section is proportional to the imaginary part of the helicity 
conserving forward Compton scattering amplitude, 

Tx.yXN,x-yXN Q^) + X^XNQiix, Q^) oc 1 + A^l^j^^'' X^iXn , (157) 
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where X^,Xn are the hehcities of the photon and target nucleon [X^ = ±1). Such a 
hehcity dependence for transverse photons emerges naturally for the axial- vector {Ai, ...) 
meson exchange. The term oc gi{x, Q^) changes the sign when the circular polarization of 
the photon is flipped and, from the t-channel exchange point of view, corresponds to the 
unnatural parity exchange, see Eq. (88). As explained in Section 1.3, the vector meson 
production amplitude derives from the Compton amplitude by analytic continuation in 
the virtuality of the photon to the vector meson pole, which should not change dramati- 
cally the asymmetry parameter in the amplitude. Consequently, if one parameterizes the 
unnatural parity exchange into the transverse vector meson production amplitude as 

^v-Aiv,A^Ajv|A^=A^ OC 1 + -A^^X^^Xn , (158) 

then the natural expectation for the transverse cross section will be ax is [270] 

Al, ^ 2Afi'^ . (159) 

There is a purposeful difference between expansions (157) and (158) because in the vector 
meson production one measures the differential cross section oc In ai-ux unseparated 
vector meson production the asymmetry is diluted for the presence of ctl, 

AX^J^^J^. (160) 
^ l + Rv l + Rv 

(here we are back to the usual notation Ai = A^j^^^^) hence the dilution factor (l + Ry) in 
(160) compared to (159). As a matter of fact, such a relationship between the longitudinal 
double-spin asymmetries for DIS and vector meson production has been suggested already 
in 1976 by Fraas on the basis of the vector dominance model [271]. 

The results of the first experimental determination of A^ for the diffractive p production in 
the HERMES experiment [272] are shown in Fig. (23) in comparison with the estimates 
from the DIS data based on Eq.(160) but without the dilution of AY by the factor (l+Rv). 
A summary of the high precision experimental data on the ratio of polarized, helicity, 
gi{x,Q'^), to unpolarized, Fi{x,Q'^), proton structure function is shown in Fig. 23, the 
important point is that the effects of unnatural parity exchange vanish at small x, for 
pQCD arguments in favor of that see [273-275]. 

5.3 The angular distribution in the SCHC approximation 

5.3.1 Theoretical expectation: angular distributions 

Within SCHC and natural parity exchange we are left with two independent helicity 
amplitudes, Tu and Tqq. Then the angular dependence of reaction 7*p Vp, with V 
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Figure 23: Upper two plots show the x-dependence of the longitudinal double- 
spin asymmetry in exclusive p° meson electroproduction on the proton (top 
box) and deuteron (bottom box). The data from HERMES [272] are compared to 
the expectations from Eq. (160) based on the DIS data from the bottom plot, which 
shows the ratio of polarized to unpolarized proton structure function from the SMC, 
EI43 and HERMES experiments. For the review and references see [276]. 
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decaying into two scalars, is: 
W {cos Oh, (ph,'^h) = 



1 3 
iV4^ 



e|7^oPcos^ Oh + ^l^iil^sin^ Oh 



1 



+ -e\Tu\^ sin^ ehCos2{<l>h - (j)h) 

- ^/2e{l + €)Re{Tn%l) sin Oh cos Oh cos($^ 



,(161) 



where = |7iip+e|7ooP and e is defined in (22). The first hne in (161) is the contribution 
of diagonal terms of the photon density matrix, i.e. with = A^^; the second line is the 
interference of transverse photons with opposite hehcities, A^ = —X'^ = ±1; the last line 
is the interference between the transverse and longitudinal photons. Note that angular 
dependence (161) involves only a single azimuthal angle ip = ^h — (ph between the (e, e') 
and decay planes. 

One can reparameterize, and analyze, the angular distribution (161) in terms r°-, which 
compose the spin-density matrix of a vector meson [67]. The generic case and the involved 
notations for the r°- are explained in Section 5.3.1, here we only notice that under the 
assumption of SCHC many the elements of r^- do vanish and (161) takes the form: 



W {cos Oh, (l)h,'^h) 



3 



00 J 



(3r°^-l)cos2^;, 



+ e cos 2$ftA/2 sin^ 6h cos 2(f)h ■ r\_^ 

— e sin 2$/i sin^ 6^ sin 20/i ■ Im{r^_]^} 

- \/ 2e(l + e) cos sin 26h cos (ph ■ V2Re{riQ} 



+ y2e(l + e) sin sin 2^^ sin 20;, ■ y2Im{r^o} • (162) 

It contains five non-zero spin-density matrix elements among which only three are inde- 
pendent due to the relations (88). Their expression via the helicity amplitudes reads 



,04 
00 



1-1 



irnP + elTooP 



-Im{r^_J 





l^lll 
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l^lll 


2 + e|Too 


2 



(163) 



Re{rl} = -Im{r?J 



RelTnToU 



2v^|rn|2 + e|TooP 
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Under the SCHC conservation, the ratio of the longitudinal to transverse cross sections, 

00' 



(26), is expressed only via matrix element 



1 ^04 

e 1 - ^00 



while the relative phase 6 between the Tu and Too amplitudes can be determined via 

cos 5 = (Re{r^J - Im{r^J) . (165) 

V nv/2 



5.3.2 Experimental results 

Figure 24 shows the results of ZEUS and HI on the five "non-zero-SCHC" (hereafter just 
"SCHC") matrix elements mentioned above. These matrix elements are placed in three 
rows. The first line corresponds to diagonal terms; the second line is the interference 
of transverse photons with opposite helicities; the last line is the interference between 
transverse and longitudinal photons. 

The matrix element Tqq is extracted from the single-differential cross section da /d cos Oh- 
The photoproduction measurements [68,71,277,278] confirm that at = the matrix 
element Tqq is zero within experimental uncertainties for p, and J/ip mesons. This 
should be expected, since in the limit ^ the longitudinal cross section must vanish. 
The behavior of Tqq as function of Q^, W and t for p mesons is shown in Figs. 24, 25. 

The steep dependence of Tqq is mainly due to the gauge invariance driven factor 
~ Q"^ /my present in the longitudinal cross section, see (81), (82). The pattern of the 
dependence is very similar for all the light vector mesons and differs in the case of 
the J/ip (not shown). The matrix element Tqq is the main source of the determination of 
the ratio Ry = cri/aT, whose Q^-behavior will be discussed in detail in the forthcoming 
Section 6.3. The t and W dependences of Tqq are consistent with being flat. This indicates 
that the energy dependence of the longitudinal, Tqq, and transverse, 7^i, amplitudes is 
very close to each other in agreement with theoretical expectations. The same holds also 
for the t-dependence of the longitudinal and transverse amplitudes. 

The Q^-behavior of the other SCHC matrix elements can be read off from Eq.(163). The 
elements r\_^ and Im{rf_^} should approach ±1/2, respectively, in the photoproduction 
limit, and are expected to decrease with increase approximately as \/{2R). This 
tendency is well observed in the data. The longitudinal-transverse (LT) interference driven 
elements Re{rfQ} and Im{r^Q} should be oc Q at small and fall off with growth 
as l/Q. The experimental data from the ZEUS collaboration do follow this expectation, 
the experimental data from HI on r^^ exhibit certain departure from the theoretical 
expectation. 
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Figure 24: The compilation of the ZEUS [15] and HI [16] results on the 
behavior of the five s-channel helicity conserving amplitudes for p meson. The first 
matrix element corresponds to the diagonal terms in photon density matrix; the 
second line is the interference of transverse photons with opposite helicities; the 
last line is the interference between transverse and longitudinal photons. The lines 
show the kffactorization predictions [54, 55]. 
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Figure 25: Upper plot. The behavior of the spin density matrix element Tqq 
for p production measured by ZEUS [83]. The solid curve is a fit of the form 
= 1/(1 + ^(m2/g2)«) yj^th ^ = 2.17 ± 0.07, k = 0.75 ± 0.3. Bottom three plots 

show the spin-density matrix element Tqq for p meson electroproduction as function 
of Q"^, W and t measured by HI and ZEUS ( [85] and references therein). 
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Figure 26: The upper plot shows the theoretical expectation from the k±- 
factorization [54, 55] for the cosine of the relative phase 5 between the transverse, 
Til, and the longitudinal, %q, amplitudes. The bottom plots show the HI determi- 
nation of cos6 as a function of Q"^ , W and t [16]. 
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As we discussed in Section 4.8.1, because of the slightly different hard scales Q^^ and Qj- 
the two helicity amplitudes Tu and 7^o can have a slightly different energy dependence, 
see Sections 3.3.1-3.3.3 and the discussion of the experimental data in Section 7. Then the 
derivative analyticity (126) would predict the nonvanishing, but very small, in the range 
of 1 -T-3°, relative phase 5 of the two amplitudes. The results of the HI determination [16] 
of cos 5 from Eq. (165) are shown in Fig. 26. The combined value of cos 5 suggests the 
statistically significant departure from 1, 

cos 6 = 0.925 ± 0.022^°;°^^ , (166) 

the sizeable relative phase 5 ~ 10 ^ 20° is difficult to accommodate in theory. The error 
bars are large, though, it must be noted that ZEUS assumes 6 = and still obtains good 
fits to the angular distribution. 



5.4 Angular dependence beyond SCHC 
5.4.1 Theoretical expectations 

Without the assumption of s-channel helicity conservation, and supposing the natural 
parity exchange, see Section 5.2.3, one is left with five independent helicity amplitudes. 
The angular dependence of the cross section (for spinless particle-antiparticle final state 
such as in the decay p tt^tt^, K^K~) is parameterized in terms of 15 spin density 
matrix elements [67]: 

1 dd 



a dcos9h d(f)h d^^ 
3 



^ (1 - ^o°o) + ^(3rS^ - 1) cos^ Oh - V2Re {r%} sin 2^, cos 0, - r?li sin^ cos : 



— e cos 2$/i (^rl^ sin^ 9h + r^Q cos^ 9h — V2Re {rjg} sin 29h cos (ph — rl^i sin^ 9h cos 2(ph^ 
— e sin ^v^Imjr^Q} sin 29^ sin (ph + Im {rl_i} sin^ 9h sin 2(phj 

1 + e) cos [t^u sin^ 9^ + Tqq cos^ 9^ — -\/2Re {^^q} sin 29 ^ cos (ph — t^I-i sin^ 9^ cos 2(ph 
+ V2e(l + e) sin<l>, (V2Im {r?o} sin2^;, sin 0;, + Im {r^i} sin^ sin20/,) . (167) 

One can develop an intuition when reading this expression. The subscripts i, k = —1,0,1 
of the matrix elements indicate the vector meson helicities Ay and Ay of the ampli- 
tudes interfering. Dependence on $/i shows the helicities of the photon of the interfering 
amplitudes. As before, the ^^-independent terms are diagonal in photon helicities; they 
are accompanied by superscript 04. Terms oc cos2$/i and oc sin2$/i originate from the 
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interference of the transverse photons with opposite hehcities; these are accompanied by 
superscripts 1 and 2. Terms oc cos$/i and oc sin$/j appear from the interference of the 
transverse and longitudinal photons; these are accompanied by superscripts 5 and 6. The 
nature of each can then be understood from its indices. For example, rj^ comes from 
interference of photons of helicities +1 and —1 and from vector mesons with both helicities 
equal to 1 or —1. Therefore, rl^ must be proportional to the double spin-flip amplitude. 

The correspondence between the matrix elements r^, the helicity amplitudes, and their 
SCHC and SCHNC properties is as follows: 



SCHC 



^04 
00 



1-1 



e\%o\^ + 
- l|T |2 



IT 12 
■iri 



Oil 
1 IT |2 

9 -!1-1 



Im{rti} = -||Tnp + i|Ti_i|2 



'168) 



strong SCHNC 
from single-flip 

7t Vl, oc %i 



Re{r%} = |Re{TnTo*} + eRe{T,o%l} + '^Re{T^_^%U} 
ReMo} = iRe{rnTo*_ J + lRe{T,^^%\} 
lm{rl,} = -lRe{Tn%U} + ^^Re{T^^^T,\} 



' 00 



V2Re{Too%\} 



(169) 



weak SCHNC 
from single-flip 

7l -> Vt, oc Tio 



rl = j=Re{r,,T,\} - j=Re{T,oT,U} 
lm{rU} = -^Re{T_,,T,\} + j.Re{T,oT*,,} 



(170) 



double SCHNC r^^t^ = -e\Tio\^ + Re{TnT*_^} 

from the double- , _ , / » n 

fl. , ^ 1 = M^i-i^ii) 171 

nip, or two 

single-flips ^i^ ^ _|7-^^|2 

It is understood that the rhs of each and every line must be also divided by e(|7oop + 

2|rion + |rnp + |ri_ip + |roi|2. 
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The classification of eacli set of density matrix elements (168)-(171) in the left column of 
this Table can be understood as follows: At large Q^, the largest SCHNC amplitude is Tqi 
(the transverse photon to longitudinal vector meson transition). Therefore, the stronger 
SCHC violation is expected in density matrix elements of the second group, in particular, 
in Tqq. At small Q^, the double-flip Ti_i is predicted to be the largest SCHNC amplitude 
as it is soft-dominated and does not require the longitudinal Fermi motion. Therefore, at 
small the largest non-SCHC spin density matrix elements is expected to be rf^-^ and 
Too- The amplitude Tiq vanishes in the photoproduction limit, is of higher twist at large 
and is expected to be always small. 

As discussed in Section 3.4.6, the small-t behavior of various amplitudes is governed solely 
by the value of the helicity flip: the matrix elements from the second, Eq. (169) and third, 
Eq. (170), groups are expected to be cx: the elements in the last group, Eq. (170), 

are cx \t'\. 



5.4.2 Experimental results: helicity properties at small t 

In experiment, the density matrix elements are obtained by minimizing the difference 
between the three-dimensional {cosOh, (ph, $/i) angular distributions of the data and those 
of the simulated events. Figure 27 shows combined ZEUS [15] and HI [16] results of this 
procedure for the p meson production. Results for (j) mesons are presented in Fig. 28. 

All the matrix elements are placed in six rows. The first three rows, are the same as 
in Fig. 24 and show the spin density matrix elements coming from SCHC transitions. 
The last three rows represent the SCHNC matrix elements, which would vanish in the 
case of strict SCHC, for the definitions and classification of SCHNC as strong, weak and 
double see Eqs. (168)-(171) and discussion in Section 5.3.1 . The fourth row shows the 
four matrix amplitudes proportional to the single-flip helicity amplitude Tqi, the fifth 
row shows the matrix elements that are proportional to the single-flip helicity amplitude 
TiQ, the last row shows matrix elements with two helicity flips, either as a double-flip 
amplitude or as a product (square) of two single-flip amplitudes. 

Among the SCHNC spin density matrix elements the best measured ones are the com- 
binations 2r|]^ + Tqq and 2rj'^ + Tqq, which describe the cross section single-differential in 
angle $/j after the integration over cosOh. 

cx 1 + v/2^(rT^cos<l>,,(ro'o + 2ri\) - e cos 2$;,(r^o + 2rl^) . (172) 

Figure 29 shows also the resent HI measurement of the Q^- dependence of these matrix 
elements for the case of p mesons [85]. The first combination is significantly non-zero. 
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Figure 27: The spin-density matrix elements measured in the reaction 7*p — > pp as 
a function ofQ^. The first three rows are the same as in Fig. 24 and show the SCHC 
matrix elements. The solid symbols present the ZEUS [15] and the open circules 
the HI [16] results. The curves represent the kt- factorization calculations [54, 55]. 
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Figure 28: The spin-density matrix elements measured in the reaction '-f*p —>■ (pp 
as a function of Q^. The notation is the same as in Fig. 27. The curves represent 
the kt-factorization calculations [54, 55]. 
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Figure 29: The HI and ZEUS measurement ( [85] and references therein) of the 
combinations of SCHNC spin density matrix elements Vqq + 2rfi and Tqq + 2rli as 
a function of Q^, W , and t. 
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Figure 30: The HI and ZEUS measurement ( [82] and references therein) of the 
t' -dependence of combinations of SCHNC spin density matrix elements Tqq + 2r\^ 
and Tqq + 2r\i compared with of the behavior of the SCHC matrix element Tqq. The 
solid curves are from pQCD calculations by Ivanov and Kirschner [212]. 



87 



while the second is compatible with zero. The results of the direct fit to the entire angular 
distribution shown in Fig. 27 indicate that it is large non-zero value of Tqq which is a source 
of departure the SCHC in the first combination. Since this matrix element is proportional 
to the Re(7oo'7^i), it indicates that the single- flip amplitude Tqi is not zero. The data 
shown in Fig. 29 give an estimate of the relative strength of this amplitude at an average 
value of \t\ ^ 0.15 GeV^ 



The dominance of Tqq among the helicity violating amplitudes at large is in agreement 
with theoretical expectations [54, 170,212]. All the other matrix elements do not differ 
significantly from zero. There is, however, some indication that the double-flip amplitude 
is non-zero at small Q^, see the last row in Fig. 27. 

The t-dependence of the s-channel helicity violating amplitudes has already been shown 
in Fig. 29 within the diffraction cone. In Fig. 30 the data [82] on helicity violating matrix 
elements 2r{^ + Tqq and 2rf^ + Tqq within extended t-region, \t'\ < 2 GeV^ are shown. 
The behavior of both of these combinations is compatible with oc \/\t^ and confirms the 
theoretical expectations. For comparison, in the same Figure we show the t-behavior of 
the helicity conserving matrix element Tqq within the same region, which is compatible 
with constant. 

5.4.3 Experimental results: helicity properties at large t 

At large-t only the proton dissociative reaction 7p — > VY has been studied. The experi- 
mental data on the spin-density matrix elements from the recent ZEUS measurement [77] 
of the large-t p, cj) and J /if) photoproduction are shown in Figs. 31, 32, 34 in combination 
with the lower-t data [76]; the HI results [79] are shown in Fig. 33. 

The dominant feature of the DIS data within the diffraction cone, |t| < 1 GeV^ was large 
Tqq driven by the dominant (Jl from the transition 7l — V^. Such a transition is absent in 
real photoproduction. Figures 31, 32 show that the matrix element rj'g is positive valued, 
definitely nonvanishing one. The results for Tqq show that the probability to produce 
longitudinal p or mesons from a transverse photon increases with \t\ up to 4 9%, 
but is nowhere close to the SCHNC dominance of ai expected in the pQCD two-gluon 
model of Ginzburg et al. [30] with the chiral-even photon-quark-antiquark vertex. The 
matrix element r^^^ is numerically substantial and gives a solid evidence for a double-flip 
contribution. All the spin density matrix elements shown in Figs. 31,32 are small, so that 
the SCHC is the empirical feature of large-t production. 
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Figure 31: The spin density matrix elements for proton-dissociative p production 
(the top panel) and production (the bottom panel) as a function of —t measured 
by ZEUS [77] (black sumbols) and [76] (open symbols). The SCHC prediction is 
shown as the dashed line. 
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Within the theoretical approaches reviewed in Section 4.10 the sole source of approximate 
SCHC is the chiral-odd photon-quark-antiquark vertex. The theoretical calculations [261] 
are shown in Figs. 33 and 34. The pQCD two-gluon approximation grossly overpredicts 
Too for all vector mesons. The BFKL model correctly reproduces the gross features of 
Tqq and r^ti, but predicts the wrong sign of r°Q. The shown theoretical results are for 
the so-called asymptotic vector meson distribution amplitude, but the sensitivity to the 
model wave function is weak and can not explain this sign conflict with the experiment. 

For the heavy quarks in the J/\E' the Fermi motion is slow and the BFKL calculations 
predict much smaller helicity flip effects compared to the results for the light vector 
mesons, cf. Fig. 33 and 34. 
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Figure 32: The spin density matrix elements for proton-dissociative J/ip produc- 
tion as a function of — t measured by ZEUS [77]. The SCHC prediction is shown 
as the dashed line. 
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Figure 33: The spin density matrix elements for proton- dissociative J/V' produc- 
tion as a function of —t measured by HI [79]. The predictions form the BFKL 
(solid curves) and pQCD two-gluon approximations (dashed curves) are shown for 
a comparison [261] . 




Figure 34: The spin density matrix elements for proton-dissociative p (left box) 
and (f) (right box) production as a function of —t measured by ZEUS [77]. The 
predictions for from the BFKL (solid curves) and pQCD two-gluon approximations 
(dashed curves) are shown for a comparison [261]. 



91 



6 The dependence 



6.1 Total cross section 7*j9 Vp 

6.1.1 Theoretical expectations: what is the correct hard scale 
for the dependence? 

As emphasized in Introduction, see Section 1.3, the vector meson production at a given 
probes the gluon content of a target at pQCD factorization scale Q ~ (Q^ + my) /A, see 
(9). Theoretical predictions of the Q^-dependence of the vector meson production cross 
section in the color dipole approach were described in detail in Sections 3.4.2 and 4.9. 
Besides, Section 4.8 shows how the leading log-Q^ approximation results are recovered, 
and improved upon, in the more consistent /c^-factorization approach at large Q^. 

The theoretical expectation for the shape of the Q^-dependence of the transverse ampli- 
tude can be read off Eqs.(78),(79) and (147). Performing the t-integration in (147) and 
combining ax and ai, one finds the non-relativistic master formula 

2 

Although the absolute normalization, the exact value of the pQCD scale Q and the 
departure of Rlt from the non-relativistic quark model expectation Rlt = 1 do depend 
on the wave function of the vector meson, and the contribution from the helicity-flip 
transitions has not been included, Eq. (174) contains all the ingredients of the full pQCD 
description. 

One can present ay either as ay = crri)- + Rv) and use the theoretical approximations 
for aT or substitute (Tt = a^/Ry and test the theoretical predictions for ai. Brodsky et 
al. [23] argued that the end-point, z ~ 0, z ~ 1 contributions are minimal in, and the 
pQCD evaluations are more reliable for, the ai and our representation (174) corresponds 
to the latter, preferred, choice. 

The longitudinal cross section is small at Q"^ "^my, but becomes the dominant feature at 
larger and tames further the decrease of the total cross section at large Q^. In view of 
Rlt ~ 1 at small Q^, and because the My term can be neglected at large Q^, the factor 

(^Q^ + '^^^ ^he master formula is not very different from (Q^ + niy). Then the Q^- 

2 

behavior of the total cross section would have been ~ (Q )~"v with the exponent ny = 3 
modulated by the gluon density squared. As we shall see below, Rlt < 1 and steadily 
decreases with Q^, so that even without the scaling violations the expected exponent 
ny < 3. 
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Finally, the dependence of the diffraction slope can not be ignored. As we argued 
in Section 3.4.2, the diffraction slope hv{Q ) decreases with growing (Q^ + m^), which 
also enhances slightly ay at large Q^. As a matter of fact, the diffraction slope too is a 
function of Q rather than Q^. 

To summarize, there are strong theoretical reasons [49, 279] for the presentation of the 
experimental data as a function of either the scanning radius or (Q^ + my). In such 
a presentation the major flavor dependence is in the explicit factor myT{y — ^ e^e~). 
There is also a hidden dependence of the absolute normalization, Rlt and of the exact 

2 

dependence of the hard scale Q on (Q^ + my) all of which depend on the wave function 
of vector mesons [21,38]. 



6.1.2 Theoretical expectations: the impact of Xg-dependent scal- 
ing violations on the dependence 

The predicted Q^-dependence is driven mainly by two phenomena: the shrinkage of the 
photon light-cone wave function with and the resulting decrease of the scanning radius, 
and the non-trivial Q^-dependence of the gluon density at fixed W . The former property 

2 

leads to the strong decrease of the transverse cross section ~ (Q at a sufficiently large 

2 2 2 2 

Q . The gluon density, G{xg,QQ), at a fixed W, depends on Q also via Xg ^ Q /W'^. 
At small to moderately large Q^, one observes the quick rise of the small-Xg gluon density 
because of the scaling violations. At larger the values of Xg are larger, the scaling 
violations are weaker, and the decrease of the gluon density towards large x„ takes over, 

2 2 

see Fig. 20. Because of such a convex Q'^ dependence of G{xg, Qq) in the log cry — logQ 
plot one must see the convex curve. 



6.1.3 The Q^-dependence: the experimental data 

The change of the character of the Q^-dependence from small to large is best seen in 
the data for light vector mesons. The experimental data on the Q^-dependence of the p° 
production are shown in Fig. 35. The upper plot in Fig. 35 shows the ZEUS data. The 
low-Q^ data, <1 GeV^, are dominated by axiQ"^) and were fitted to the form 

"'"^ ' " W- + mg„)"» ''''' 

where Rp{Q^) = cl/o't and the exponent rip = 2 as dictated by VDM. The dotted 
line shows the result of the fit including the real photoproduction point. The fit yields 
rrieff = 0.66 ± 0.11 GeV which is close to the mass of the p-meson as expected in VDM. 
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Figure 35: The total cross section of diffractive p meson production as a function 
of Q"^. Upper plot shows ZEUS data from [81] fitted together with photoproduction 
point [68]. The curves represent fits using the function cr{Q'^) oc (1 + -R(Q^))/(Q^ + 
mlfjr)'^ (dotted line) and function cr{Q'^) oc 1/ [Q"^ + m?pY (dashed line). The left 
bottom plot shows recent HI data [85] compared with published data [16, 81]. 
The right bottom plot presents the kt-factorization predictions [54, 55] based on the 
oscillator (solid line) and Coulomb (dashed line) wave functions compared with the 
published HI and ZEUS points. 
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Figure 36: The published HI [87] and ZEUS [278, 280] data on the total cross 
section of diffractive meson production as a function of Q^. The kf factorization 
predictions [54, 55] for the oscillator (solid line) and Coulomb (dashed line) wave 
functions are shown for comparison. 
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Figure 37: The preliminary ZEUS data [88] on the (p meson production cross 
section are shown together with the published data shown in Fig. 36. 
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Figure 38: Upper plot shows recent ZEUS results (solid symbols) on exclusive 
J/i/j electroproduction cross section as a function of Q"^ at (W) = 90 GeV [90]. 
ZEUS photoproduction [71] and HI electroproduction (open symbols) [89] cross 
sections are also shown. The full line is a fit to the ZEUS data of the form (176). 
The inner error bars represent the statistical uncertainties, the outer bars are the 
statistical and systematic uncertainties added in quadrature. The bottom plot shows 
the kt- factorization predictions [54,55] based on oscillator (solid lines) and Coulomb 
(dashed lines) wave functions compared with the published data. The solid curve 
is the k±- factorization result upon the normalization to the photoproduction point 
with the rescaling factors shown in the figure, the curves for two wave functions 
merge within the thickness of lines 
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The VDM value of the exponent Up = 2 for the \ow-Q data must be compared to the 
large-Q^ pQCD expectation, n = 4, without scahng violations, see Eq. (81). Indeed, at 
higher Q^, the cross section of p production steadily departs from the dotted line. The 
parameterization of the form 



1 



yielded Up = 2.32 ± 0.10 for ZEUS 95 data [81] at > 5 GeV^, Up = 2.24 ± 0.09 for 
HI 96 data [16] and Up = 2.60 ± 0.04 for preliminary HI data at > g GeV^ [85]. 
The effective exponent Up of the Q^-decrease rises steadily with the region being fitted, 
from 2 in the soft region to 2.6 at the highest Q^, in good agreement with the theoretical 
expectations of the convex cry(Q^). 

Now we comment on the comparison with theoretical predictions from the /c_L-factorization 

[54,55]. As we cautioned in Section 1.4, the onset of the truly pQCD hard regime for the 

light p-mesons requires large Q^, see Eq. (13). At smaller x the applicability of hard pQCD 

description somewhat improves, still the soft-hard decomposition of the gluon density 

shown in Fig. 33 indicates that the nonperturbative soft component of the gluon density 

remains substantial at gluon momentum squared < (1 ^ 2) GeV^. According to the 

2 

discussion in Section 4.8, see especially Fig. 19, in order to have the hard scale Q > (1^2) 
GeV^ one needs > 10 GeV^. The p-production at smaller is strongly affected by 
nonperturbative physics. Still we notice that the fc^-factorization predictions [54, 55] 
reproduce within the overall normalization factor ~ 2 the measured cross section which 
drops by nearly four orders in the magnitude from the real photoproduction to the largest 
value of Q^. One must not be jubilant, though, since the unintegrated gluon density in 
the soft region has been adjusted to reproduce real photoabsorption and DIS at small-Q^. 
At > 10 GeV^ there is a very good agreement between experimentally observed and 
predicted Q^-dependence. In the nonperturbative region the Coulomb WF is doing better 
job at small Q^, but at higher the experimental data deviate from the curve for the 
Coluomb WF and agree better with the results for the oscillator WF. 

The similar pattern is seen in the 0-production shown in Figs. 36, 37. The preliminary 
— 2 

high-Q data shown in Fig. 37 do agree better with the /c_|_-factorization results for the 
oscillator WF. 

Finally, the experimental results for J/ip, are shown in Fig. 38. These data correspond 

2 

to a sufficiently large hard scale Q . Correspondingly, once the theoretical curves are 
normalized to the photoproduction data as shown in the bottom plot of Fig. 38 - here we 
chose a normalization to the ZEUS point, - the curves for the oscillator and Coulomb wave 
functions become indistinguishable and the resulting description of the experimentally 
observed Q^-dependence is very good. 



97 




1 10 10^ 

(Q'+m'^) [GeV'] 

Figure 39: HI and ZEUS measurements [16, 68, 70, 75, 80, 81, 89, 278, 280-283] 
of the total cross sections a{'j*p — > Vp) as a function of [Q"^ + My) for elastic 
p,uj,(f), J/ip and T meson production, at the fixed value W = 75 GeV. The cross 
sections were scaled by SU(5) charge-isospin factors. The curve is a fit by formula 
(177) to the HI and ZEUS p data, and the ratio D of the scaled u, and J/ip 
cross sections to this parameterization is presented in the insert. 
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A parameterization of the dependence of the recent ZEUS data [90] on the J/\l/ cross 
section in the same form as (176) with an appropriate change of the meson mass yielded 

2 

nji^ = 2.44 ± 0.08 [84]. The impact of the mass term in Q on the exponent ny is 
substantial which is well illustrated by the HI analysis [87] of the combined HI and 
ZEUS data on different vector mesons. A fit performed on the HI and ZEUS p data using 
the paramtrization 

with the resuh ai = 10689 ±165 nb, as = 0.42 ±0.09 GeV^ and ag = 2.37 ±0.10 is shown 
in a curve in Fig. 39. Notice, that the so found exponent 03 is close to the recent ZEUS 
result rij/^ = 2.44 ± 0.08. The + My dependence of the HI data on and of the 
combined HI and ZEUS data on J/tjj were found to follow the same parameterization as 
shown in Fig. 39. 



6.1.4 The vector meson production as a probe of the gluon den- 
sity in the proton 

The deceptively simple Eq. (174) suggests that vector meson production cross section 

2 

discriminating among the different models for the gluon density G{xg,QQ). Ryskin was 
the first to make this point [22], for early discussion see [173]. As we saw above, in all 
the cases the color dipole/Zcj-factorization model with the unintegrated glue adjusted to 
the proton structure function data is doing a good job on the Q^-dependence of ayiQ"^)- 
The /c_|_-factorization results for o"y(Q^) for the two extreme parameterizations of the 
wave function are shown in Figs. 35, 36, 38. The choice of the wave function has a 
marginal impact on the predicted Q^- dependence. The two curves do typically envelop 
the experimental data points, but a mismatch of the factor ~ 2 in the overall normalization 
between the theory and experiment can not be eliminated at the moment. 

Here we illustrate the model dependence with the following example. Recent ZEUS data 
on J/\l/ electroproduction were compared in [90] with predictions of the three pQCD 
models by (i) the extended Bloom-Gilman duality [216] based estimates by Martin et al. 
(MKT, [221]) for different NLO DGLAP parameterizations of the gluon density (making 
use of the NLO gluon densities with leading order impact factors is somewhat inconsistent, 
though): MRST02 [284], CTEQ6M [285] and ZEUS-S [286], (ii) leading logQ^ estimates 
with nonrelativistic J/^ wave function by Frankfurt et al. (FKS, [239,240]) for CTEQ4L 
gluon density [142] and (iii) the color dipole model with unitarity corrections by Gotsman 
et al. (GLLMN, [287]). For the heavy J/^' the pQCD scale is sufficiently large 
already for real photoproduction. The results are shown in Fig. 40. After the absolute 
normalization is adjusted to the photoproduction point, the resulting dependence is 
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Figure 40: The recent ZEUS data on Q"^ -dependence of J electroproduction 
at (W) = 90 GeV [90]. The curves in the top box represent the predictions of the 
QCD models MRT FKS and GLLMN (see text), the curves in the bottom plot are 
the MRT results based on different gluon densities in the proton. All theoretical 
curves are rescaled as indicated to fit the ZEUS photoproduction point at {W) = 90 
GeV. The inner error bars represent the statistical uncertainties, the outer bars 
are the statistical and systematic uncertainties added in quadrature. An overall 
normalization uncertainty of ^g^° was not included. 
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not much different from tlie fcx-factorization results shown in Fig. 38 and is consistent 
with the experimental data. The renormalization factors vary from 0.9 for the GLLMN 
model to 2.98 for the MRT model with the MRST02 gluon density. One would conclude 
that the predictions for Q^-dependence from pQCD master formula (174) are to a large 
extent model independent ones. Fortuitous models with good reproduction of the absolute 
value of o"y(Q^) would not be a wonder, but one should not rejoice with that and in 
general must be content if the theoretical and experimental values of cTy(Q^) agree within 
~ 50%. Hereafter we only shall show the /c_L-factorization results [54,55] for the oscillator 
parameterization of the wave functions. 



6.2 The flavour dependence: ratios (Jy/cTf. 



— 2 

The mass term in the scanning radius (8) and the corresponding hard scale Q of Eq. (9) 
change dramatically from the p,uj to (p to J/\E'. For this reason a comparison of 
for the different vector mesons as a function of makes no sense. This is clearly seen 
from Fig. 41 which shows the ratios of cross sections crv{Q'^)/c'p{Q'^) as a function of Q^. 
They exhibit a very steep dependence on from real photoproduction to large Q^: even 
for the the ratio rises by a factor of ~ 3, whereas for the J/"^ it rises by more than 
two orders of magnitude. The experimentally observed Q^- dependence is well reproduced 
within the fc^-factorization approach [55]. 

We reiterate the point of Section 6.1 about the approximate restoration of flavour sym- 
metry if the cross sections are scaled by the non-relativistic factor rnvTiV e"*"e~) and 
compared at equal hard scale Q , i.e., at equal (Q^ + rriy). According to the Review of 
Particle Properties [102] 

^ ^ ' = p : uj : (f) : J/4j = 0.32 : 0.029 : 0.077 :1, (178) 



(the uncertainties in the rhs of (178) from the decay widths vary from several to ~ 7 per 
cent for light to heavy mesons and are not shown). 

The J/ip has been chosen as the reference point as it is the best approximation to a 
nonrelativistic quarkonium, for light vector mesons the non-relativistic approximation 
is evidently poor. As a test of theory one must rather compare the predictions for the 
flavour dependence at a sufficiently large {Q^ + rriy) ~ 20 GeV^. Here the /c_L-factorization 
predictions can be summarized as 

1 



Vv 



= 0.68 : 0.068 : 0.155 : 1. (179) 

fact 



There is a marginal change from the oscillator to Coulomb parameterization for wave 
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Figure 41: The ratios of the uj [86], (f) [87] (the PHP point is calculated using 
[68, 278]) and J /ip (calculated using only recent ZEUS data [68, 71, 81, 90], for HI 
results see Figs. 42, 77) to the p° cross sections as a function ofQ^. The horizontal 
solid lines correspond to the SU(4) predictions, while the horizontal dashed lines 
correspond to the pQCD predictions in the non-relativiatic approximation given by 
Eq.(178). As shadowed band the corresponding predictions from the kt-model for 
the oscillator wave functions [55] are shown. The width of the band is taken just 
arbitrary and does not correspond to the theoretical uncertainties. The kt-model 
prediction for the J/ip/ p ratio at Q^=100 GeV^ is shown separately as a dashed- 
dotted line. 
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Figure 42: A compilation of the flavor-res caled using k±-factorization Eq.( 179) 

total cross sections rjy^aviQ'^) f^f the p, u, 0, J/\l/ vector mesons [16, 68, 11, 81, 86, 
87, 89, 90, 278, 280, 288] as a function of -\-my. Shown are also the corresponding 
predictions from the kt-model for the oscillator wave functions. [55]. 
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functions, the theoretical uncertainty in the rhs of (179) is about 10 per cent. In Fig. 42 we 
plot the fiavor-rescaled (Tv{Q'^)'nv^ \k^-ia.c\. vector mesons as a function of (Q^+my). 

The universal (Q^ + my)-dependence is clearly seen and the overall agreement between 
the experimentally observed and theoreticaly predicted flavor dependence is good. 

For uncertain reasons, perhaps by misinterpretation of the discussion in [239], one often 
discusses the restoration of the SU{A) flavour symmetry, 

p° : : : J/?/^ = 1 : 1/9 : 2/9 : 8/9 = 1.125 : 0.125 : 0.22(2) : 1. (180) 

Incidentally it is not much different from (179), except of J/ip., see Fig. 41. We caution 
that there are no sound reasons for such an SU(4) ratio of cross sections even at very 
large Q^. The S'?7(4)-factor weighted empirical cross sections were shown in Fig. 39 and 
can indeed be fitted by the flavor-independent universal curve, but this agreement with 
SU (4) rations must be regarded as an accidental one. 

6.3 The ratio Ry = ctl/ctt 
6.3.1 Theoretical expectations 

For heavy vector mesons treated as a nonrelativistic quarkonium, the pQCD predicts 
Rv{Q'^) = ^Rlt with Rlt = 1, see (148). Despite all the uncertainties with the wave 
functions, important point is that for a meaningful evaluation of Rv{Q^) the transverse 
and longitudinal vector meson must be related by the rotation-invariance. The Fermi mo- 
tion effects in the longitudinal and transverse amplitudes are different, though, the latter 
being more sensitive to the end point contributions from 2 ~ 0, 2; ~ 1. Such corrections 
are automatically incorporated in the color dipole and /c^-factorization calculations. First, 

as discussed in Section 4.8, the hard pQCD scales for transverse and longitudinal cross 

2 2 

sections do slightly differ, Qqj^ < Qgl^ which already leads to a substantial reduction of 
Rlt at large Q^, as it was found in [21]. Second, the SCHNC transitions contribute for 
the most part to the transverse cross section, especially for light vector mesons, further 
lowering Rlt- Third, as discussed in Section 4.10, the predictions for Ry^Q"^) are poten- 
tially very sensitive to the presence of a hard, quasi-pointlike qq component in the vector 
meson wave function. As a test of this sensitivity we evaluate RviQ"^) for a squeezed 
p-meson as discussed in Section 4.10. 

One way to circumvent problems with the poorly known wave function of vector mesons is 
to resort to the Bloom-Gilman duality for diffractive DIS [216]. Because of the strikingly 
different dependence of the duality integrals for 0"^, ctt on the duality interval, see Eqs. 
(95), (95), the ability [220,221] to reproduce both cty{Q'^) and Ry^Q^) is not surprising. 
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It is not clear, though, whether such an evaluation is entirely consistent with the rotation 
invariance constraints or not. 

In the color evaporation model discussed in Section 3.4.7 the J/\E' is a part of the open 
charm cross section and, arguably, the ratio 0^1 (Jt must be close to that for the open 
charm, Ry ^ Rdis ^ which would conflict the large-Q^ data on J/ip production to be 
shown below. 



6.3.2 Experimental results 




Figure 43: A compilation of the experimental data on the ratio Rp = ctl/ctj' 
and R^j: = Rpm'^p/Q'^ for the p meson production as a function of Q"^. The recent 
preliminary data from HI [85] and ZEUS [83] measurements are just scanned 
from the plots of the correspondent papers and shown together with the published 
results [16, 68, 81, 277]. The theoretical predictions are as follows: the dotted curve 
shows the estimates by Martin et al. based on the Bloom- Oilman duality [220], 
the solid curve and dashed curves show the results from the k±- factorization for the 
conservative radius of the p-meson and the squeezed p-meson, respectively [54, 55]. 



Figure 43 shows the ratio 

-Kplv j - — r~* V ' 

o-t(7 P PP) 

as a function of for the p meson production. The right plot shows the same data 



105 



presented in the form of 

2 

TTl 

At the small values of the ratio is an approximately linear function of and Rlt{0) ~ 
1, but Rlt{Q^) decreases steadily with rising and the growth of Rp{Q'^) slows down. 
Whereas some of the published data sets even suggest flattening of Rpi^Q"^), the new 
preliminary data [83,85] indicate the steady large-Q^ rise of RpiCf') 
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Figure 44: A compilation of the experimental data on the ratio Rj/^ = Ui/aT 
and RjP^ = Rj/^m'^j^^/Q'^ for the J/i/j meson production as a function of Q^. The 
open circles represent the HI [89] and solid symbols the ZEUS J/ip measurements 
[90]. The solid curve is a prediction from the k±- factorization approach [55]. 

The dotted curve shows an evaluation Rp{Q'^) [220, 221] based on the duality approach 
[216]. Whereas in the k±_ factorization one starts with the qq continuum production 
amplitudes, projects it onto the = 1~ state and averages over the masses M ~ my 
of the qq state with the vector meson wave function as the weight function. In the 
duality approach [216] one rather calculates |Tp for a fixed diffractive mass M and then 
integrates over certain range of M around my. For the p meson Martin et al. take the 
duality interval M = [Mi, Ms] = [600, 1050] MeV. Because of the different sensitivity of 
duality integrals for ctl^t to the duality interval, see Eqs. (93), (95) one readily finds a 
good description of the dependence of Rp{Q'^) and the expense of a possible conflict 
with the rotation invariance. Although the absolute value of Rp{Q'^) is essentially adjusted 
to the data, the steady rise of Rp{Q'^) at large must be regarded as the genuine pQCD 
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prediction. For the uncertainties with the duahty description of the absolute cross section 
see Section 3.4.7. 

The sohd curve shows the result of the /c_|_-factorization calculation with the conservative 
radius of the p-meson. Although at largest ~ 20GeV^ the predicted Rlt ~ 0.3 <^ 1, as 
found already in 1994 within the color dipole model [21], there is an obvious disagreement 
with the large-Q^ data. A comparison between the data and the /cf-factorization calcu- 
lations in terms of and cr-r separately shows [54] that the source of this discrepancy 
is the transverse cross section ar, whose value at high is significantly underestimated 
by the model. A strong sensitivity of predictions for Rp{Q^) to the short-distance prop- 
erties of the vector meson is illustrated by the dotted curve which is the fc^-factorization 
result for the squeezed p- meson. As anticipated in Section 4.10, the prime effect of higher 
short-distance density in the p-meson is an enhancement of ax and suppression of Rp{Q^). 

Figure 44 shows the ratio OhjaT for the J j^l) mesons compared with the prediction of the 
/cj-factorization model. The recent ZEUS result averaged over is Ri^t = 0.52±0.16 [90]. 
The recent high accuracy experimental results for large-Q^ production are shown in 
Fig. 45. If compared against equal Q^/m^, the theoretical results and the experimental 
data for Ri^^ show similar behavior for all three vector mesons: p, </> and J/^/'. For 
instance, the fit to experimental data on Tqq for the p-production shown in Fig. 25 can be 
reinterpreted, via Eq. (164) with e «i 1, as -Rp = a((5^/m^)^ with a = 1/^ = 0.46 ± 0.015 
and b = k = 0.75 ± 0.3 which agrees perfectly with the recent ZEUS parameterization for 
the (j) production shown in Fig. 45: a = 0.51 ± 0.07 and b = 0.86 ± 0.11. 
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Figure 45: A compilation of the experimental data on the ratio R(p = ci/aT 
for the (j) meson production as a function of . Shown are the recent preliminary 
data from ZEUS [88] together with the published results from ZEUS [278, 280] 
and HI [87]. The solid curve shows the fit of the form R^ = a{Q'^ /m?^^ with the 
parameters cited in the figure. 
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7 The energy dependence and Regge properties of 
difFractive vector meson production 



7.1 Theoretical expectations 

The local energy dependence of o-'y*p~,vpiW^,Q'^) is usually parameterized as 

a^*p^Vp{W^, Q^) oc . (181) 

The exponent 6 for the t-integrated total cross section is controlled by the x-dependence 
of the integrated gluon density, see Eq. (145), and the shrinkage of the diffraction cone. 
Within the standard exponential approximation for the t-dependence 

and 

^ = - ■ (183) 

KQg) 

The microscopic pQCD models of the Pomeron exchange, see Sects. 3.2, 3.3 and Eq.(60), 
suggest that the Pomeron is a more complex object than an isolated single pole - it is 
either the branching point or a sequence of Regge poles. For this reason the exponent 6 
can depend on Q^, my, the energy range, helicities etc. Eq. (183) emphasizes the crucial 
point [21] that within the color dipole and /c^-factorization approaches 6 only depends 
on Qq ~ (Q^ + my). The corrections to the (Q^ + my)-scaling are well understood and 

2 2 

marginal: a slight departure of the hard scale Qq from Q has been discussed in Sect. 
4.8. Specifically, at equal values of + my, the p production is in a somewhat softer 
regime than the J/ip production and one can expect that Sj/^ should be slightly larger 
than for the 6 p. 

The /c^-factorization phenomenology of DIS structure functions described in Sect. 3.3.3 
strongly suggests that, for purely numerical reasons, within the kinematical range of the 
HERA experiments the vacuum exchange can be approximated by the two, soft and 
hard, Pomeron poles with approximately Q^-independent intercepts. The values of the 
intercepts anticipated in this approximation were given in Sect. 3.1.4 and Sect. 3.3. They 
translate into 6 ~ 0.2 ^ 0.4 in the soft region and much larger 5 > 1 ^ 1.5 when the 

2 

interaction becomes sufficiently hard. Since the pQCD hard scale Q increases not only 
with Q^, but also with my, the theory predicts the corresponding hierarchy 6j/^ > 6^ > 6p 
at any Q^, including photoproduction; for exact values see Fig. 46. On the other hand, 
the variations of 6 from p,uj to (j) to J/\l/ to T are weak against the variable + my. 
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Figure 46: The kf -factorization predictions for the exponent 5 of the fits a oc 
to the W dependence for different vector mesons [54, 55]. Notice a strong flavor 
dependence in the left box where 5 is plotted as a function of . The flavor 
dependence is much weaker when 5 is plotted against + My (the right box). 

The energy dependence of the ^ji^S) production should be discussed separately. As 
pointed out in Section 3.4.3, the node of the 25* radial wave function leads to a partial 
compensation between contributions from dipoles with the size below and above the node. 
As a result, the typical dipole sizes that contribute to ip{2S) production are smaller 
than for J/ip production. This makes the pQCD scale for ip{2S) production somewhat 
harder than for J/ip, which results in 6^^2S) > ^j/tp and the rise of the cross section ratio 
a{ip{2S)) / a{J / ip) with growing energy is predicted [26,54]. 

7.2 Experimental results: real photoproduction 
7.2.1 Ground state vector mesons 

The summary of the experimental data on the energy dependence of all vector meson 
photoproduction cross sections measured at HERA is shown in Fig. 47 together with the 
total cross section cr^jj and the results from fixed target experiments. There is a clear 
pattern of variation of the exponent 5 from light to heavy mesons: For the light vector 
mesons, cu, p and 0, the rise of the cross sections with energy is well described by W^'"^^ 
behavior. Already the first measurement of the J/ip photoproduction at HERA showed 
a much steeper rise with W than that observed for the light vector mesons. The energy 
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Figure 47: The total vector-meson photoproduction cross section as a function 
of W for all vector mesons measured at HERA shown together with the results of 
fixed target experiments and compared with the total photoah sorption cross section 



7P 
tot- 
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dependence was found to be W^^''^~^'^\ which imphes the effective Pomeron intercept 
around aip{0) ~ 1.2. The still larger 6 has been found for \E'(25'), the experimental 
data on T production are not conclusive yet. Figure 48 summarises of the experimental 
data on the p photoproduction and Fig. 49 shows the recent ZEUS photoproduction 
measurement [71] of the J/tp photoproduction in a wide W range (20 -h 300 GeV) together 
with earlier data from HI and fixed target experiments. The result of the fit in the form 
a ocWHo the ZEUS data with > 30 GeV yielded 6 = 0.69 ± 0.02(stat.) ± Om{syst.). 




Figure 48: The W-dependence of exclusive p photoproduction cross section mea- 
sured by the HI [277] and ZEUS [68] experiments compared with predictions of 
the kf-factorization model for the oscillator and Coulomb wave functions of the p 
meson. 



The fcf-factorization and color dipole calculations, shown in Fig. 46, predict precisely 
such a variation of the energy dependence with the mass of the vector meson. The real 
photoproduction of light vector mesons is dominated by non-perturbative component of 
the unintegrated gluon density of the proton and/or soft dipoles, and 6 is small, but the 
contribution from hard gluons rises from light to heavy vector mesons, which is a universal 
prediction from all pQCD motivated models. The energy dependence from the specific 
/c_L-factorization model [54,55] comes out right for both p (Fig. 48) and J/\l/ mesons 
(Fig. 49). 
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Figure 49: A compilation of the experimental data on exclusive J/ip photopro- 
duction cross section as a Junction of W measured. The black points show the 
recent ZEUS data for the Jj'^ —>■ and J/ip e+e~ decay channels [71]. 

The inner bars indicate the statistical uncertainties; the outer bars are the statis- 
tical and systematic uncertainties added in quadrature. Results from the HI [70] 
and fixed target experiments [289, 290] are shown by open symbols. Recently HI 
reported about correction of the published cross sections and they are now much 
closer to the ZEUS ones. The kf factorization predictions [55] for the oscillator 
and Coulomb wave functions are shown by the solid and dashed lines, respectively. 
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7.2.2 Radially excited vector mesons 



The \1/(2S') photoproduction was studied by HI [73] using tagged and untagged data 
samples of \1/(2S') II and \I/(2S') lliriT decay channels, where // stands for e~^e~ or 
fi'^fi". The measurement of the ratio R of "^{23) to J/ip photoproduction cross sections 
as a function of W is shown in Fig. 50. The suppression of this ratio, i? -C 1, is due to 
the node effect, the resulting strong cancellations between the contribution from the large 
size and small size components of the "^{28) make the predicted ratio strongly dependent 
on the model for the wave function. The experimental data agree with /cj-factorization 
predictions based on the oscillator wave function. The overall ratio for the data gives the 
value R = 0.166 ± 0.007(stat.) ± 0.008{sys.) ± 0.007{BR) in a good agreement with the 
previous measurements [72] (for the branching error calculation see [73]). A fit to this 
ratio of the form R oc W^^ yields A6 = 0.24 ± 0.17, which indicates that the energy 
dependence of the ^^(25) photoproduction cross section is slightly steeper than that for 
the J/ip meson. This agrees with the color dipole model expectations [26,49]. 




Figure 50: The ratio of photoproduction cross sections \&(2S') over J/ip as a 
function ofW [73]. The kt-factorization predictions [55] are shown for comparison. 

The definitive p'{2S) assignment of the excited p states is still pending and the experimen- 
tal data on the p'{2S) production are not available yet. Here we simply mention that all 
color dipole [178,181,182] and/or /c^-factorization [54,55] calculations based on the node 
effect predict a{p' {2S)) / a{p{lS)) < 1 with very steep dependence up to > (3 -i- 5) 



113 



GeV^, whereas the duahty relation estimates by Martin, Ryskin and Teubner [220] gives 
(j{p' {2S)) / (j{p{lS)) > 1. This simply shows that one must be careful with application of 
duality to radial excitations: because the mass spectrum of the qq pairs which enters the 
duality integral does not exhibit any non- monotonous Q^-dependence [3,217], the steep 
Q^-dependence driven by the node effect can easily be missed. 



7.2.3 Test of the vector dominance model 

In the color dipole language, the success of the vector-dominance model (VDM) for real 
photoproduction derives from the proximity of the (quark flavor dependent) color dipole 
distributions in the photon and light vector mesons. The straightforward extension of the 
VDM approximation (67) to the jp elastic scattering amplitude reads 



T(7P-7P)= ' y^ Ti^p^Vp). (184) 

Assuming the pure imaginary amplitudes, which is a good approximation because in 
real photoproduction ^5 <^ 1, one can extract the photoproduction amplitudes from the 
forward differential cross sections, whereas the jp elastic scattering amplitude is related 
by optical theorem to crtotiip)- This leads to the Stodolsky sum rule [291] 



(^tot{lp) = Wl67r ■ 



dt 




lOTT ■ 



1=0 ,/_„o,„^ V JV 



3 dt 



(185) 

i=0 



The test of VDM sum rule has been reported by the ZEUS collaboration [292]. The 
VDM analysis of the low-energy data gave /y/47r = 2.20, 23.6 and 18.4 [42] for p°, u 
and 0, respectively. Then, based on the photoproduction data at 70 GeV, the VDM 
sum rule gives a value of 111 ± 13(exp.) ph for the photon-proton total cross section at 
W^p = 70 GeV. The meson contributes about 85% of this value. The interpolation of 
photon-proton total cross section at a center-of-mass energy of W-yp = 70 GeV, obtained 
by interpolation between the present measurement and the lower energy measurements 
using the Regge model fits, is 139 ± 4 /xb. The two numbers are close to each other; the 
point that simplified VDM model does not saturate the sum rule is well known, for the 
review see [42]. 



7.3 Experimental results: vector mesons in DIS 

7.3.1 The impact of hard scale on the energy dependence 

The cross section for the exclusive p° electroproduction measured by ZEUS [83] and HI [85] 
as a function of W for different values of is presented in Fig. 51. The curves represent 
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Figure 51: The W dependence of the p electroproduction cross section for different 
values of as measured by ZEUS [83] and HI [85]. The lines represent the results 
of fitting a oc with the 6 values indicated in the figures. The shaded area 
in the ZEUS case indicates additional normalization uncertainties due to proton 
dissociation background. 



the result of the a oc fits to the data. The exponent of the energy dependence increases 
with the growth from about 0.2 at low up to 0.8^ 1.0 at high Q^, in agreement with 
theoretical expectations form the color dipole approach shown in Fig. 46. In contrast to the 
p production, the J/ip production cross-section exhibits almost the same ly-dependence 
for all measured values of Q^, including the photoproduction limit, see Fig. 52, in good 
agreement with the theoretical prediction from the color dipole approach shown in Fig. 46. 

The point that the correct hard scale is oc {Q^ + my) is clear from Fig. 53 which shows 
the dependence of the exponent 6 on (the Ihs. box) and (Q^ + my) (the rhs box), 
respectively. For the p-mesons the both plots give a clear indication of the rise of S with 
and Q . At largest the exponent 6 rises up to 1 -^ 1.2, which implies the rise of 
the effective intercept of the Pomeron up to aip(O) ~ 1.3. If it were to be plotted as a 
function of Q^, the exponent S for the J/\l/ would be completely out of the observed trend 
for the p-meson, see the Ihs box where only the PHP point is shown. When plotted as a 
function of {Q^ + my), the same result for the J/\l/ is perfectly consistent with the results 
for the p at the same value of (Q^ + fny), confirming with the theoretical expectation of 
the (approximate) fiavor symmetry in this variable. The theoretical values of 6 shown in 
Fig. 53 were evaluated for the range W = 50 ^ 100 GeV. The overall agreement between 
the experiment and /c^-factorization approach is good. The lower box of Fig. 52 shows 
shows the /ct-factorization predictions [54, 55] normalized to the ZEUS photoproduction 
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Figure 52: T/ie to'p plot shows the recent ZEUS results on exclusive J/ip elec- 
troproduction cross section as a function of W for four values of [90]. ZEUS 
photoproduction [71] and HI electroproduction [89] cross sections are also shown. 
The full lines are fits to the ZEUS data of the form a oc . The inner error bars 
represent the statistical uncertainties, the outer bars are the statistical and system- 
atic uncertainties added in quadrature. An overall normalization uncertainty of^^^ 
was not included. The bottom plot shows the kt- factorization predictions [54, 55] 
normalized to the ZEUS photoproduction cross section at {W) = 90 GeV and com- 
pared to the recent ZEUS data. 
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Figure 53: The exponent 5 of the fits a oc to the W dependence for the p 
production as summarized in [83] including the PHP J/ijj point from [70] presented 
as a function of Q"^ (left plot) and combined p and J/ip data [71, 89, 90] presented 
as a function ofQ'^ + My in the right plot. The predictions from the kt-factorization 
approach [54, 55] are also shown. 

cross section at (W) = 90 GeV; tlie quality of the theoretical description of the data 
points is good. 



7.3.2 Discriminating the models for gluon density 

The early discussion of the J/ \E' production was centered around the idea [22] of stringent 
tests of models for the gluon density Q^)- An example of such a test based on the 
experimental data which were available at that time is illustrated by the Ihs plot of Fig. 54. 
The shown theoretical curves [173] are based on the 1995 updates of parameterizations for 
the gluon density adjusted to the HERA data on the small-x growth of F2p{x,Q'^). The 
energy behavior of each vector meson production can also be reproduced by a number of 
simple parameterizations within models described in Sects. 3.1.4 and 3.4.8. An example 
of a fit to the J/\l/ photoproduction cross section in a model by Fiore et al. [229] in which 
Pomeron is the double-pole at ajp(O) = 1 is shown in the rhs plot of Fig. 54. 

The past decade brought great improvements in our understanding of the small-x gluon 
densities and in the vector meson productuon data compared to the 1995 situation shown 
in Fig. 54. The discrimination of the models remains weak, though: for the reasons 
discussed in Section 6.4, there are substantial uncertainties in the absolute normalization 
which can not be eliminated at the moment. This point is further illustrated by the 



117 




Figure 54: The left box: tests of models for the small-x gluon density using the 
early HERA data on the energy dependence of J/'^ photoproduction together with 
data from fixed target experiments. The plot taken from [173] shows the pQCD 
calculations for the GRV95 [293] and MRS95 [294] parameterizations of the proton 
gluon density based on the pre-1995 data on F2p{x,Q'^) and illustrates the status 
of the subject at that time. The right box: the parameterization of the energy 
dependence of J ftp production in the Fiore et al. model of a double-pole Pomeron 
with intercept ajpifS) = 1 [229]. 
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Figure 55: The tests of the pQCD predictions for various models for the small- 
X gluon density of the proton vs. the recent ZEUS data on J/ip production as a 
function ofW for different values ofQ^ [71, 90]. The curves in the top box represent 
the predictions of models MRTxl.49 (solid curve) FKSxl.7 (dashed curve) and 
GLLMNxO.9 (dotted curve), the curves in the bottom plot are the MRT results 
based on different gluon densities in the proton: the solid curve - ZEUS-Sxl.49, the 
dashed curve - CTEQ6Mx2.22, the dotted curve - MRST02x2.98. All theoretical 
curves are rescaled as indicated above to fit the ZEUS photoproduction point at 
(W) = 90 GeV. The inserts show the exponent S of the parameterization a oc 
for these models. 
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recent ZEUS data on J/\[^ electroproduction shown in Fig. 40 in comparisons with the 
predictions from the more recent MRT, FKS and GLLMN models for the gluon density 
described in Sect. 6.4. All models need large rescaling to adjust the normalization to the 
photoproduction at (W) = 90 GeV. At the expense of such rescaling all model calculations 
reasonably reproduce the gross features of the energy dependence. For comparison, the 
fcj-factorization predictions [54, 55] shown in the bottom plot of Fig. 52 give an equally 
good description of the same data using a rescaling factor 0.86. One exception is the MRT 
results for the MRST02 gluon density for small Q^; the likely reason for the found strong 
discrepancy is that the corresponding values of Q are too close to the lower Q^-boundary 
of applicability of this particular set of parton densities [284]. 



7.3.3 Comparison of vector meson production and inclusive DIS 
from Regge model viewpoint 

We repeatedly made a point that the pQCD vacuum exchange is not an isolated Regge 
pole. Section 3.3 and Fig. 11 show how the effective exponent of the ^-dependence changes 
from r(K^) for the unintegrated glue J-'{x, h?) to A(Q^) for the integrated glue Q"^) 
to A(Q^) for ^2(0;, Q'^) or aU{x, Q"^). This entails the failure of naive Regge factorization 
in a comparison of the energy dependence of vector meson production and inclusive DIS, 
which was nicely demonstrated by ZEUS collaboration [295] in their study of the ratio 

a,.^y,{W\Q^) ^ (186) 
The argument is as follows. According to Eqs. (181), (183), 

(187) 

where the last form in terms of a^p holds if the QCD Pomeron were an isolated Regge- 
pole. Notice, that A(Q^) as defined in Sect. 3.3.1 is different from the exponent A((5^) 
in 

see a comparison in Fig. 10. Such a difference between the intercepts A(Q^) and A((5^) 
and their substantial dependence on shown in Fig. 10 already go beyond the rigors of 
the Regge theory. Still one can try to probe the vacuum exchange in the numerator and 
denominator of ry equalizing the relevant hard scales, i.e., evaluating the ratio r]^^ with 
^^7o?(W^^Q^) taken at = (Q^ + m^)/4. In Fig. 56 we present AL(Q^) and A(Q^) 
plotted vs. common hard scale. Only if one ignored the difference between \{Q^) and 
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A((5^), the dependence in (187) and (188) would be controlled by one and the same 
ajp and one would expect the TV-dependence of the ratio r]^^ of the form 

rlt^{W')\ ^ = «^-7;^-l- (189) 
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Figure 56: A compilation of the exponents X^g^{Q'^) = ^6{Q'^) and A{Q ) of Regge 
fits to the energy dependence of the p photoproduction cross section, CTpiQ"^) oc 
(14/2^2Af„,(Q2)^ anc? of the proton structure function [296, 297], F2p{x, = i(Q^ + 

m?p)) oc a;~^^*^ \ plotted against {Q"^ + m^), in comparison with the results from the 
k ^-factorization model [54, 55]. 



The experimental results on rYot for the p and J/\l/ are shown in the Ihs plots of Fig. 57. 
They are consistent with little or no ly-dependence for the p production. Notice, that an 
approximate constancy of ry for the p-production at small is very much reminiscent of 
the familiar very weak energy dependence of the ratio (Jei/<ytot in ttA^, KN., NN interaction, 
see plots in [102]. Here the smallness of the exponent r/ is to a large extent due to the 
term a'jp/hiW'^) from the shrinkage of the diffraction cone. The ly-dependence of r]^^ for 
the J/\l/ production is substantial. This hard-scale and process-dependence of rj has been 
interpreted as an evidence for the breaking of the Regge factorization. 

The energy dependence of r]^j expected from /cj-factorization is shown in the rhs boxes 
of Fig. 57. It includes the effect of shrinkage of the diffraction cone. The theoretical 
results do correctly reproduces the trend of the experimental data shown in Ihs boxes of 
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Figure 57: The ratio ry = <Jv{Q'^)/o'tot{Q") for the p production as function of 
energy for several values ofQ = \{Q'^ + my): the Ihs plots show (a) the experimen- 
tal data on the p and J ftp production [295], the rhs plots show the corresponding 
theoretical expectations for the energy dependence of ry [55]. 



122 



Fig. 57. Although the ZEUS experimental results for the p mesons are consistent with 
^Yot — const, within the error bars they do not exclude the theoretically expected weak 
energy dependence shown in the rhs plot. The /c^-factorization correctly describes the 
change of the energy dependence of r^ot from the light p to heavy J/\E'. 
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8 The t-dependence and properties of diffractive cone 



8.1 Low-i: proton-elastic and proton-dissociative events 
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Figure 58: Schematic diagram of proton- elastic (left) and proton- dissociative 
(right) vector-meson production in ep interaction: ep — > eVp and ep — > eVY . 

The vector meson production differential cross section da/dt exliibits a pronounced for- 
ward diffraction cone, wliicli spans up to \t\ < 1.0 GeV^. Sucli a diffraction cone is familiar 
from liadronic scattering processes. Witliin diffraction cone, tlie differential cross section 
falls off with |t| approximately exponentially, see Sect, 3.1.3 and Sect. 3.1.4. In close sim- 
ilarity to hadronic scattering, the dominant process here is the proton-elastic (hereafter 
just "elastic") production of ground state vector mesons 7p —>■ Vp, see Fig. 58 left, where 
the label "elastic" is a reminder of the VDM relationship between the photoproduction 
and VP elastic scattering amplitude, recall Section 3.4.1. 

At larger t the elastic production dies out and the proton-dissociative reaction 7p —>■ 
VY takes over. One can argue [114] that the relative importance of the elastic and 
proton-dissociative reaction is precisely the same as in proton-nucleus and proton-proton 
scattering described in Section 3.1.5. At small t within the diffraction cone the proton- 
dissociative production will be smaller than the elastic production but still the biggest 
background contribution to the elastic vector meson production. 

Experimentally, the direct separation of the two processes is only possible if the leading 
proton is measured in the Leading (Forward) Proton Spectrometer (LPS, FPS) or if the 
hadrons from the proton-dissociative system Y with a sufficiently high mass (> 3 ^ 4 
GeV) are observed in the forward part of the detector. Because the forward part of each 
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detector has a beam-pipe hole the smaller mass states Y will just escape undetected. In 
the case when the proton or its excitation escapes through the beam-pipe hole undetected, 
one needs to estimate the portion of the proton-dissociative events and subtract it from 
the visible cross section based on the Monte Carlo modeling. Such a procedure leads to 
sizeable systematic uncertainties. 

If the Regge factorization decomposition (35) were exact, then the t-dependence from the 
7 — > transition vertex and from the t-channel exchange will cancel in the ratio 



Rat^o{el/d^ss) = ^ /da{^p-.VN) ^^^^^ 

dt / dt 



which must be flavor independent and controlled by the change of the t-dependence from 
p — > p to p — > y transition [114]. Furthermore, the Regge factorization predicts that 
Ratio{el / diss) will not change from photoproduction of vector mesons to pp scattering, 



Rat^o{el/d^ss) = ^-iPP -^PP) / MpP pN) ^ ^,,,,(,^/^,,,) . 



In Fig. 59 the experimentally measured ratio Ratio{el / diss) is shown as a function of t 
for p, and J/ ip mesons photoproduced at HERA. The same ratio for the pp reaction is 
shown for a comparison. All ratios coincide within the errors, which is consistent with 
the hypothesis of Regge factorization and lends a support to the separation procedure. 
The elastic reaction production dominates at |t| < 0.4 — 0.6 GeV^, the proton-dissociative 
production takes over at larger Notice a close similarity to a comparison of the elas- 
tic and nucleus-dissociative pA scattering shown in Fig. 8. Figure 59 shows the typical 
interplay of the t dependence of the elastic and proton-dissociative production in electro- 
production of J/ijj as measured by HI [89]. One sees that the relative contribution of the 
proton-dissociative events is increasing rapidly with t making the elastic measurement 
possible only at rather low t values (typically |t| < 0.6 ^ 1.0 GeV^). 

At still larger momentum transfers, |t| ^ 1 GeV^, the vector meson production is dom- 
inated by proton-dissociative reaction. Besides that, the production mechanism changes 
and the exponential fall-off of the differential cross section is superseded by the inverse 
powers of \t\ or {my + \t\). 



125 




0.5 1 1.5 2 

Itl [GeV'] 



Figure 59: Upper plot shows the ratio of elastic to proton-dissociative differential 
cross sections as a function of —t for vector-meson photoproduction [76], together 
with data from pp reactions [298] at ^/s = 23.4 and 38.3 GeV. The curve is result 
of combined fit to all the data. Bottom plot presents the \t\ distribution for the 
J/i[ sample as measured by HI [89]. The dashed line is result of a fit taking 
the background contributions into account. The full line corresponds to the elastic 
contribution assuming an exponential distribution. The contributions from proton 
dissociation and non-resonant background are shown separately. 
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8.2 Diffraction cone at low-i 



8.2.1 The diffraction slope versus hard scale: theoretical expec- 
tations 

The t-dependence of the differential cross section as small-t is usually parameterized in 
terms of the diffraction slope by, see Eq. (34). The fitted values of by depend slightly on 
the range of t, the more refined parameterization 

d(T / 7 I I I i2\ 

— OC exp(— Oy|t| — c\t\ ) 

(JjO 

with the curvature parameter c allows one to extend the fits of the experimental data 
up to |t| ~ 1 GeV^. The allowance for the curvature c does not shift the value of fitted 
by significantly. Several other definitions of the effective slope can be encountered in 
literature, e.g., the derivative of the logarithm of the differential cross section at t = 0, 



, 1 da 
bv = -- — 

a at 



t=o 



or, less often, the definition in terms of the average momentum transfer squared, 



they can differ from by defined by fits to Eq. (34) by 1 ^ 2 GeV^^. 

The decomposition (35) of the diffraction slope into the target transition, beam transi- 
tion and the t-channel exchange components is exact in the simple Regge model. Simi- 
lar decomposition holds also in the /c^-factorization and color dipole approaches despite 
breaking of the strict Regge factorization. In the photoproduction the key point is the 
shrinkage of the qq state of the photon from light to heavy quarks accompanied by the 
related decrease of the radius of the ground state vector meson, see Section 3.4.2, which 
entail the hierarchy of diffraction slopes 

bj)^ < bt < bf . (192) 

The Q^-dependence of the diffraction slope is driven by the decrease of the scanning radius 
rs with Q'^, see Section 3.4.2, Eq. (80), and the Q^-dependence of the Regge shrinkage 
term through the Regge parameter W^/IQ"^ + i^v), see Section 4.5, Eq.(128): 

W^)=.„ + 2a'rlog(^^)+^^ (193) 
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For xq = 8.3 ■ 10"^ as defined in Section 4.5, ^ Bj^i is the approximately flavor- 
independent contribution from tlie proton. Tlie slirinkage of tfie diffraction cone is for the 
most part a property of the x^-dependent skewed gluon density, see Section 4.5. Besides 
that, because the scahng violations do depend on x„, the scanning radius and hard scale 

2 

Q will exhibit slight energy dependence, which shall affect the energy dependence of 
the differential cross section [26,299]. The net effect of such non-Regge corrections is a 
substantial reduction of the observed a'p from the input a'sFKL — 0-25 GeV~^ in the 
parameterization (129). 

According to Eq. (80), the diffraction slopes for different vector mesons should be close 
if taken at equal values of the scanning radius rs or equal values of (Q^ + my) [25]. 
More detailed analysis in [26] has shown that for the J/ijj the diffraction slope is slightly 
smaller, by ~ 0.5 GeV~^, than for the p at the same {Q'^+rriy). In addition, the arguments 
presented in Section 4.8 suggest a somewhat larger scanning radius and larger diffraction 
slope for the transverse amplitude Tn compared with the longitudinal amplitude %o. The 
effect has been suggested and evaluated for the flrst time in [26]. Even for light vector 
mesons at small the expected difference is small, ~ 0.5 ^ 1 GeV~^, for heavy vector 
mesons the contribution to the diffraction slope from the scanning radius is small and the 
effect is negligible one. 

The arguments of Section 3.1.5 are fully applicable to vector meson production. For the 
proton dissociative photoproduction one expects a substantially smaller diffraction slope 



The principal point is that the difference of diffraction slopes for the elastic and proton- 
dissociative reaction must be approximately and flavor independent. 

8.2.2 Experimental results: real photoproduction 

The experimental results on the diffraction slope b measured in photoproduction at HERA 
are summarized in Tab. 4. 

The values of diffraction slopes depend on the t-region, where the flts are performed, which 
is familiar from 7Tp,pp scattering [111,300]. The flavor dependence is consistent with the 
theoretical expectations described in Section 8.2.1: the diffraction slope rises with the size 
of the meson with the exception of the \1/(2S'). The charmonium models [47,48] suggest 
the radius of the ^^(25) as large as the radius of the 0, so that the diffraction slope of 
elastic ^^(25")^ scattering would be the same as in (pp elastic scattering. Then in the naive 
VDM one would expect 6^(25) ~ In contrast to that, the counterintuitive 
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Process 


Value of slope b, GeV ^ 


W, GeV 


\t\,GeV^ 


Exp. 


Ref. 


7P pp 


10.31 ± 0.77{stat.) ± 0.52{sys.) 
10.9±0.3{stat.)tliisys.) 
6.0 ± 0.3{stat.)toi{sys.) ± OA{mod.) 


25-70 
50-100 
85-105 


0.073-0.45 
0.-0.5 
0.4-1.2 


HI 
ZEUS 
ZEUS 


[69] 
[68] 
[76] 


IP ~^ pY 


5.8±0.3{stat.) ±0.5{sys.) 
2.4 ± 0.2{stat.)+°ol{sys.) ± 0.3(moc/.) 


50-100 
85-105 


0.025-0.5 
0.4-1.2 


ZEUS 
ZEUS 


[68] 
[76] 


7P (pp 


7.3± 1.0(stat.) ±0.8(syst.) 
6.3 ± 0.7(stat) ± 0.6(sys.) ± 0.3(morf.) 


60-80 
85-105 


0.1-0.5 
0.4-1.2 


ZEUS 
ZEUS 


[278] 
[76] 


'-fp —>■ (pY 


2.1 ± 0.5{stat.) ± 0.3{sys.) ± 0.4(morf.) 


85-105 


0.4-1.2 


ZEUS 


[76] 


7P J / ipp 


4.73±0.25(stol)lJ^J^(sys.) 
4.99±0.13(stol)±°-^5(sys.) 
4.15±0.05(stat)+°:?°(s2/s) 
4.0 ± 1.2(stat.)+?:J(s|/s)l|^:e(mod.) 


40-150 
40-150 
20-290 
85-105 


0-1.2 
0.07-0.9 

0-1.2 
0.4-1.2 


HI 

HI 
ZEUS 
ZEUS 


[70] 
[73] 
[71] 
[76] 


IP ~^ j/ipy 


0.7 ± OA{stat.) ± 0.2(si/s.)„o:3(''^oa.) 
1.07±0.03(stat.) ±0.11(si/s.) 


85-105 
40-150 


0.4-1.2 
0.15-3 


ZEUS 
HI 


[76] 
[73] 


-fP^^{2S)p 


4.69 ± 0.57(stot.) ± 0.46 


40-150 


0.07-0.9 


HI 


[73] 


7P ^ ^{2S)Y 


0.59±0.13(stat.) ±0.12(sys.) 


40-150 


0.15-3 


HI 


[73] 



Table 4: Recent HERA measurements of the b slope for PHP elastic and proton- 
dissociative vector-meson production. 



was predicted in [26] on the basis of node effect [18]. The HI results for both elastic and 
proton-dissociative ^(25') production [73] are consistent with this prediction. Within the 
experimental error bars Ab = h^^ — h'^'^^^ ^ (4.5 -r- 5) GeV~^ is flavor independent which 
agrees perfectly with the theoretical expectation (56), (194) and must be contrasted to a 
strong flavor dependence of hei- 

Notice a very small hdissilP J/'^Y), which is consistent with the equally small slope for 
the double-dissociative hadronic diffraction pp — > XY, for the compilation of the hadronic 
data from the fixed target to ISR energies see [122, 123]. 



8.2.3 Experimental results: Q^-dependence of the diffraction 
slope 

Figure 60 shows the Q^-dependence of the diffraction slopes fcp(Q^) and bj/^{Q^). In the 
case of the p-meson, a strong dependence of the bp is observed: b ^ 11 GeV~^ in the 
photoproduction hmit and drops down to 4 5 GeV"^ in the hard electroproduction. 
In striking contrast, in the case of J/ijj production the Q^- dependence is very weak or 
absent, bj/^ ~ 4.5 GeV~^ was found throughout the entire range. 

Figure 60 shows also the comparison between the elastic and proton-dissociative slopes 
in p meson production. One clearly sees an approximately Q^-independent difference 
Ab = b^p{Q'^) — b'^p^^ 4 5 GeV~^ between the two data sets, which coincides with 
the photoproduction value shown in Tab. 4. This is an important confirmation of the 
theoretical expectation (194). 

Figure 61 shows the above results for elastic slopes &p (Q^) and &j/^(Q^) plotted as a 
function of (Q^ + ^y)- The strong flavor dependence obvious if one compares the pho- 
toproduction values bp\Q'^ = 0) at bj^^{Q'^ = 0) is dramatically reduced in the variable 
(Q^ + my). The theoretical curves are from the /c^-factorization model [55]. Because of 

2 

flavor dependent departure of the hard scale Q from the simple estimate (9), the theo- 
retical predictions for the J/ip are slightly lower than for the p. Similar results were found 
earlier in the color dipole model [26,301]. The origin of the weak Q^-dependence for the 
J/\l/ is in a very small contribution to &j/^(Q^) from b^^viQ"^) oc r| oc 1/(Q^ + my), see 
Eqs. (80), (193). 

The similar approximate {Q"^ + my) scaling holds for the proton-dissociative reaction: the 
results for 6^*^* for the largest values of in Fig. 60 are perfectly consistent with the 
photoproduction value for 6j*f^ from Tab. 4. 
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Figure 60: The preliminary ZEUS data on diffraction slope b for elastic and 
proton- dissociative production of the p [84]) upper plot. Bottom plot shows the 
recent ZEUS [90] and HI [89] data for elastic production of the J /ip mesons as a 
function ofQ^. 
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Figure 61: The diffraction slope b of elastic production of the p (HI: [16], ZEUS: 
[68,81], the data from [84] are not shown) and J/ip (HI: [89], ZEUS: [90]) mesons 
as a function of + My. The predictions of from the kt-factorization approach 
[54, 55] are shown. 

8.3 Shrinkage of the diffraction cone and the Pomeron trajec- 
tory 

8.3.1 The VK-dependence in photoproduction 

Figures 62 and 63 sliow tlie energy dependence of tlie diffraction slope b for the p and J/ ip 
pfiotoproduction, respectively. One sees a steady growth of the diffraction slope, i.e., the 
shrinkage of diffractive cone. When parameterized in terms of (37), the p photoproduction 
data yield a'jp = 0.23±0.15{stat.) +Q lj{syst.) GeV"^ [68]. This result for the soft reaction 
is consistent with a'p = 0.25 GeV~^ found in the simple Regge pole description of elastic 
pp scattering. The J/"^ photoproduction must be regarded as a hard reaction, here the 
evidence for nonvanishing a'jp is even stronger. Fitting to the form b = bQ{90GeV) + 
Aa'jplni^) gave a'p = 0.116 ± 0.026(stat.)+°;°^°(s|/st.) GeV-^ [71]. 

8.3.2 The Pomeron trajectory 

In addition to fitting the differential cross section da/dt at given W and evaluating the 
diffraction slope 6 as a function of W, one can study the W^-dependence of the differential 
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Figure 62: The diffraction slope bp as a function of W for p photoproduction 
[69]. The continuous line shows the result of the fit of the form (37) to the recent 
HI measurement, other HERA and fixed target measurements are shown for com- 
parison; the extrapolation of the fit to the low W region is indicated by the dashed 
line. 
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Figure 63: The slope bj/^ as a function ofW for J /ip photoproduction as measured 
by ZEUS [71]. The results of the kt-factorization calculations [54,55] based on the 
oscillator (solid line) and Coulomb (dashed line) wave functions are compared with 
the data. 
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cross sections at fixed t. According to Eq. (28) the ly-dependence of tlie differential cross 
section is oc H^^('^Jp(*)~i), so that one can extract ajp versus t , i.e. measure the effective 
trajectory ajpit) of the t-channel vacuum exchange. 
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Figure 64: Determination of the Pomeron trajectory from the reactions p^p 
and '-•IP (pp. The dots are results of ZEUS measurements [76]. The solid 
lines are results of the linear fit. The DL parameterization for the soft Pomeron 
trajectory [108] is shown for comparison as a dashed line. 

Figures 64, 65 demonstrate the ajpit) measured by ZEUS in the photoproduction of the 
p, and J mesons. The obtained parameters of the effective Pomeron trajectories are: 

aip{p- 1) = (1.096 ± 0.021) + (0.125 ± 0.038) ■ t ; 

ajp{^; t) = (1.081 ± 0.010) + (0.158 ± 0.028) ■ t ; (196) 
aip{J/^; t) = (1.200 ± 0.009 ± 1!]:[!;^) + (0.115 ± 0.018 ± tol^oil) " t ■ 

which must be compared with the DL parameterization (50). For all the vector mesons 
the determined slope of the vacuum trajectory is significantly non-zero. Within the error 
bars, the value of a'jp appears to be insensitive to the particular type of the vector meson. 

The dependence of the extracted Pomeron trajectory on has been studied by ZEUS for 
the case of J/\& production. Figure 66 shows comparison between the Pomeron trajectory 
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Figure 65: Pomeron trajectory as obtained from J/ip photoproduction measure- 
ments by ZEUS [71]. The results of HI measurements [70] are shown for compar- 
ison. Linear fits to the ZEUS and HI data are shown and compared with the DL 
parameterization for the soft Pomeron trajectory [108]. Bottom. The predictions 
from the kf-factorization approach [54, 55] are compared with the ZEUS data. 
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Figure 66: A comparison of the Pomeron trajectory extracted form the ZEUS 
data on real photoproduction and electroproduction of the J/\l/ [71, 90]. The solid 
lines are results of the linear fit. The Donnachie-LandshofJ parameterization [108] 
of the soft Pomeron trajectory is shown for comparison as a dashed line. 
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at = and = 6.8 GeV^. The A;_|_-factorization predicts a slight increase of the 
effective intercept with Q^, see Fig. 46. Even if O-'^fkl parameterization (129) were 

a constant — the solutions of the color dipole BFKL equation [248] for the diffraction slope 
exhibit weak /t^-dependence of (y'sFKLy though [25], — the above described extraction will 
yield a weakly Q^- and flavor-dependence values of a'jp. Such a non-Regge effects in the 
extracted a'jp have been evaluated in [55,299] and demonstrated in Fig. 67. Because of 
the same non-Regge effects, the theoretical results for the effective vacuum trajectory are 
sensitive to the wave function of the vector meson. The effect can be seen from at the 
bottom plot of Fig. 65, it is negligible for all the practical purposes. 
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Figure 67: The kt-factorization model predictions [54, 55] for the Q'^ -dependence 
of the slope of the effective Pomeron trajectory, a'p, for the p [68, 76] and J/ip 
[71, 90] production. The results for a'jp obtained from the study of cone shrinkage 
h{W) and from the trajectory analysis a]p{t) are shown separately by open and filled 
dots, correspondingly. 



8.4 The gluon-probed radius of the proton and the Pomeron 
exchange radius 

The above presented experimental data on the diffraction slope for the J/"^ production 
give 

Bn = b2G ^ (4 ^ 4.5) GeV-2 . (197) 

It is the non-perturbative parameter which defines the form factor of the proton probed 
by color-singlet two-gluon current, G2g(^) ~ exp(i62G^)- K must be compared to the elec- 
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tromagnetic form factor Gem{t), the familiar dipole parameterization for which, Gem{t) 
1/(A2 - t)2 with PS 0.7 GeV^ amounts to 



em 



4 

A2 



(5.5 ^ 6) GeV 



-2 



(198) 



The departure of 62G from B^m niust be regarded as substantial one. The photon only 
couples to charged partons in the proton, whereas the gluonic form factor is sensitive to 
the distribution of all color charges - the (anti) quarks and gluons - in the proton. Should 
one interpret the observed inequality 62G < Bem. as an evidence for the gluonic lump in 
the center of the proton? 

The same problem is present in irp^pp scattering. For instance, Schiz et al. find that 
at 200 GeV the observed t-dependence of scattering amplitudes is very well reproduced 
by the product of the charge form factors of the beam and target particles [300]. One 
would interpret that as an equality of the strong interaction and electromagnetic radii, 
62G ~ -Bern; thls particular energy. However, in hadronic scattering the shrinkage of the 
diffraction cone is rapid, and at lower energies ~ 10 GeV, still in the applicability domain 
of the Regge model, both and hpp are substantially lower than at 200 GeV, and one 
runs into the same problem with 62G < -Bem- 

A very small value of the gluon radius of the proton, as a matter of fact, 62G ^ -Bem? 
has been found by Braun et al. from the QCD sum rules [302]. The QCD sum rule 
results depend on the interpolating field for the nucleon, though. Braun et al. take the 
qqqg operator, whether the radius evaluated for the hybrid higher Fock component of the 
proton applies to the whole proton or not remains an open issue. 

Another small parameter which emerges in small-t diffraction is the Pomeron exchange 
radius which is probed either in the double double proton-dissociative pp scattering or in 
proton-dissociative vector meson production at small scanning radius rs such as 7p — > 
J/^Y of the p production in large-Q^ DIS, 7*^ — > pY: 



The latter two reactions have a hard scale on the vector meson side, and one may link 
the anomalously small bjp to the small propagation radius of perturbative gluons -Rc ~ 1 
GeV~^. The pp double diffraction reaction probes the soft Pomeron, why do the soft and 
hard Pomeron exchanges have equally small exchange radius is an open issue. 



bp ^ h{pp XY) ~ h{^p J/^Y) 



b{-f*p ^ pF) ~ 1 GeV 



-2 



(199) 
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8.5 Beyond the diffraction cone: large |^| as a hard scale 

8.5.1 Large-t vector mesons as a Mueller-Navelet isolation of 
the hard BFKL exchange 

The basics of the pQCD treatment, and hehcity properties, of large-t vector meson produc- 
tion were reviewed in Sections 4.10 and 5.4. Here we focus on the t- and ly ^-dependence 
of the large-t cross section. 

On the theoretical side the large-|t| production of vector mesons is a very promising testing 
ground for ideas on the BFKL Pomeron, because the large momentum transfer A flows 
along the BFKL Pomeron from the target to the projectile. There is a close analogy to 
the long sought Mueller-Navelet isolation of the hard BFKL exchange by selecting DIS 
events with one hard jet in the target fragmentation region and the second hard jet in 
the photon fragmentation region [303]. Within the effective parton description of the 
proton-dissociative 7p — > VY in terms of the elastic scattering 7g Vq' the recoil quark 
(antiquark, gluon) with large transverse momentum A gives rise to the target hard jet 
of Mueller & Navelet, whereas the large-t vector meson is a substitute for the Mueller- 
Navelet hard jet in the photon fragmentation region. The vector meson production is 
even more advantageous because one has an access to larger rapidity gaps between the 
vector meson and the system Y than it would be possible in the Mueller-Navelet two-jet 
process. 

The pQCD two-gluon approximation misses the dependence on the rapidity gap and the 
total normalization must be adjusted to the experimental cross section. Otherwise, as we 
shall see below, it is doing a reasonable job on the t-dependence. The state of the art 
BFKL based calculations use Lipatov's leading order fixed-af approximation for the 
unintegrated off-forward gluon density. By the logics of the calculation one may expect 

ar^'^-c^dtl) (200) 

and expect the Regge energy dependence (154) with the BFKL trajectory given by 
Eq. (59). 

8.5.2 Theoretical expectations for flavor dependence at large-t 

A crude reinterpretation of very involved theoretical calculations [46,53,261,304,305] 
starting with the pQCD subprocess 7*g — Vq' is as follows: 

First, the relationship between the cross section of the theoretical partonic subprocess 
7g — > Vq' and the experimentally observed 7p — > VY involves the effective number of 
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partons in the proton, 



1 



N,{t)= I dx\x')'^--^'^^'U^g{x',t) + y][q{x',t)+q{x',\t\)]\, (201) 




where the Regge dependence on x' is reabsorbed into the flux of equivalent partons. After 
Np{t) is factored out, one obtains the cross section of the partonic subprocess at a fixed 
energy W^g-. 



da{'yq — ^ Vq') 



dt 



MlP^VYJ^ (202) 



The cut Xmin = 0.01 used by ZEUS collaboration [77] gives Np{t) shown in Fig. 68. HI 
collaboration [79] imposes the cut My < M^„^ = 900 GeV^, which translates into 





t 




M2 

max 


+ 


t\ 



•^min , ^-o I Ul ■ (203) 



Evidently, for ajp{t) = 1 the number of partons diverg and at x^m = 0.01 

it exhibits rise with \t\ because of the scaling violations in the gluon density, see Fig. 68. 
In the opposite to that, for the HI cut (203) the decrease of Np with |t| is driven by the 
rise of Xmin, see Fig. 68. The sensitivity of Np to the trajectory apit) is very strong: for 
Ajp = ajp(t) — 1 = it is a true number of partons, is very large and is a steep function of 
Xmin, for Ajp = 0.5 the integral (201) is reminiscent of the momentum sum rule integral 
and yields weakly t-dependent Np ^ 2.5. see Fig. 68. 

Second, in view of an approximate SCHC with dominant ar the cross section must be 
proportional to myT(y — >• e~^e~). That does not exhaust the flavor dependence because 
the onset of the hard regime does obviously depend on the mass of the heavy quark. The 
pQCD two-gluon exchange approximation suggests that for slow Fermi motion in vector 
mesons the appropriate hard scale is 

g'^(m^ + |t|) (204) 

( [261,305] and references therein). Although the Fermi motion can change the coefficient 
in front of rriy, one must conclude that for the J/\l/ the large t means |t| ^ "^j/t/i ~ 
10 GeV^. The numerical studies by Poludniowski et al. [261] show that even for light 
mesons the variation of the constituent quark mass from nip/ 4 to mp/2 to rrtp changes 
the predicted cross section by a factor of ~ 2 even at |t| as large as 10 GeV^. see Fig. 69. 

Third, in the pQCD two-gluon approximation the target quark is regarded as pointlike 

one and the t dependence is entirely due to the 7 — s> transition vertex with the hard 
2 

scale Qf-. Then, upon the x-integration in (151), 

da2G{.iq V(l') 2m\ 2n,|N":4r(^:l^e^e^ /ont;^ 
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X . =ltl/(ltl+900) 

min ^ ' 
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Figure 68: The t-dependence of the effective number of partons Np{t) for the 
ZEUS (top box) and the HI (bottom box) kinematical cuts for different values of 
the intercept a = 1 + A [306]. Shown also are the curves for the t-dependent 
trajectory aip{t) parameterized by Eq. 209 and shown in Fig. 71. 
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\t\ range 
(GeV^) 


{\t\) 
( GeV2) 


6 


a{t) 


2-5 


3.06 


0.77 ±0.14 ±0.10 


1.193 ±0.035 ±0.025 


5-10 


6.93 


1.29 ±0.23 ±0.16 


1.322 ±0.057 ±0.040 


10-30 


16.5 


1.28 ±0.39 ±0.36 


1.322 ±0.097 ±0.090 



Table 5: The value of 6 obtained when applying a fit to the data of the form 
a{W) oc for each \t\ range, together with the corresponding vacuum trajectory 
a{t) obtained from a{t) = [6 + 4)/4. The first uncertainty is statistical and the 
second is systematic. 

where we indicated the natural choice of the running couplings of gluons to the target 

2 

quark, a^dtl), and to quarks in the vector meson, as{Qt)- Based on the experimental 
data on vector meson decays [102] the prediction for flavor dependence of dap"^^ / Np{t) at 

2 

identical values of the hard scale is 



p : uj : (f) : J/ip 



0.8 X 



9 



2.1 X 



9 



56 X 



8 
9 



(206) 



Recall that in the studies of the dependence one had to compare the cross sections at 
at identical {Q^ + my). 

Fourth, in the BFKL approximation the target quark becomes effectively non-pointlike one 
and introduces the approximately flavor independent factor oc l/|t| to the cross section, 
whereas the t-dependence from the 7 — > ^ transition vertex will be weaker: 



da 



BFKL 



iiq ^ Vq') 



dt 



e 




(207) 



^2 



Consequently, the flavor dependence (206) at identical values of the hard scale must 
be tested for \t\N-\t){dap'' /dt) rather than dap''/dt. Notice an extra suppression 

oc l/((5t)^^^^^^^ coming from the Regge parameter. To run the strong couplings in (207) 
is to go beyond the accuracy of the scaling BFKL approximation. 



8.5.3 The experimental results: measuring the trajectory of the 
hard BFKL Pomeron 

The VT-dependence of the proton-dissociative vector meson production has been measured 
by both ZEUS and HI collaborations. 

The absolute W dependence of the J/ip production has been measured by HI Collab. [79]. 
The experimental data are shown in Fig. 70, the results of the Regge fits are presented 
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[GeV] 

Figure 70: The HI results for the J/\l/ photoproduction cross section as a function 
of W^p in three bins of \t\ [79]. The inner error bars correspond to the statistical 
error and the outer error bars are the statistical and systematic errors added in 
quadrature. The solid lines show the predictions from the BFKL model [307], the 
dashed-dotted curve is the result from the DGLAP model [308]. 
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Figure 71: A compilation of the ZEUS and HI results on the determination of the 
Pomeron trajectory aip(t) (A{t) = a]p{t) — 1) from the J/\l/ photoproduction. The 
solid curve shows a possible interpolation from the regime of shrinking diffraction 
cone at small t to the hard BFKL regime at large t. 

in Tab. 5. Tliese results from HI give a solid evidence for ajpit) > 1 at large The 
found values of Ajp = aip(t) — 1 are close to the leading order BFKL prediction with 
(0.1-^0.15), a comparison with the NLO intercept is unwarranted at the moment. 
On the other hand, the large-t extrapolation of the ZEUS results (197) suggests aip(t) < 1 
for |t| > 2 GeV^. (For the p and mesons the similar crossover takes place at \t\ ~ 1 
GeV^, one would readily attribute that to the process being still soft.) Now recall that 
the shrinkage at small-t is driven by the infrared growth of as by which the low-t BFKL 
evolution becomes sensitive to the infrared region around the finite propagation radius Rc 
for perturbative gluons. In contrast to that, in the Mueller-Navelet large-t regime the large 
momentum transfer ~ A flows through propagators of all t-channel gluons, the infrared 
contribution will be suppressed and gross features of the fixed a^^^^ ~ oisif), leading 
log^, BFKL evolution will be recovered. One can fancy the nonlinear dependence of 
the vacuum trajectory of the form 



shown in Fig. 71 for = 2.5 GeV^. 

The Ansatz (208) for the Pomeron trajectory turns over at \t\ ~ 1-^2 GeV^, which is close 
to the natural scale ~ (0.5 ^ 1) GeV^. The change of the sign of the derivative a'{t) 
from small to large t is supported by the ZEUS experimental data shown in Figs. 72 and 



a{t) = 1.2 + (0.16 GeV-^)t ■ 



A6 



(208) 



t|3+A6 



+ 0.16 



t^ + AS 
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73. Although the ly-dependence of the efficiency of the photoproduction tagger hinders 
the direct measurement of the absolute ly-dependence and the determination of aip{t) 
with the present ZEUS data [77], the large-t slope of the vacuum exchange trajectory can 
be measured in the tagger independent manner, 

,'tiw'/f cx W-'('-« (209) 

The experimental data from ZEUS for this ratio and the found values of the slope of the 
vacuum trajectory a' are shown in Fig. 72. The summary of the low-t and high-t results 
for a' is presented in Fig. 73, where we also show a'{t) for the parameterization (208). 

8.5.4 The experimental results: the ^-dependence for a nucleon 
target 

The t-dependence of the p, and J/ ip meson proton-dissociative production cross section 
at high-t is shown in Figs. 74 and 75. It is much slower than the exponential one typical 
of the diffiaction cone and is in broad agreement with the inverse power law as discussed 
in Section 8.2. 

The quantitative interpretation of the experimental data taken at the moderately large t 
depends on the choice of the hard scale. A fit to the observed t-dependence in the form 
da/d\t\ oc \t\~"'^ yielded for the ZEUS data [77] the exponents 

Up = 3.21±0M{stat.)±0.15{syst.) (ZEUS, 1.2 < |t| < 10 GeV^) ; (210) 
= 2.7±0.1{stat.)±0.2{syst.) (ZEUS, 1.2 < |t| < 6.5 GeV^) ; (211) 
nj/^ = 1.7±0.2{stat.)±0.3{syst.) (ZEUS, 1.2 < |t| < 6.5 GeV^) . (212) 

which must be compared to ~ 4 expected from theory. The values of |t| < 6.5 GeV^ 
in the ZEUS data on the J/\& production are arguably too small for the onset of the true 
large-t behaviour. The J/"^ production data from HI [79] extend to \t\ < 30 GeV^, see 
fig. 75, and give the exponent 

nj/^ = 3.00±0.08(stal) ±0.05(s?/sl) {HI, 4.46 < |t| < 30 GeV^) . (213) 



A comparison of the results (212) and (213) shows an importance of the finite mass 
effects in the t-dependence, see also Fig. 69. The theoretical calculations by Poludniowski 
et al. [261] within the scaling BFKL approximation are shown in Figs. 75.76, and 
clearly show an improvement from the pQCD two-gluon to BFKL approximation. The 
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Figure 72: The da/dt ratios for p(left) and (p (right) production cross sections as 
a function of W in five(four) t intervals. The lines represent the result of the fit 
with Eq. 208. 
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Figure 73: HERA results on a' for the elastic and proton-dissociative (at \t\ > 
1.3 GeV^^ vector meson production [68, 71, 76, 77, 79] compared with the a'{t) 
from the parameterization (208). The reference value for soft hadronic interactions 
^'soft ~ 0-25 GeV~^ is shown as a dashed line. The points are put in the center of 
\t\ range in which a' is measured. The vertical inner bars indicate the statistical 
uncertainty and the outer bars represent the statistical and systematic uncertainties 
added in quadrature. 
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Figure 74: The ZEUS results for the differential cross sections da^p^vY / d\t\ in 
the range 80 < < 120 GeV and x > 0.01 for p°, and J/ip [77] production. 
The lines are results of the fit to the data with the function t)"". The shaded 
bands represent the correlated uncertainties due to the modeling of the hadronic 
system Y. 
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Figure 75: The differential cross section da^'-fp VY)/dt for the J/ip photopro- 
duction at large t m the range 80 < < 120 GeV and x' > |t|/(900(GeV2) + \t\) 
from HI collaboration [79]. The theoretical results for the pQCD two-gluon ex- 
change with fixed and running as and scaling BFKL approximations are from 
Poludniowski et al. [261]. 

impact of tlie running strong coupling on tlie pQCD two-gluon results for the t-dependence 
is substantial, the BFKL calculations are for the fixed coupling. The results for the 
and J/\l/ are based on the parameters of the model which were adjusted to the p 
photoproduction. 



8.5.5 The experimental results: the flavor and t-dependence for 
a partonic subprocess q'Y 

The most direct test of the BFKL approach to large-t vector mesons is provided by 
the reanalysis [309] of the experimental data in terms of the cross section of partonic 
subprocess 



rfa(7g Vq') 



dt 



"^'^'"-""'oc.^^, (214) 



Np{t) dt (|t|+m^)"^ 



Here the residual t-dependence in the rhs probes the true dynamics of hard '-/V transition, 
see Eq. (207). In such a representation the scale- invariant BFKL approximation predicts 
the flavor independent 

ny = 3 + 2ABFi^L ~ (3.2^3.6). (215) 

The major problem with the extraction of the partonic cross section is that the absolute 
value and t-dependence of the number of effective partons Np{t) exhibits a strong sensi- 
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Figure 76: The differential cross section da^'-fp —* VY)/dt for the p (upper plot) 
and (j) (lower plot) photoproduction at large t in the range 80 < < 120 GeV 
and x' > 0.01 from ZEUS collaboration [77]. The theoretical results for the pQCD 
two-gluon exchange with fixed and running as and scaling BFKL approximations 
are from Poludniowski et al. [261]. 
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tivity the the pomeron trajectory, see Fig. 68. The resulting uncertainty propagates into 
the magnitude and t-dependence of the cross section of |t|(i(T(7g — > Vq')/dt in (214) and, 
consequently, into the determination of the exponents ny from the fit of the partonic cross 
section to the parameterization (214). 

For the consistency with direct experimental measurements of the pomeron trajectory one 
must use the parameterization (208) which correctly reproduces all the features of the HI 
and ZEUS data shown in Figs. 71 and 73. The resulting fits to the ZEUS data [77] on the 
p,0 and J/^ production at (W^p) = lOOGeV yield Up = 2.08 ± 0.06, = 1.83 ± 0.13, 
nj/^{ZEUS) = 0.78±0.64. In the HI data [79] on the J/^ production the energy {W^pit)) 
slightly rises with t which introduces a certain bias into the t-dependence and enhances 
nj/^. Neglecting that bias and excluding the point at lowest t yields nj/^{Hl) = 2.55±0.2. 

The above cited error bars do not include the theoretical uncertainties of Np{t) connected 
to the parameterization of the t-dependence of aip{t). For the sake of illustration, we cite 
here the results found if Np(t) is evaluated for fixed aip{t) = 1 + Abfkl ~ 1.25, although 
such a flat aip{t) is inconsistent with the HI and ZEUS data shown in Figs. 71 and 73. 
In this case the t-dependence of Np{t) will be much weaker for both the ZEUS and HI 
cuts, see Fig. 68, it doesn't change substantially with the further increase of Abfkl- The 
partonic cross sections extracted from the same data will have much steeper t-dependence 
and the fitted exponents ny will be substantially larger than for the pomeron trajectory 
of Eq. (208): Up = 2.86 ± 0.05, = 2.66 ± 0.12, nj/^{ZEUS) = 3.88 ± 0.62 ~ 3 for the 
ZEUS data and n,j/^{H\) 3.86 ± 0.26 for the HI data. The statistical error bars in fitted 
values of ny for the two choices of Np{t) are misleading because the model dependence in 
the extraction of the exponent ny of the t-dependence is much larger than the statistical 
error bars. In their scale-invariant BFKL calculations shown in Fig. 75, Poludniowski et 
al. [261] use fixed as = 0.25, which amounts to even larger Abfkl = 0.66 and ny 4.3. 
Within those uncertainties, the data on different vector mesons do not exclude the flavor 
independent ny and the observed t-dependence does not conflict the BFKL expectation 
(215). 

A very large a]p{t) = 1 + Abfkl = 1-66 used in the theoretical calculations [261] conflicts 
the experimental data on the Pomeron trajectory shown in Figs. 71 and 73. Furthermore, 
the effective number of partons Np{t) evaluated with Abfkl = 0.66 is by the factor 
~(5-10) smaller than for the experimentally suggested trajectory shown in Fig. 71. This 
uncertainty is not discussed in [261] and casts a shadow on the agreement between 
the theory and experiment in the magnitude of the cross section. Before drawing firm 
conclusions on the status of the BFKL approach one needs much better understanding 
of the Pomeron trajectory and incorporation of the realistic Pomeron trajectory into the 
theoretical formalism. 

The same model-dependence of extraction of the partonic observable \t\da{'~fq — * Vq') /dt = 
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Figure 77: (a) The ratio of the to p° cross sections as a function of —t or 
. The (p/p^ results as a function of —t for proton-dissociative photoproduction 
from this analysis are shown with solid circles and those from the ZEUS 1995 [76] 
measurement with the solid squares. The shaded bands represent the size of the cor- 
related uncertainties due to the modeling of the dissociative system, Y. Open trian- 
gles atQ^ ^0 GeV^ [278], = rCeV^ [86] and = 12.3GeV2 [280] represent 
the (f)/p^ ratio of the elastic cross sections as a function ofQ^ from ZEUS, while the 
open squares represent those from HI [87]. (b) The ratio of the J/ip to cross sec- 
tions as a function of —t or Q^. The same convention for symbols as for (p/ p^ ratio 
is used. Open triangles at Q"^ k, OGeV^ [282] and Q"^ = 3.5, 13GeV^ [81] represent 
the ZEUS measurements, while the open squares represent those of HI [16,89]. The 
dashed lines correspond to the SU(4) predictions, while the dotted and dashed-dotted 
correspond to the pQCD values given by Eqs. 206 and 178, respectively. 
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\t\{da{'-fp — > VY)/dt)/Np{t) affects the discussion of the ffavor- dependence of large t cross 
sections. First, the point that neither nor \t\ are the correct hard scales to compare dif- 
ferent vector mesons is illustrated by Fig. 77. Whereas for light meson the t-dependence, 
as well as the Q^-dependence of the ratio a^/ap is weak, the ratio crjj^/ap changes by 
more than two orders of magnitude. The apparent approach to the SU(4) ratios at largest 
measured value |t| is misleading - no true asymptotics can be reached at \t\ < rn'^j/^- 

Strong departure from the SU (4) ratios is evident from Fig. 78, where we show the ZEUS 
and HI data in the form of the partonic subprocess observable \t\da{'~fq — > Vq')/dt = 
\t\{da{jp — * VY) / dt) / Np{t) plotted as a function of {\t\+my). Although the experimental 
data on the p, and J/\l/ vector mesons don't have an overlap in {\t\ +my), it seems safe 
to extrapolate the ZEUS data on the p and production to + = 12.5 GeV^ 

typical of the ZEUS data on the J/ip production. The so extrapolated p and cross 
sections have the factor of ~ 2 uncertainty. 

We recite from (206) the flavor dependence of {1 / Np{t)){day^^ / dt) at identical (Q^ 
based on the vector meson decay properties: 



p : a; : : J/?/^ = 1 : 0.8 X ^ : 2.1 X ^ : 56 X ^ . (216) 

If the number of effective partons Np{t) is evaluated with the pomeron trajectory (208) 
which correctly reproduces the HI and ZEUS data shown in Figs. 71 and 73, then the 
extrapolation of the p and cross sections to {\t\ + rri^j/^) = 12.5 GeV^ gives the cross 
section ratios (within the factor ~ 2 extrapolation uncertainty) 

1 2 1 8 

p : : J/^p =1 : - x 2.1 x - : — x 56 x - . (217) 

2 9 15 9 

If Np{t) is evaluated for Abfkl = 0.25 = const, then the same extrapolation gives slightly 
different cross section ratios 

p : : JM ^ 1 : - x 2.1 x - : - x 56 x - (218) 
H /y- 3 9 7 9 ^ ' 

For the both choices of the Np{t) the principal effect is an enhancement of the light vector 
meson production with respect to the J/\l/ production. Such an enhancement due to the 
chiral-odd ^qq transitions [242] is present in calculations of Poludniowski et al [261]. 
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Figure 78: The t dependence of the large-t vector meson production measured 
by ZEUS [77] and HI [79] presented as a differential cross section of the partonic 
subprocess \t\da{'yq —>■ Vq')/dt = \t\{da{'jp — > VY) / dt) / Np(t) . The straight lines 
are results of the fit by Eq. 214, the lowest \t\ data point of HI has been excluded 
from the fit. The left box shows the results for the number of effective partons Np{t) 
evaluated for the Pomeron trajectory a]p{t) of Eq. (208) shown in Fig. 71, the right 
box is for Np{t) evaluated for the fixed Ajp = 0.25. 
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9 Summary and conclusions 



Slightly more than a decade ago, quite unexpectedly, HERA has become a unique facility 
for exploring the diffractive physics in an entirely new domain of many different hard 
scales. Compared to the fixed target data, the center of mass energy W, the explored 
regions of and t were extended by one order of magnitude. By now the principal 
features of the flavor-, t-, Q^- and H^^-dependences of the observed cross sections, polar- 
ization properties of produced vector mesons are well established. The high statistics of 
the data from HERA allowed for the first ovservation of SCHNC in high energy diffrac- 
tion. Regarding the experimental situation, there is an overall consistency between the 
experimental data from HI and ZEUS Collaborations and the early data from fixed target 
experiments. 

What did change in our understanding of high-energy diffractive scattering after that 
decade of amassing high-precision experimental data by HI and ZEUS Collaborations at 
HERA? What was the impact of these data on theoretical ideas on high-energy vacuum 
exchange in the t-channel - the Pomeron? How strongly is the post-HERA pQCD Pomeron 
different from the pre-HERA Pomeron approximated by an isolated Regge pole with an 
intercept a^o/t ~ 1.1? Is our theoretical understanding of diffraction sufficient to get all 
the information we can, and would like to, from the available experimental data? 

The scaling violation in inclusive DIS is a classical short-distance phenomenon dominated 
by scales ~ l/Q [124-126]. Thanks to the 1970's groundbreaking works by Fadin, Kuraev, 
Lipatov and Balitsky [35-37] the theory was well prepared to the observed departure from 
the DGLAP evolution and the BFKL Pomeron reinterpretation of a steep small-x rise of 
structure functions discovered at HERA. 

The principal virtue of diffractive vector meson production in DIS at HERA is that it is a 
short distance dominated process. At small to moderate t within the diffraction cone, the 
short- distance property is quantified by a scanning radius [17-20], at large-t the short 
scale is set by 1/ ^/\t\- The early discussion of importance of large-t photoproduction of 
vector mesons as a testing ground of ideas on the pQCD Pomeron goes back to mid-80's 
work by Ginzburg, Panfil and Serbo [30]. The idea of unified color dipole description of 
inclusive DIS and diffractive vector mesons, the concept of the scanning radius [17-20] 
and the importance of Q = ^iQ"^ + rriy) as a hard scale at small t [21,22], and the 
possibility of vector meson production as a testing ground for models of gluon density in 
the proton emerged in early 90's in works by Kopeliovich, Nemchik, Zakharov et al. [20] 
and Ryskin [22] and Brodsky et al. [23]. Ever since then the collective effort by many 
groups has led to a fairly refined treatment of vector meson production at small to large 
|t| and to an understanding of an accuracy and limitations of the leading log^ pQCD 
approaches. The wealth of the experimental data collected by HI and ZEUS confirmed 
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all the gross features of the pQCD based description of the process. 

An approximate (Q^ + my )-scaling of all observables - total cross section [21], the diffrac- 
tion slope [25], the exponent of the energy dependence [21] — alias the Pomeron intercept 
— has been the recurrent theme in our discussion of the experimental data, its experi- 
mental confirmation at HERA must be regarded as a major discovery and an undoubted 
success of the pQCD approaches. The pQCD dictated relationship between the energy 
dependence of vector meson production and inclusive DIS encoded in terms of the energy 
dependence of the gluon density [20, 22, 23] has been confirmed experimentally beyond 
reasonable doubt and is still another major discovery at HERA. 

Success with theoretical predictions of the absolute cross sections is modest one. The 
(Q^ + TTiy) as a hard scale for small-t diffractive vector mesons is entirely analogous to 
as a hard scale in inclusive DIS. The pQCD only predicts the dependence on those 
hard scales starting from certain soft input. Within the color dipole and /c^-factorization 
approaches this input is universal for vector meson production and inclusive DIS, still it 
does not come from first principles of pQCD. Only one leg of the pomeron, which couples 
to the 'j*V transition, rests on the hard pQCD ground, the second leg which couples to the 
proton is always in the soft regime. Once the fc^-factorization and other related pQCD 
model predictions are normalized to the J/\i/ photoproduction data, the description of 
the observed Q'^- and energy dependence of <7j/q,{Q'^) is close to a perfect one, which must 
be regarded as a real success of pQCD in the domain of hard diffractive scattering. 

Nonetheless a factor of ~ 2 sensitivity of pQCD model predictions to the wave function of 
vector mesons is a well established limitation of the leading order log ^ formalism and can 
not be eliminated at the moment. The ratio of longitudinal and transverse cross sections, 
Rv = o'l/o't, is an example of an observable which exhibits especially strong sensitivity 
to the wave function of light vector vector mesons [253]. A very good demonstration of 
the sensitivity to the wave function is a node effect in the \E'(2S') production [17-19] which 
suppresses the cross section, makes it grow with energy faster than the J/\E' production and 
leads to counterintuitive inequality of diffraction slopes for the "^{23) and J/\& production 
[26]. The further tests of pQCD predictions for diffractive vector meson production call 
upon the development of NLO /c^-factorization formalism. This includes the theoretical 
understanding of the effect of higher Fock states in vector mesons and the derivation of 
the /c_L-factorization impact factors for the 'y*V transition to NLO in log^. 

The works by D.Ivanov and Kirschner [212] and Kuraev et al. [170], in conjunction with 
B.Zakharov's [211] early pQCD motivated discussion of helicity flip in hadronic scattering, 
have led to understanding of SCHNC as a generic property of high-energy scattering the 
origin of which does not require an applicability of pQCD. The emerging phenomenology 
of spin properties of diffractive vector mesons has been very successful at small t, the 
SCHNC in high energy small-t diffractive scattering is sill another major discovery at 
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HERA. On the theoretical side the chapter has not been closed - further studies of the 
impact of chiral-odd qq components in the photon suggested by D.Ivanov et al. [242] are in 
order. They were found to be crucial [261] for theoretical explanation of an approximate 
SCHC in the large-t photoproduction data. Here the BFKL based phenomenology of 
the absolute normalization and t-dependence is in a reasonably good shape, a fly in the 
ointment is a sign discrepancy in the helicity-fiip amplitude found by Poludniowski et 
al [261]. 

What new did we learn about the hard Pomeron trajectory? In inclusive DIS the pQCD 
Pomeron can only be probed at t = and one can not tell a difference between the 
fixed branching point and moving pole options for the Pomeron. In diffractive vector 
meson production the full t-dependence of the Pomeron exchange can be probed. The 
experimental observation by ZEUS CoUab. of Gribov's shrinkage of the diffraction cone 
in J/\l/ production is an important evidence for the hard Pomeron being a moving j- 
plane singularity — this is definitely the first important new finding on the hard Pomeron 
trajectory beyond the reach of inclusive DIS. One option is that the pQCD Pomeron is 
a sequence of isolated moving Regge poles as advocated in the pioneering Fadin-Kuraev- 
Lipatov work on the BFKL Pomeron [35] . Unfortunately, the accuracy of the combined set 
of the proton structure function and vector meson production data does not yet allow one 
to resolve the fine structure of those poles. The theoretical discussion has been confined 
to the evaluation of the slope of trajectories of these poles [25], the large-t behavior of 
trajectories remains an open issue. The color dipole approach predicts an approximately 
-independent shrinkage of the diffraction cone, there is a weak evidence for that from 
ZEUS measurements, but the experimental situation is not conclusive yet. 

Large t as a hard scale brings in new opportunities. In this case the pQCD Pomeron 
is expected to be in hard regime all the way through from the target to 7*1^ transition. 
The experimental information on large-t vector mesons is very exciting. The HI data on 
large-t J/\l/ mesons gave a very interesting evidence, supported also by ZEUS data, for 
the antishrinkage, a'p < 0, and emergence of the hard BFKL Pomeron exchange with 
intercept ap ~ 1.3 at \t\ > 3 GeV^. This is the second important new finding on the 
hard Pomeron trajectory beyond the reach of inclusive DIS. Such a transition from the 
shrinkage to the antishrinkage is plausible, but has not yet been explored theoretically. 
Here important issues for future theoretical studies are the sensitivity of the turn over 
from shrinkage to antishrinakge to the infrared regularization of pQCD and its (<5^ + ?7iy)- 
dependence. To this end, an experimental study of the interplay of two hard scales 
— (Q^ + '^y) ^iid \t\ — would be most interesting. The flavour and t-dependence of 
large-t cross sections does not conflict the estimates based on the leading order BFKL 
approach [261], but higher precision data are needed for more definitive conclusions. 

To summarize, the program of diffractive vector meson studies at HERA was exceptionally 
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fruitful one. The matching theoretical development followed, still more work is needed: 
the pressing issues include the t-dependence of the Pomeron trajectory from small to large 
t, understaning the role of higher Fock states in vector mesons and derivation of NLO 
/c_L-factorization, the further studies of helicity properties of large-t vector mesons. On 
the experimental side, new results on vector mesons are expected form several more years 
of run of HERA. 
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Note added 

After the main body of this review has been completed, D.Yu. Ivanov and his collaborators 
reported a long waited NLO analysis of 'j*V production [310-312]. It is a very involved 
calculation and not yet a full fledged NLO /c^-factorization analysis because the vector 
mesons have still been treated in the coUinear approximation. The major expectation 
was that the NLO calculations would fix more reliably the magnitude of the production 
amplitude, specifically, the NLO amplitudes must have a stability window as a function 
of the so-called factorization scale. To this end, the NLO results exhibit a discouraging 
instability of the pQCD expansion. First, the NLO corrections are twice as large in 
the magnitude, and of the opposite sign, than the LO amplitude of photoproduction of 
the J/\I^ [311]- Second, the NLO amplitude for the electroproduction of the p lacks an 
expected stability window vs. the factorization scale [312]. A further analysis of NLO 
correction, for instance, studies of the sensitivity of the stability window to models for 
skewed parton densities, and an independent rederivation are called upon to clarify this 
important issue. 
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